Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 74,046 A→G 18.4% L159P (CTC→CCC)  thiP ← thiamine/thiamine pyrophosphate ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091374,0460AG18.4% 188.1 / 4.3 81L159P (CTC→CCC) thiPthiamine/thiamine pyrophosphate ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (16/46);  new base G (14/2);  total (33/48)
Fisher's exact test for biased strand distribution p-value = 9.98e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCG  >  NC_000913/73810‑74283
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
aGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCGAGTTGGCGCAGGTCGCCG‑GGGATGTTTTCCAGGGCCTGGGGt                                                                                                                                                                                                                                                >  1:454307/1‑239 (MQ=255)
          cAGCGCCTGAAAGATTGCCAGCGCGAAAGTGGTCGCCTGCGGACCGCCCCGCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGGGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGc                                                                                                                                                                                                                                      <  1:68699/239‑1 (MQ=255)
           aGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTCACCACCGCCAATCGCCGACACCGAAAAAATGCAAGCTACGCCTCCCAACCTAGGCGGCAAGTTGACGCTGTTCGCCGCGGGAGTTTTTCCAG‑GCCTGGGGgaataaaaggc                                                                                                                                                                                                                                      >  2:541477/1‑233 (MQ=255)
           aGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGACGGGATTTGTCGCCGCAACCACGGCCATTCGACGAAGCGGAAAAAATACCAGCTCCGCATCCCAAGCAGGGCGGCAAGGTGACGCCATTCGCCG‑GGGGTGGTTTCCAGGGCAGGGAGtacaaaggggat                                                                                                                                                                                                                                     >  2:68699/1‑229 (MQ=255)
                   aTAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACCCGGTTCCCCG‑GGGAGGTTTTCCAGTGCCTGGGGTAATAAACGGCCCCCCATCg                                                                                                                                                                                                                             >  1:336906/1‑239 (MQ=255)
                     aGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCgg                                                                                                                                                                                                                            <  2:919930/238‑1 (MQ=255)
                      gATTGCCAGCTCGATAGGGGTCCCCGGGCGACCGCCCCCCAGCGATAGCCCGGTGGCGAAGCTGGCAAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCa                                                                                                                                                                                                                          <  1:52812/239‑1 (MQ=255)
                      agggacaggCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCa                                                                                                                                                                                                                          <  1:841594/231‑1 (MQ=255)
                         tGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTTTCGCCGCAACCACGGCCATTCGACGAAGCAGAAAAAATGCCAGCTACGCATCCCCAGCCGGGCGGCAAAGTGGAGGTGGTGGCCG‑GGGATGTTTTCCAGGCGCGGGAGtaaaaaaaaaaggcgcaacggcaaat                                                                                                                                                                                                                       >  2:90226/1‑219 (MQ=255)
                              gCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCCTAACCACGGCCATTCGCCGAAGCGGCAAAAATGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTGTTTCCCC‑GGGGTGTTTTCCCGTGCCCGGGGAAAAAAGCGGCGGCCCACCGGCA‑GAGGaaaa                                                                                                                                                                                                                 >  2:52812/1‑239 (MQ=255)
                              gCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGGCGCCGTAACAACGGCCCTTCGACGAAGCGGACAAAATGCCAACTACGCATCCCCAGCTAGGCGGGCAGTTGACGCTGTTCGCCG‑GGGGTGTTTTCCCATGCCTGGGGCAGTAAGCGGCCCGCCaacagccggataaaa                                                                                                                                                                                                                  >  2:742919/1‑225 (MQ=255)
                              gCTCGATAGGGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGGGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTaaaa                                                                                                                                                                                                                 <  1:1085449/239‑1 (MQ=255)
                               cTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCCAGTTGACGCTGTTCGCCG‑GGGAGGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCA‑GATTaaaaa                                                                                                                                                                                                                >  1:179273/1‑239 (MQ=255)
                                     tggggcgcccTGCGGGCCCCCCCCCCGCGATGGCCCGGTGGCGAAGCCGGGTAAAAGGAGCAAAAAGTATAGCGCAGCAACCGGCGGGATTTGGCGCCGTAACCCCGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                           <  2:142454/230‑1 (MQ=255)
                                     aggggccggcTGCGGCCCCCCCCCCAGCGATAGCACGGTGGCGAAGGTGGCGAAACAGAGCATAAAGATAAGCCCAGCAACCGGCGGGATTTGTCGCCGTAACCCCGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                           <  1:1092837/229‑1 (MQ=255)
                                     aGTGGCCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                           <  1:406512/238‑1 (MQ=255)
                                         gtccccgcgggACCGCCCCCCGGCAGTAGCCCAGTGGCGAAGTTGGCGAACCAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGTCCATTCGACGAAGCGGAAAAAATGCCCGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                      <  2:372471/230‑1 (MQ=255)
                                         gTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGGTTTCCAGTGCCTGGTGTCATAAGCGGCCCGCCATCGGCC‑GCTTTAAAAAaacactggc                                                                                                                                                                                                      >  1:435122/1‑234 (MQ=255)
                                         gTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGGAACCACGGCCATTCGACGAAGCGGAAAAAAAGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTGTTTCCCG‑GGGAGGTTTTCCAGTGCCTGGGGAGAAAAACGGCCCCCCACCGGCA‑GGTTAAAAAaaacgtctct                                                                                                                                                                                                      >  2:841594/1‑233 (MQ=255)
                                          tCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                      <  1:493083/238‑1 (MQ=255)
                                          gcGCCGGCGGCCCGCCCCCCAGCGATAACACGGTGGCGAAGGCGGTGAAACAGACGATAAAGAGAAGCGCAGCACACGGGGGGATTTGTCGCCGTAACCACGGCCATGTGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGTCGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                      <  2:611590/237‑1 (MQ=255)
                                                 ccgACCGCCCCCCAGGGATAGCCCGGTGGCGAAGCTGGGGAAACAGAGACAAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑aaaa                                                                                                                                                                                            <  2:59781/237‑1 (MQ=255)
                                                 cGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCCGCTACGCATCCCAAGCTGGGCGGCAAGTTGACCCTGGTCGCCG‑GGGAGGTTTTCCCGAGCATGGGGTAATAAGCGGCCCCCCCTCCGCAGGTTTAAAAAACACA‑GGGCCCG‑A‑aaaa                                                                                                                                                                                            >  2:1035525/1‑239 (MQ=255)
                                                   gACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAACGGAAAAAATGCCAGCTAGGCATCCCAAGCTGGGCGGCAACTTGACGCTGTTGGCCG‑GGGAGGTTTTCCCGGGCCTGGGGTAAAAAGCGCGTCGCCACCGCCA‑GATTaaaaaaacaacaggacggacaaaaa                                                                                                                                                                                            >  2:284452/1‑220 (MQ=255)
                                                      cGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGCTGTTTTCCAGTGCCTGGAGGAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCt                                                                                                                                                                                       >  1:967334/1‑239 (MQ=255)
                                                         ccccgcgCGAGAGAACGGTGGCAAAGTGGGCTACACATAGCATAAAGAGACGCGCAGCACCGGGCGGGATTTGCCGCCGTAACCCGGGCCAGTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCGGGGGGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATGCCTTGc                                                                                                                                                                                    <  2:937788/233‑1 (MQ=255)
                                                                   tAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCA‑GCTTAAAAAACACATGGGCCAG‑C‑AAAAAACCTTTCAGGCCGTAgg                                                                                                                                                                          >  1:238295/1‑239 (MQ=255)
                                                                    tccgcGGGGGGGAAGGTGGCGAACCAGAGCAAAAAGATAAGCGCAGCGCCGGGGGGGGTTTGGAGCCCTACCCACTGGCAATCGACAAAGCGGAAAAAATGGCAGCGACGCATCCCAAGCTGGGGGGCAAGTGGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAT‑C‑AAAATACCGTGCAGGCCGTAgg                                                                                                                                                                          <  2:224511/234‑1 (MQ=255)
                                                                     gCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCA‑GGTTAAAAAAAACACGGGCCAG‑C‑AAAAAACCTTGCAGGCCGTAaggg                                                                                                                                                                        >  1:522274/1‑238 (MQ=255)
                                                                          gTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGACTTGTCGCCGCCACCACGGCCATTCGACGACGAGGAAAAAATGCCCGCTACGCATCACAACCTGGTCCGCAAAATGATGCTGTTCGCCG‑GGGCTGGTTTCCAGAGACTGGCGTAATAAGCGGCCCCCCACCCGCA‑GGTTAATAAACAAATGGCCCCC‑A‑AAAAAACCTTCCaagagccagggggcaaa                                                                                                                                                                   >  2:547671/1‑223 (MQ=255)
                                                                              ggaggCTGGGAAAACAGGGCATAAAGGAAAGCGCAGGAACCGGGGGGTTTTGTCGGCGTAACCCCGGCCATTCGACGAAGGGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGGCGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  2:546229/234‑1 (MQ=255)
                                                                              ggAAGCTGGCGAAACAGGGAAAAAAGATAAGCGAAGAACCCGGCGGGATTTGTCGCCGTAACCCCGGCCATGCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  2:1060563/237‑1 (MQ=255)
                                                                              ctaaGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:192428/236‑1 (MQ=255)
                                                                              cGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGTATTTGTCGCCTTAACCACGGCCTTTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:90226/238‑1 (MQ=255)
                                                                              cGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGt                                                                                                                                                                <  1:819234/238‑1 (MQ=255)
                                                                                aaGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAAGGCCAGCTACGCATCCCCAGCTGGGCGGCAAGTTGACGCTTTTCGCCG‑GGGATGTTTTCCCGGGCCTGGCGGGATAAGCGGCCCGCCATCGGGC‑GATTAAAAAAAACATCGGCGAG‑C‑AAAATACCTTGCAGGCCGAGGGGGGCAAAgaaca                                                                                                                                                              >  2:192428/1‑234 (MQ=255)
                                                                                 aGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCACGCAGGGCGGCCAGTTGACGCCGTTCGCCG‑GGGATGTTTTTCAGTGCCTGGGGTAATAACCGGCTCGCCACCGGCC‑GGTTAAAAAACACCAGGGCCAGCC‑AAAAAACATTCCGCGCC‑TAGGGGGAAAAggcacac                                                                                                                                                            >  2:1038925/1‑234 (MQ=255)
                                                                                     ggCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTTCAGGCCGGAAGGCGAAAAAGTCCACCCCa                                                                                                                                                        >  1:1056249/1‑239 (MQ=255)
                                                                                           aCAGAGCATAAAGATAAGCGAAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGa                                                                                                                                                  <  1:638188/239‑1 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCC‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGGGAAAAGGTCCACTCCAGAccagccgat                                                                                                                                             >  1:305805/1‑232 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                           >  1:576557/1‑179 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                           <  2:576557/179‑1 (MQ=255)
                                                                                                gCATAAAGATAAGCGCAGAAGCCGGCGGGATTGGTCGCCGTACCCACGGCCATTCGACGAAGCGGAAAAAGAGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGc                                                                                                                                                <  1:1036678/236‑1 (MQ=255)
                                                                                                    aaagaatagcggcggcaCCGGCGGGATTTGGCCGCGGCACCACCGCCATTCGACGAAGGGGAAAAAAGGCCAGCGACGCATCCCAAGCTGGGGCGCAAGTTGACGGTTTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGc                                                                                                                                         <  2:238295/223‑1 (MQ=255)
                                                                                                    aaaGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGc                                                                                                                                         <  1:1004320/239‑1 (MQ=255)
                                                                                                           aGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAAAgggt                                                                                                                                  >  1:606889/1‑239 (MQ=255)
                                                                                                           aGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGGAATAAGCGGCTCGCCACCGGCA‑GATTAAAAAACACATGGGCCCG‑CAAAAAAACCTTCCA‑GCCGGAGGGCGAAAAAGGCCCCTcccacaccgaccatttgccaggggg                                                                                                                                  >  2:1004320/1‑216 (MQ=255)
                                                                                                                    aCCGGGGGGATTTGGCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTgccagcca                                                                                                                         <  1:284452/239‑1 (MQ=255)
                                                                                                                         cggGTTTTCTCGCCGTAACCACGGCCATTCGCCGAAGCGGAAAAAATGCCAGCTGCGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCt                                                                                                                    <  2:1106928/239‑1 (MQ=255)
                                                                                                                          gggATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGAAGGGCGAAAAGGTCCACCCCAAACCGAGCGATTGGCAGAGTGTTGCCAGCCAGcccgg                                                                                                                   >  1:745890/1‑237 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCg                                                                                                                <  1:217155/238‑1 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCg                                                                                                                <  1:435599/238‑1 (MQ=255)
                                                                                                                              tttGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGGTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGGGTAATAAGCGGCTCGCCATCGGCC‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGGACACTCCAAACCACACGACTGGCAGAGGGTTTCCACCCAACCCCgggc                                                                                                                >  2:45298/1‑236 (MQ=255)
                                                                                                                                   cGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATa                                                                                                          <  1:630110/239‑1 (MQ=255)
                                                                                                                                   cGCCCTAACCCCGGCCTTTCGGCGAGGCGGAAAAAATGCCAGCTACGCTTCCCAAGCGGGGGGGCGAGTGGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACCCATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATa                                                                                                          <  2:454307/239‑1 (MQ=255)
                                                                                                                                                gCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                      >  2:302483/1‑136 (MQ=255)
                                                                                                                                                gCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGc                                                                                                                                                                                                      <  1:302483/136‑1 (MQ=255)
                                                                                                                                                    ttCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGa                                                                                         >  1:27558/1‑239 (MQ=255)
                                                                                                                                                     tCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTTCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAACCATCGGCCAG‑C‑AAAATACCTTGACGGCCGCAGGGCGAAAAGGGCCACTCCAGACCGAGCGATTGGACGAGGGTTGTCCCCCCGCCCTGGCGACCACAGACACCAAGaagaccaaa                                                                                        >  2:585493/1‑232 (MQ=255)
                                                                                                                                                      cgacgaAGCGAAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGaaa                                                                                       >  1:776337/1‑239 (MQ=255)
                                                                                                                                                         cgaAGGGGAAAAAATGCCGGCTACGATCCCCAAGGGGGGGGGAAAGTTTGCGGTGTTCGCCG‑GGGAAGTTTACCAGTGCCTGGAGTAATAAGCGGCTCGCCAGCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACa                                                                                    <  2:435122/239‑1 (MQ=255)
                                                                                                                                                          agatggttataagagacAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGGTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACa                                                                                    <  1:236545/222‑1 (MQ=255)
                                                                                                                                                               ggAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGAc                                                                               <  1:45298/238‑1 (MQ=255)
                                                                                                                                                               ggAAAAAAGGCCAGCTACCAATCCCAAGTGGGGCGGGAAGTTGACGGTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGAc                                                                               <  2:606889/238‑1 (MQ=255)
                                                                                                                                                                  aaaaaGGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGAGGATGTTCGCCG‑GGGATGTTTTCCAGTGCCGGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCaa                                                                            <  1:1038925/238‑1 (MQ=255)
                                                                                                                                                                   aaatggccagttcccgttCCCAAGGTTGGCGGAAAGTTGACGCTGTTCGCCG‑GGGAGGTTGCCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAaca                                                                          <  2:27558/222‑1 (MQ=255)
                                                                                                                                                                        gCCAGTTACGCATCCCAAGTGGGGCGGAAAGTTGACGTGGTTCGCCG‑GGGATGTTTTCCAGGGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCggg                                                                     <  2:776337/239‑1 (MQ=255)
                                                                                                                                                                          ctggtACGCATCCCAAGCTGGGCGGCAAGTGGACGCTGTTGGCCG‑GGAATGCTTTCCAGTGCCTGGATTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAGACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGgga                                                                    <  1:1064436/234‑1 (MQ=255)
                                                                                                                                                                          cAGCTCCGCATCCAAAACTGGGTGTCAAGGTGACCCTGTTCGCCG‑GGGATGTTTTCCAGTGACTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACTTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGgga                                                                    <  1:913933/238‑1 (MQ=255)
                                                                                                                                                                              tACGCCTCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATTCCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGAt                                                                <  1:1103626/238‑1 (MQ=255)
                                                                                                                                                                                            tGGGGGGCAAGTGAGCGCTGTTCGCCG‑GGGATGTTTTCACGGGCCGGGGAGAATAAGGGGCTCGCCATCGGGA‑GGTTAAAAAACAAATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGcac                                                  <  2:336906/238‑1 (MQ=255)
                                                                                                                                                                                              gggggcAAGTGGGCGCTGTTCGCCG‑GGGTTTTTTTACAGTGCCTGGAGTAAGAACCTGCCCGCCAGCGGCA‑GATTAAAAAACCCATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGcacac                                                <  2:305805/235‑1 (MQ=255)
                                                                                                                                                                                                                ttCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                           <  2:905321/67‑1 (MQ=255)
                                                                                                                                                                                                                ttCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAAcaca                                                                                                                                                                                                           >  1:905321/1‑67 (MQ=255)
                                                                                                                                                                                                                ttCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGc                                                                                                                 >  1:497201/1‑155 (MQ=255)
                                                                                                                                                                                                                ttCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGc                                                                                                                 <  2:497201/155‑1 (MQ=255)
                                                                                                                                                                                                                         ggATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGGCCAACACCGGGAGGTACCAGGTCATTGACCAAAAACGCaaaaagcccaccccaaccgga                     >  2:571897/1‑220 (MQ=255)
                                                                                                                                                                                                                               tttCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAgcgc               <  1:895024/239‑1 (MQ=255)
                                                                                                                                                                                                                               tttCCAGTGCCGGGGGTAAAAAACGGCTCGCCATGGGAA‑AATTAAAAAACCCAGGGGCCAG‑C‑AAAATTCCTGGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAgcgc               <  2:522274/239‑1 (MQ=255)
                                                                                                                                                                                                                                ttcaggggCCTGGAGTAAAAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGGGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGcc              <  1:585493/232‑1 (MQ=255)
                                                                                                                                                                                                                                           ggAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAAACGCGACCAACACCGGGAGGAGCAAGGTCACTGCACCCAAACCAAAAAGCGCCACCCGCCCCCGGAAGCGCCTGCGaacaagc   >  2:304624/1‑233 (MQ=255)
                                                                                                                                                                                                                                             aGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATagaacc  >  1:725313/1‑235 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
AGTCGTAACTCAGCGCCTGATAGATTGCCAGCTCGATAGTGGTCGCCTGCGGACCGCCCCCCAGCGATAGCACGGTGGCGAAGCTGGCGAAACAGAGCATAAAGATAAGCGCAGCAACCGGCGGGATTTGTCGCCGTAACCACGGCCATTCGACGAAGCGGAAAAAATGCCAGCTACGCATCCCAAGCTGGGCGGCAAGTTGACGCTGTTCGCCG‑GGGATGTTTTCCAGTGCCTGGAGTAATAAGCGGCTCGCCATCGGCA‑GATTAAAAAACACATGGGCCAG‑C‑AAAATACCTTGCAGGCCGTAGGGCGAAAAGGTCCACTCCAGACCGAGCGATTGGCAGAGTGTTGCCAGCCAGCCCTGGCGACCATAGACGCTAAGAATGCCGAAAACAGCGACCAACACCGGGAGGATCAAGGTCATTGCACACAGACGCAACAGCGCCAGCCGACCCGGAAAGCGCCTGCGATAGAGCG  >  NC_000913/73810‑74283

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