Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,084,903 T→G N307T (AAC→ACC)  fdoG ← formate dehydrogenase O subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,084,9030TG97.4% 391.2 / ‑0.3 114N307T (AAC→ACC) fdoGformate dehydrogenase O subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base G (52/59);  total (52/62)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

Only 100 of 114 total aligned reads displayed.

AGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACG  >  NC_000913/4084660‑4085148
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           
aGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGTGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGtat                                                                                                                                                                                                                                                       >  1:704448/1‑250 (MQ=255)
 ggAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTATCAGCAAGTTCCACACGCTGCGCGGATGTTGCAGCGTGGTATCGCGGTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAg                                                                                                                                                                                                                                                     <  2:711123/251‑1 (MQ=255)
      gcggcTTTTGGCGTACCACAGATGTTGTCAACCACGACTGGCGGGTTGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTGCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTc                                                                                                                                                                                                                                                <  2:772447/247‑1 (MQ=255)
      gcgTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTTGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTc                                                                                                                                                                                                                                                <  1:762009/251‑1 (MQ=255)
      gcgTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTc                                                                                                                                                                                                                                                >  1:747399/1‑251 (MQ=255)
        gTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCgg                                                                                                                                                                                                                                              >  1:37520/1‑251 (MQ=255)
            tttGGCGTACCTCCTATGTTTTCAACCCCATCGGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCATCGCGGATGTGGCAGCGTGGTATCGCGTTTGGCGAAACCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTATGTGCTTTTTTCCGCGTCGTAGCCGGTGATCAGGCCATCTTCAAAGCCGTACTCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGgtgt                                                                                                                                                                                                                                          <  2:716385/250‑1 (MQ=255)
            tttGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGgtgt                                                                                                                                                                                                                                          >  1:440878/1‑250 (MQ=255)
              tGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTAt                                                                                                                                                                                                                                        <  1:339429/250‑1 (MQ=255)
               ggCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGGGCTGTTTCAGCAAGTTCCAAACGCAGCGCGGATGTTGCAGCGTGGTATAGCGTTTGGCGAAGCCGTTTTCGTACAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAAAAGGCCATCTTCAAAGCCGTACTCCTCACGCCCGATCAGGCTGGCGGTGGGCTAGGCTTCGGTGAAttt                                                                                                                                                                                                                                      >  2:851702/1‑250 (MQ=255)
               ggCGTACCACAGATGGTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTc                                                                                                                                                                                                                                     <  2:531180/251‑1 (MQ=255)
                gCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTACACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTcg                                                                                                                                                                                                                                    <  1:259807/251‑1 (MQ=255)
                     ccACAGATTGTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTGCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGGCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGtt                                                                                                                                                                                                                               <  2:187805/251‑1 (MQ=255)
                       acaGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTACACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGTTgg                                                                                                                                                                                                                             >  1:426107/1‑250 (MQ=255)
                         aGATGTTTACAACACCATCGGGCGTGTCGCGGGAATCGTGCTGTGTCAGCAAGTTCTACACGCAGCGTGGATGGTGCAGCGGGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCTTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGTTGTATTCGCGGTTGAAt                                                                                                                                                                                                                           <  2:414893/251‑1 (MQ=255)
                           aTGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCTCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGGAttt                                                                                                                                                                                                                         >  1:156267/1‑251 (MQ=255)
                                ttCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGAGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCAt                                                                                                                                                                                                                    >  1:781483/1‑251 (MQ=255)
                                  cAACCCCATCTGGCTTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATTTTGCAGCGTGGTAGCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTg                                                                                                                                                                                                                  <  2:782239/251‑1 (MQ=255)
                                   aaCCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAATTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGAGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTCTTACGGGGTGAATTTTTTATTGt                                                                                                                                                                                                                 >  2:155612/1‑251 (MQ=255)
                                    aCCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGtt                                                                                                                                                                                                                <  1:262101/251‑1 (MQ=255)
                                      cacaGCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGTAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGGCCAGGAGGATTTATCGTACTTGCGTTTTGCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTca                                                                                                                                                                                                              <  1:742463/251‑1 (MQ=255)
                                        caTCTGGCGTGTGTCGGGAAACGTGCTTTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGATGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCTTCTTTAAAGCCGGAATCCTCACGCGCGATCAGTCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTcagc                                                                                                                                                                                                            <  2:702259/251‑1 (MQ=255)
                                           cgggCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGAGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTcagcag                                                                                                                                                                                                          <  2:641135/248‑1 (MQ=255)
                                             ggggtggaGCGGGAAACGGGCTGATGCAGCAAGTTCCACCAGCAGCGCGGATGTTGCAGCGTGGTATCGGGTTTGGCGAAGCCGTGTTCGGCCATGTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTAGAGGCTTCGGTGTAGTCGCGGTTGAATATTTCATTGTgcagca                                                                                                                                                                                                           <  2:576013/240‑1 (MQ=255)
                                             ggCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTGGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAgg                                                                                                                                                                                                         <  2:113801/249‑1 (MQ=255)
                                                gtgtAGCGGTAAACGTTCTGTTTCAGCAAGTTCCACACGCAGTGCGGATGTTGCATCGTGGTATTGCGTTTGGCGATGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACTGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCTGTGGTATAGGCTTCGGTGTATTCTCGGTTGAATTTTTCATTGTTCAGCAGGTACAg                                                                                                                                                                                                    <  2:185513/251‑1 (MQ=255)
                                                gtgtAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGc                                                                                                                                                                                                   >  1:175062/1‑251 (MQ=255)
                                                  gtAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCa                                                                                                                                                                                                  <  1:303446/251‑1 (MQ=255)
                                                    aGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGCAGACGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTGGAATTTTTCCTTGTTtagcaaggtcaggaa                                                                                                                                                                                                 >  2:830994/1‑240 (MQ=255)
                                                     gCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATa                                                                                                                                                                                               <  1:830994/251‑1 (MQ=255)
                                                        ggAAACGTGCTGTTTCAGCAAGTCCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGc                                                                                                                                                                                            <  1:738253/251‑1 (MQ=255)
                                                           aaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGACCAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGGATTCGCGGTTGAATTTTTCATTGTTTAGCAGGGAAAGCAATAcacccg                                                                                                                                                                                         >  2:736719/1‑246 (MQ=255)
                                                           aaCGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGATGGCGGTGGGTTAAGATTCGGGGTATTAGCGGTTGAATTTTTCATTGTTCAGCAGGATCAGCAATACACCTg                                                                                                                                                                                         >  2:10131/1‑251 (MQ=255)
                                                             cGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCCCGCTCAGGCTGGCGGCGGTCTAGGCTTCGGTGTATTCGCGGTTGGATTTTTCATTGTTCAGTACGACCAGCCATACGCCTGAc                                                                                                                                                                                       >  1:786482/1‑251 (MQ=255)
                                                                   gTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGa                                                                                                                                                                                  <  1:378242/249‑1 (MQ=255)
                                                                      tCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGc                                                                                                                                                                              >  1:536878/1‑250 (MQ=255)
                                                                       cAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCa                                                                                                                                                                             <  1:774437/250‑1 (MQ=255)
                                                                       cAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTCACGAAGCCGTTTTCGTCCACTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCAGTGGTCTAGGCTTCGGTGTATTCGCCGTTGGATTTTTCATTGTTCAGCAGGAACAGCAATACGCCTGACAGGAAAGCa                                                                                                                                                                             >  1:507765/1‑250 (MQ=255)
                                                                        aGCAAGTTACACACGCAGCGAGGATGTTGCAGCGTGGAATCGCGATGGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGCTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTACACGCCGTAATCCTCACGCAAGATCCAGCTGGCTGTGGTATGGGATTCGGTTAATTAGCGGGTGGCTTTTTCCTTGTTCAGCAGGTACataaagacaacagataggcaagaaat                                                                                                                                                                            >  2:694104/1‑226 (MQ=255)
                                                                             gTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCTTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAg                                                                                                                                                                      <  1:875230/251‑1 (MQ=255)
                                                                             gTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCATTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTTAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAg                                                                                                                                                                      <  2:463613/251‑1 (MQ=255)
                                                                             gTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAg                                                                                                                                                                      <  1:203618/251‑1 (MQ=255)
                                                                                ccACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATTGCAGTAc                                                                                                                                                                   >  1:607923/1‑251 (MQ=255)
                                                                                          cgcgGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAAt                                                                                                                                                         <  1:715011/251‑1 (MQ=255)
                                                                                           gcgGATGTTGCAGCGTGGTATCGCGTTTGGCGATGCCGTTTTCGTCCAGATCATAAGACCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCTCGCACGACCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCggta                                                                                                                                                              <  1:284010/245‑3 (MQ=255)
                                                                                            cgGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTACAGTTCATAAGTCAAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGCGGTAAAGGCTTTGGGGTATTCCCGGTGGTATTTTTCCTGTTTCAGCAGGTACAGCAATTAGCCTGGCAGGGACGCAACGTCAGTACCggcaccaaa                                                                                                                                                         >  2:236574/1‑242 (MQ=255)
                                                                                              gATGTTGCAGCGTGGTATCGCGTTTGGCGACGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATAAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTTATTGTTCAGCCGGTACAGCCATACGCCTGCACGGAAAGCAATGTCAGTACCGGAACGAATAggg                                                                                                                                                     >  2:819051/1‑251 (MQ=255)
                                                                                               aTGTTGCAGCTGGGTATCGCGTTTGGCGACTCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAgggg                                                                                                                                                    <  2:21370/251‑1 (MQ=255)
                                                                                                      aGCGTGGTATAGCGTTTGGCGAAGTCGTTTTCGTCCAGCTCATACGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTACTCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGTGATATTCGCGGTTGAACTTTTCCTTGTTCACCAAGGACCGGAATTCGCCTGCACGGAAAGCAATGTCAGTTCCGGAAAGAAGAGGGGCaacgg                                                                                                                                              >  2:401373/1‑246 (MQ=255)
                                                                                                        cGTGGTTTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATagtagt                                                                                                                                          <  2:37520/251‑1 (MQ=255)
                                                                                                        cGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGGTGAATTTTTTATTGTTCAGCAGGTACAGCAATACGCCTGACAGGGAAGCAAATTCCGTACCGGAACGAATAAGGGCAtagtagg                                                                                                                                          >  2:167881/1‑250 (MQ=255)
                                                                                                              agcgcgTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTGTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCt                                                                                                                                     <  2:156332/248‑1 (MQ=255)
                                                                                                              aTCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGGCCAGGAGGATTTAGCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCt                                                                                                                                     <  2:838381/250‑1 (MQ=255)
                                                                                                                cgcgTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAcc                                                                                                                                  <  1:835574/251‑1 (MQ=255)
                                                                                                                 gcgTTTGGCGAAGCCGTTTTCGGCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccg                                                                                                                                 <  1:722928/251‑1 (MQ=255)
                                                                                                                  cgTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgc                                                                                                                                <  1:307577/251‑1 (MQ=255)
                                                                                                                   gTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTAccgcc                                                                                                                               >  1:826112/1‑251 (MQ=255)
                                                                                                                      tGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGCTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCAATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGATGACTTTTTCATTGTTCTGCAGGGACAGCACTACGCCTGCCAGGAAAGCATTGTCAGTACCGGAACGAATAGGGGCATAAGTGGCAGCTAccgccgcc                                                                                                                            >  2:162336/1‑251 (MQ=255)
                                                                                                                       ggCGAAGCCGTTTTCGTCCAGTTCATAAGGCCAGGAGGATGTATCTTCCTTGCTTTTTTCCGCGTCGTATCCGGGGAACAGGCCATCTTCAAAGCCGTAATCTTCACGCACGACCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGTGGCATAGTAGTCAGCTAccgccgcc                                                                                                                            <  2:868815/250‑1 (MQ=255)
                                                                                                                          gAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGAATAGGCTTCGGTGTCTTCGCGGTTGAATTTTTCATTGTTCAGCAAGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCAAAGTAGTCAGCTACCGCCGCCGTAc                                                                                                                        >  1:652116/1‑251 (MQ=255)
                                                                                                                             gccgtgttctgccggTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAACGCCGTAATCCTCACGCACGAGCAGGCTGGCGGTGGTATAGGCTCTGGTGTAGTCGCGGTTGAATTTTTCATTGTTCAGCAGGTAGAGCAATACGCCTGACAGGACAGCAATGCCAGTGCCGGAACGAATAGGGGCATAGTAGTTAGCTACCTCCGCCGTAcgcg                                                                                                                     <  2:318866/237‑1 (MQ=255)
                                                                                                                              ccGTTTTCGTCCAGTTCAGAAGTCCTGGAGGATTTATCGTACTTGCGTGTTGCCGCGCTGTAGCCGGTGAAAAGGCCATCTTCAAAGCCGTATTCCTCACGCACGATCAGGCTGGGGGTGGTATAGGCTTCGGGGTTGTCGCGGTTGAATTTTTCATTGTTCAGCAGTTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTGCCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGt                                                                                                                    <  2:432369/251‑1 (MQ=255)
                                                                                                                                  tttCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTaa                                                                                                                  >  1:602426/1‑249 (MQ=255)
                                                                                                                                  tttCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAACGGAACGAATAGGGGCCTAGTAGTCAGCTACCGCCGCCGTAcgcggaa                                                                                                                  >  1:147557/1‑246 (MQ=255)
                                                                                                                                      gTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTAc                                                                                                                        <  1:529123/239‑1 (MQ=255)
                                                                                                                                      gTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCCGGTACAGCCATCCGCCTGCCAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGtcc                                                                                                                        >  2:529123/1‑237 (MQ=255)
                                                                                                                                        ccAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGga                                                                                                          <  2:612257/251‑1 (MQ=255)
                                                                                                                                             tCATAAGTCCAGGAGGATTGATCGTCCTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCTTAATCCTCACGCAGGATCAGGCTTGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGATAGCAATGTCAGTCCCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGgatcga                                                                                                      <  1:572018/250‑1 (MQ=255)
                                                                                                                                                      cAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATCGGCTTCGGGGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACCGCCATACGCCTGACAGGAAAGCAATGTCAGTACCGGGACGGATAGGGGAATAGTAGTCACCTACCGCCGCCGTACGCGTAAAGCGAggtgcgcgcccaacc                                                                                              >  2:219041/1‑236 (MQ=255)
                                                                                                                                                       aggaggATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGAGCGATCACAATCAGc                                                                                           >  1:623703/1‑251 (MQ=255)
                                                                                                                                                        ggaggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCt                                                                                          >  1:134427/1‑251 (MQ=255)
                                                                                                                                                         gaggaTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCAGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCtt                                                                                         <  1:226342/251‑1 (MQ=255)
                                                                                                                                                           ggaTTTATCGTACTTGCGTTTTTCCGCTTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcg                                                                                       <  1:370332/251‑1 (MQ=255)
                                                                                                                                                             aTTTATCGTACTTTCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgcg                                                                                     <  1:851702/251‑1 (MQ=255)
                                                                                                                                                             aTTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTcgcg                                                                                     >  1:502994/1‑251 (MQ=255)
                                                                                                                                                                        ttGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGACGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGGATTTTTCATTGTTCAGCAGGTCCAGCAATACGCCTGACAGGACAGCAATGTCAGTACCGGAACGAATAGGGGCCTAGTAGTCAGCTACCGCGGCCGGACGCGGAAAGCCAGGAACGATCTCACTCAGATTCGCGACGTtgagtat                                                                          >  1:682363/1‑246 (MQ=255)
                                                                                                                                                                          gcttttttCCGCGTCGTAGCCGGTTAACAGGCCTTCTTCAAAGCCTTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTTCCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAAtt                                                                         <  1:16468/247‑1 (MQ=255)
                                                                                                                                                                            gTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGGATAAGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCCGCAGGTACCGGAAAACGCCGGACAGGAAAACAATTTCAGGACCGGAACGAAGAAGGGCATAGTAAACAACTACCGCCTCCGTCCCGGCAAAGCGACGAGAGGTCACCATCAGCTTCGCCCCGTTGTGAGTTTTg                                                                      >  2:418023/1‑251 (MQ=255)
                                                                                                                                                                             tttttCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCAAGATCAGGCTGGCGGTGGTATAGGCTTCGGGGTATTCGCGGGTGAATTTTTCATTGTTCAGCAGGTACAGAAATACGCCTGAACGGACAGCAATATCAGTGCCGGAAAGACTAGGGGCCTAGTAGTCAGCTACCGCCGCCGTACGCGAAAAGCGAGGATCGCACACCATACGCTTCGCGCCGTCGTGGCTTTTgg                                                                     >  2:844277/1‑251 (MQ=255)
                                                                                                                                                                                gtcctcgcataagcgggTGAACAGGCCATCTTCAAAGCCGAAATCCACACGCACGATCAGGCTGGCGGTGGAATAGGCTTCTGTGTATTTGCGGTGGAATTTTTCATTGTGCAGCAGGCACAGCAATACGCCTGACAGGAAAGCAATGACAGTACCGGAAAGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTTCGCGTAAAGCGTGGATCGAGCACAATCAGCGTCGCGCCGTTGTTAATGTTGGCtt                                                                  <  2:121101/236‑1 (MQ=255)
                                                                                                                                                                                         gtAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTAAAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGAACCGGAACGAATAGGGGCCTAGTAGTCAACTACCGCCGCCGTACGCGTAAAACAGGGATCGATCACAAACAGCTTCGCGCCGTTGGGAATTTTGGCTTCCATCAccc                                                         >  1:802726/1‑251 (MQ=255)
                                                                                                                                                                                           aGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAg                                                       <  1:243546/251‑1 (MQ=255)
                                                                                                                                                                                               ggTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGa                                                   <  1:389756/251‑1 (MQ=255)
                                                                                                                                                                                                    aCAGGCCTTCCTCAAAGCCGTAATCCTCACGCACGATCAGGCGGGGGGTGGTATAGGTTTCGGTGTATTCGCGGTTAAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGGCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCg                                              <  2:858548/251‑1 (MQ=255)
                                                                                                                                                                                                        gCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATTTCAGTACCGGAACGAATAGGGGCCTAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCCAAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCCTCGCCCAGCGGAACCCCAccc                                          >  2:182284/1‑250 (MQ=255)
                                                                                                                                                                                                         ccATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCgg                                         <  1:868595/251‑1 (MQ=255)
                                                                                                                                                                                                           aTCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGAAAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTAAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACCATCAGCTTCGCGCCGCTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCgggg                                       >  1:179445/1‑250 (MQ=255)
                                                                                                                                                                                                               taaaaGCCGTAATCCTCACTCCCGATCAGGCGGGCGGGGGTATAGGCTTCGGTGTATTCGCGGTTGAAGTTTTCATTGTTCAGCAGGTACAGCAATGCGCCGGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTATTCAGCTCCCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATGTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGt                                       <  2:339664/245‑1 (MQ=255)
                                                                                                                                                                                                               tCAAAGCCGTAATCCTTACTCACGTATAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGAAAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAg                                   <  2:652116/250‑1 (MQ=255)
                                                                                                                                                                                                                cAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCt                                 >  1:184751/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaaGCCGTAAGCCTCACGCACGATCAGGCTGGCGGTGTTATAGGCTTCGGGTTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCtt                                <  1:28246/251‑1 (MQ=255)
                                                                                                                                                                                                                   aGCCGTAATCCTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTAATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCa                              <  2:175062/251‑1 (MQ=255)
                                                                                                                                                                                                                             cTCACGCACGATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTAACGGAACGAATAGGGGCATAGAAGTCAGCTTCCGCCGCCGTACGCGTAAAGCGGGGATTGATCAAACTCAGCTTCCCGCCGTTGTGAAATTTGGCTTCCCTCGCCCAGCGGCACCCGACCGGGGGCGCTTCCGAGGCGttt                    >  2:270123/1‑249 (MQ=255)
                                                                                                                                                                                                                                   cacgATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCa             <  1:736719/251‑1 (MQ=255)
                                                                                                                                                                                                                                    acgATCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTCCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCAt            <  1:730022/251‑1 (MQ=255)
                                                                                                                                                                                                                                        aTCAGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGAACAGCAAAACGCCTGACAGGAAAGCAATGTCCGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAACGAGGATCGATCACAATCCGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCACAACGGAACCCGACCGGGTGGGCTTCAGCGGCAGTTCCGCCAATcac         >  2:295286/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aTCAGGCTGGCGGTGGTATAAGCTTCGGGGTATTCGCGGTTTAATTTTTCATTGTTCAGCAGGTAAAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTATAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGAGAGCTTCAGCGGCGTTTCCGCCCATcac         >  2:428928/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aGGCTGGCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATcaccac      <  1:494314/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aGGCTGGCGGAGGTAAAGGCTTCGGTGTAGTGGCGGGTGAATATTTCGTTGTTCAGCAGGTACAGCAATACAACAGAAAGGAAGGAAATGTCAATACCGGGACGTAAAGGGGCATAGAAGTCAACTAACGACGCAGTACGAGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGGAATTTGGCATCCATCGCCCAAACGAACCCGCCCGGGGACGGTTCAGCCGCGGTTCCGCCCATcacacc      >  2:381640/1‑248 (MQ=255)
                                                                                                                                                                                                                                                ggCGGTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACAacgacg  >  1:860680/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                           
AGGAACGCGTCTTTTGGCGTACCACAGATGTTTTCAACCACATCTGGCGTGTAGCGGGAAACGTGCTGTTTCAGCAAGTTCCACACGCAGCGCGGATGTTGCAGCGTGGTATCGCGTTTGGCGAAGCCGTTTTCGTCCAGTTCATAAGTCCAGGAGGATTTATCGTACTTGCGTTTTTCCGCGTCGTAGCCGGTGAACAGGCCATCTTCAAAGCCGTAATCCTCACGCACGATCAGGCTGGCGTTGGTATAGGCTTCGGTGTATTCGCGGTTGAATTTTTCATTGTTCAGCAGGTACAGCAATACGCCTGACAGGAAAGCAATGTCAGTACCGGAACGAATAGGGGCATAGTAGTCAGCTACCGCCGCCGTACGCGTAAAGCGAGGATCGATCACAATCAGCTTCGCGCCGTTGTGAATTTTGGCTTCCATCGCCCAGCGGAACCCGACCGGGTGAGCTTCAGCGGCGTTACCGCCCATCACCACGACG  >  NC_000913/4084660‑4085148

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: