Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,141,415 T→G T75T (ACA→ACC)  unknown ← Minor pilin of type IV secretion complex, VirB5

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,141,4150TG100.0% 220.7 / NA 77T75T (ACA→ACC) unknownMinor pilin of type IV secretion complex, VirB5
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (40/37);  total (40/37)

AGCTCATACTTTCCACGGTGGAAACAATATTTGCTTCACAGGAAGTAATTTCACCTTGACGTAAATTCTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTAGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAGCCATTTATTTCACCATACAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCTGTCTGTGAGGCTAATTGGTCGGCATACGCTTTCAGCTCATTTTTATAATGCGAAACTGTATCCGTCCACTGTTTAGCCTCTTTTGCCCATTGTTGCACCGTTTGAAGCATGCTCATTTGACGTTGACCGGCAGCTGCGGCATCAACAACCGGAATACCGCTGGCTTGTACCAATACCGGCATAGATAAAACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAAC  >  1/1141174‑1141660
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
agCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCGCCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAg                                                                                                                                                                                                                                              >  1:325610/1‑251 (MQ=255)
agCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAg                                                                                                                                                                                                                                              >  1:603477/1‑251 (MQ=255)
agCTCATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAg                                                                                                                                                                                                                                              >  1:306939/1‑251 (MQ=255)
 gCTCATACTCTCCACGGTGGAAACAATATTTGCTTCGCAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAgg                                                                                                                                                                                                                                             >  1:32216/1‑251 (MQ=255)
    cATACTCTCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTa                                                                                                                                                                                                                                          >  1:495735/1‑251 (MQ=255)
        ctctCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTg                                                                                                                                                                                                                                      <  1:59247/248‑1 (MQ=255)
        ctctCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTg                                                                                                                                                                                                                                      <  1:350774/248‑1 (MQ=255)
          ctCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCcggtc                                                                                                                                                                                                                                                    >  1:830797/2‑235 (MQ=255)
          ctCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGt                                                                                                                                                                                                                                    >  1:728930/2‑251 (MQ=255)
          ctCCACGGTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGt                                                                                                                                                                                                                                    >  1:724952/2‑251 (MQ=255)
                 gTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTgg                                                                                                                                                                                                                                     >  1:816596/1‑243 (MQ=255)
                 gTGGAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATa                                                                                                                                                                                                                             >  1:854013/1‑251 (MQ=255)
                   ggAAACAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACg                                                                                                                                                                                                                           >  1:92052/1‑251 (MQ=255)
                     aaaCAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGcc                                                                                                                                                                                                                         <  1:234735/251‑2 (MQ=255)
                     aaaCAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGcc                                                                                                                                                                                                                         <  1:639592/251‑2 (MQ=255)
                     aaaCAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGc                                                                                                                                                                                                                          <  1:515667/250‑1 (MQ=255)
                     aaaCAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGc                                                                                                                                                                                                                          >  1:42943/1‑250 (MQ=255)
                        cAATATTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGccttg                                                                                                                                                                                                                      >  1:656982/1‑247 (MQ=255)
                          atatTTGCTTCACAAGAAGTGATTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGccttgag                                                                                                                                                                                                                    <  1:358342/251‑7 (MQ=255)
                                            gtgaTTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATgaga                                                                                                                                                                                                  >  1:142724/4‑248 (MQ=255)
                                              gaTTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATgagata                                                                                                                                                                                                >  1:239945/2‑246 (MQ=255)
                                              gaTTTCACCTTGACGCAAATTTTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATgagata                                                                                                                                                                                                >  1:350373/2‑246 (MQ=255)
                                                                            gTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTg                                                                                                                                                                  >  1:897788/1‑251 (MQ=255)
                                                                                     gTTGTATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCATTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgt                                                                                                                                                          <  1:492531/250‑1 (MQ=255)
                                                                                         tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                     <  1:253734/251‑4 (MQ=255)
                                                                                         tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                     <  1:808919/251‑4 (MQ=255)
                                                                                         tATTTTTCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                     <  1:422696/251‑4 (MQ=255)
                                                                                           tttttCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTGCATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaacag                                                                                                                                                   <  1:889162/251‑6 (MQ=255)
                                                                                            ttttCAAACAGGCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaacagt                                                                                                                                                  >  1:740101/1‑245 (MQ=255)
                                                                                                      ggCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGttttgctt                                                                                                                                                                        >  1:827483/1‑214 (MQ=255)
                                                                                                       gCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATg                                                                                                                                                                                                     <  1:393345/189‑1 (MQ=255)
                                                                                                       gCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                           <  1:85713/247‑5 (MQ=255)
                                                                                                       gCTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                           <  1:707896/247‑5 (MQ=255)
                                                                                                        cTTTTGGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacata                                                                                                                                        >  1:737084/1‑242 (MQ=255)
                                                                                                            tggCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatactgatt                                                                                                                                  <  1:889212/249‑14 (MQ=255)
                                                                                                               cATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGttt                                                                                                                                                >  1:341961/1‑234 (MQ=255)
                                                                                                                     cgcTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                           <  1:848010/233‑5 (MQ=255)
                                                                                                                     cgcTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGa                                                                                                                         >  1:703525/1‑251 (MQ=255)
                                                                                                                       cTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGcc                                                                                                                                                                                                                         >  1:49025/1‑152 (MQ=255)
                                                                                                                                               aCTAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTCGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCgg                                                                                               <  1:447852/251‑1 (MQ=255)
                                                                                                                                                 tAAATCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGaa                                                                                             <  1:803527/251‑1 (MQ=255)
                                                                                                                                                    aTCGATTACCGGGGTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATAc                                                                                          <  1:110855/251‑1 (MQ=255)
                                                                                                                                                              ggggTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGcc                                                                                                                                                                                                                         >  1:387467/1‑113 (MQ=255)
                                                                                                                                                               gggTAAACTCCGATTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTga                                                                               >  1:245680/1‑250 (MQ=255)
                                                                                                                                                                           aTTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATcaactac                                                                                                  >  1:620605/1‑217 (MQ=255)
                                                                                                                                                                           aTTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCggt                                                                   >  1:687568/1‑250 (MQ=255)
                                                                                                                                                                           aTTTTAAACCATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCggt                                                                   >  1:244693/1‑250 (MQ=255)
                                                                                                                                                                                  accATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGaca                                                            <  1:25655/250‑3 (MQ=255)
                                                                                                                                                                                   ccATTTATTTCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                           <  1:93089/171‑5 (MQ=255)
                                                                                                                                                                                       ttatttCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatactga                                                                                                                                    <  1:341936/174‑12 (MQ=255)
                                                                                                                                                                                         atttCACCATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATcaact                                                                                                    >  1:695142/1‑203 (MQ=255)
                                                                                                                                                                                               ccATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatac                                                                                                                                       >  1:707682/1‑155 (MQ=255)
                                                                                                                                                                                                cATATAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGcc                                                                                                                                                                                                                         >  1:405549/1‑79 (MQ=255)
                                                                                                                                                                                                cATATAAACTATGCGCTCCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGcagc                                                                                                              <  1:435093/187‑4 (MQ=255)
                                                                                                                                                                                                 atataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAg                                             <  1:490461/247‑1 (MQ=255)
                                                                                                                                                                                                 atataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAg                                             <  1:260477/247‑1 (MQ=255)
                                                                                                                                                                                                  tataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTgaactaatgc                                                                       <  1:809715/221‑10 (MQ=255)
                                                                                                                                                                                                  tataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTgaac                                                                             <  1:809530/215‑4 (MQ=255)
                                                                                                                                                                                                  tataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGt                                            >  1:735545/4‑251 (MQ=255)
                                                                                                                                                                                                   ataaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                           >  1:783842/3‑151 (MQ=255)
                                                                                                                                                                                                     aaaCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatac                                                                                                                                       >  1:463491/1‑149 (MQ=255)
                                                                                                                                                                                                       aCTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatac                                                                                                                                       <  1:305237/155‑9 (MQ=255)
                                                                                                                                                                                                            gcgcTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTa                                  >  1:655244/1‑251 (MQ=255)
                                                                                                                                                                                                               cTTCACCAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatactgatt                                                                                                                                  <  1:114637/152‑14 (MQ=255)
                                                                                                                                                                                                                    ccAATAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGAtttt                          >  1:828561/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaTAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCActgttgaac                                                                                                                                                     >  1:113527/1‑123 (MQ=255)
                                                                                                                                                                                                                      aaTAAAGGCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaac                                                                                                                                           >  1:296677/1‑132 (MQ=255)
                                                                                                                                                                                                                            ggCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTgtaacatac                                                                                                                                       >  1:66591/1‑126 (MQ=255)
                                                                                                                                                                                                                             gCTTGTACATCACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTGTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGt                  <  1:592353/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCtt           <  1:142079/247‑1 (MQ=255)
                                                                                                                                                                                                                                       cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac       <  1:725014/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac       <  1:691574/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cACGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATac       <  1:294045/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aCGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATaca      <  1:337627/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cGAATTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATacac     <  1:58906/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aaTTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGacaagac                                                        >  1:177883/1‑192 (MQ=255)
                                                                                                                                                                                                                                            aTTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACaa   >  1:342920/1‑250 (MQ=255)
                                                                                                                                                                                                                                            aTTCCGGTCTGTGAGGCTAATTGGTCGGCATACGCCTTGAGTTCGTTTTTGTAATGAGATACGGTTTCTGTCCATTGTTTTGCTTCTTTTGCCCACTGTTGAACAGTTTGTAACATACTGATTTGCCGTTGACCTGCAGCAGCACCATCAACTACCGGAATACCGCTAGCTTGAACTAATGCCGGTATAGACAAGACAGCAATTAAAGTGAGTTTTTTAAAGATTTTCATTACGTTCTCCTTATACACAAc  <  1:799367/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     
AGCTCATACTTTCCACGGTGGAAACAATATTTGCTTCACAGGAAGTAATTTCACCTTGACGTAAATTCTTACAACGGTCAAAAAGGTTGTATTTTTCAAACAGGCTTTTAGCATTCGCGCTTAATGCGTTTTTCCCACCGGAAACTAAATCGATTACCGGGGTAAACTCCGATTTTAAGCCATTTATTTCACCATACAAACTATGCGCTTCACCAATAAAGGCTTGTACATCACGAATTCCTGTCTGTGAGGCTAATTGGTCGGCATACGCTTTCAGCTCATTTTTATAATGCGAAACTGTATCCGTCCACTGTTTAGCCTCTTTTGCCCATTGTTGCACCGTTTGAAGCATGCTCATTTGACGTTGACCGGCAGCTGCGGCATCAACAACCGGAATACCGCTGGCTTGTACCAATACCGGCATAGATAAAACAGCAATTAAAGTGAGTTTTTTAAACATTTTCATCACGTTCTCCTTATACACAAC  >  1/1141174‑1141660

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑