Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,889,365 C→T K245K (AAG→AAA)  cysI ← sulfite reductase, hemoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,889,3650CT97.3% 204.6 / ‑1.6 73K245K (AAG→AAA) cysIsulfite reductase, hemoprotein subunit
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (31/40);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.68e-01

GACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGAC  >  NC_000913/2889124‑2889611
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
gACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGttt                                                                                                                                                                                                                                                     >  1:100880/1‑251 (MQ=255)
 aCGCTCCTCTTCCGCTTTCAACGTCTATACCCCCTCGCTCTCCAGCTGGTTTTTGGTTGTGGCATTTTTACGATCGGTTCTGTTACCCCAGTCTCGCTGAGTTGTCCCGA‑CGTC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCTCCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGtttt                                                                                                                                                                                                                                                    <  2:642335/251‑1 (MQ=255)
 aCGCTCCACTTCCCCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTCCGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGtttt                                                                                                                                                                                                                                                    <  1:867323/251‑1 (MQ=255)
  cgtcccACTGCCGCTTTAAACGTCCCAACCCCCACGCGTTCCAGCGTGTATTGGGGTTGGGCATTTTTACGATCGGTATGGTTCCCCCAGTTACGCTGATTTGTCGCGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTc                                                                                                                                                                                                                                                   <  1:857750/247‑1 (MQ=255)
   gcgccACTTCCGCTTTAAACGTCACACCCCCCACTCGCTCAAGCGTGTATTTGGTTTGGGCAGTTGTACGTTCGGGTCGGTTACCCCAGTCCCGTGGAGTTGTC‑CGAACGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCGCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCg                                                                                                                                                                                                                                                  <  1:283167/248‑1 (MQ=255)
      ccACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCg                                                                                                                                                                                                                                               >  1:165527/1‑251 (MQ=255)
           tCCGCTTTAAACGTCTCAACCCCCCCGCGCTCCGGCGTGTATTTGGTTGTGGCATTTTTACGGTCGGTTCGGTTCCCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCTTCAGATAGCCAAACTCGC‑TCGCGGGGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCg                                                                                                                                                                                                                                           <  1:858910/250‑1 (MQ=255)
           tCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAGGCCCCCCAGTTTGCCGTTTTCGGCGATCGc                                                                                                                                                                                                                                          >  1:724451/1‑251 (MQ=255)
               cTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCACCAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGGTAGCCACg                                                                                                                                                                                                                                      >  2:268677/1‑251 (MQ=255)
                    aaCGTCTAAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGGTGGTTTTCTTGTTGACGTGTTCGATGGAAAGCCCACGGCCCCCCAACAGGTTTAAGACCACAAGTTTGCCGTTTTCGGCGATCGCCACGa                                                                                                                                                                                                                                     >  2:180372/1‑247 (MQ=255)
                       gTCTCCACCCCCAGTCGCTCCAGCGTGTATTGGGTTTGGGTATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCTCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCAt                                                                                                                                                                                                                             <  1:569637/251‑1 (MQ=255)
                           cAACCCCCACGCGCTCCAGGGTGTATTTGGTTTTGCCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGACACGA‑AGTTCATGT‑Cg                                                                                                                                                                                                                        <  1:217643/251‑1 (MQ=255)
                                       gcTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAgg                                                                                                                                                                                                           >  1:723549/1‑250 (MQ=255)
                                          cccgCGTGTTTTTGTTTTTGGCATCTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGTCGTGTTCGATGGAAAGCCCTCCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCg                                                                                                                                                                                                        <  1:449517/248‑1 (MQ=255)
                                                   aTTTGGTCTTGGCTTTTTTACGATCGGTTCGGGTAACCCAGTCACGCTGAGTTGTTACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGACCAACTCAC‑TCGCGGTGCGGGCGGAGGTTTTCTTGTTGCCGAGTTCGATTGAAAGCCCACCGCCCACCACCAGGTTAAAGCCCAACAGTTTGCCGTTTTCGGCGATCGCCACGG‑AGTTAATGT‑CGTTGGCG‑TGcgtatagcaatcgttt                                                                                                                                                                                               >  2:368291/1‑236 (MQ=255)
                                                      tggtttggtctttttttCGATGGGTTCGGTTCCCCCATGCCCGCTGAGTTGTCCCGA‑CGGT‑GTCGGCCGCCGCCAGCGTAGGCTCCTGCGGCATATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCTCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCtg                                                                                                                                                                                             <  2:734702/239‑1 (MQ=255)
                                                       ggTTTTGTCTTTTTTTCGTTGGGTGCGGTTCCCCCAGTCCCGCTGAGTTGGCACG‑‑CGGCCTTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCtgtg                                                                                                                                                                                           <  2:81670/251‑1 (MQ=255)
                                                       ggTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGA‑TTCGGCCACCGCCAGCGAATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACAAACAGGTTAAAGACCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGATCA‑GTGCGTTGGCG‑TGccaataagtatcgttctt                                                                                                                                                                                             >  2:181066/1‑232 (MQ=255)
                                                        gTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCTTTCTGTgg                                                                                                                                                                                          >  2:35057/1‑251 (MQ=255)
                                                          tttGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCTCGA‑CGGC‑TTCGGACACCGCCAGCGTCTGCTCAAGCGGCAGTTAGCCAAACT‑ACATCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGACCACCGCCCACCAACAGGCTAAAGCCCACCAGTTTGCCGTTTACGGCGCTCGCCACGACAGT‑CATGT‑AGTAGGCGCTG‑ACATCGATATAGTTCTGTGGCg                                                                                                                                                                                        >  2:363806/1‑251 (MQ=255)
                                                          tttGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCAGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGC‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTgcggcg                                                                                                                                                                                        >  1:521246/1‑251 (MQ=255)
                                                             ggCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCgggg                                                                                                                                                                                     >  1:309822/1‑250 (MQ=255)
                                                                      aCGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑CGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCTTTTTGTGGCGGG‑agttcacccc                                                                                                                                                                           >  2:437248/1‑242 (MQ=255)
                                                                        gATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGC‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCCAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCCACAGGTTACAGCCCACCAGTTTGCCGTTTTCGGCGGTCGCCCCGA‑CGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑GTCACC‑ACCGTg                                                                                                                                                                        >  2:338220/1‑251 (MQ=255)
                                                                         aTCGGTTCGGTTACCCCAGCCCCGCTGAGTTGTCACGA‑CGGC‑TTCGGCCTCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTgg                                                                                                                                                                       <  1:810563/251‑1 (MQ=255)
                                                                          tcggttcggtTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCgggggtcatcccctggg                                                                                                                                                                        >  2:228784/1‑237 (MQ=255)
                                                                            gggtcggttACCCCAGTCACGTGGATTTGCCATGA‑GGGC‑TTCGGCACCCGCCAGCGGATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTGCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑ttt                                                                                                                                                                   <  2:612494/248‑1 (MQ=255)
                                                                              ttcggttACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTCGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACCGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGC‑AGGTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCggggagcaaaaccggagtttg                                                                                                                                                                    >  2:557371/1‑233 (MQ=255)
                                                                               tcggttACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑AACACttccctgggtttt                                                                                                                                                                   >  1:60563/1‑238 (MQ=255)
                                                                               ccggttACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCCGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCACACCAACCGGTTCAAGCCCACCCGTTTGCCGATTTCGGCGATCGCCACGC‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTTTGGCggaggcggaacctgggggtttg                                                                                                                                                                   >  2:368773/2‑231 (MQ=255)
                                                                                         cAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACAAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGCGG‑TTTTGTATTTATGCgg                                                                                                                                                      >  1:270900/1‑251 (MQ=255)
                                                                                              aCGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGAAGCTAGCCAAACTCAC‑TCACGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGACAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCAATCACCG‑GG‑GTTTGAATTTAcgccgaaggg                                                                                                                                                 >  1:780192/1‑244 (MQ=255)
                                                                                                    aGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑GTCAAT‑ACCCTGGGTTTTGA‑TTTTCGGGGCaggtggggct                                                                                                                                           >  1:468788/1‑245 (MQ=255)
                                                                                                       gtggtcaggccgtctccGGCCACCGCCAGCGTATGCGCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTg                                                                                                                                          <  1:673691/235‑1 (MQ=255)
                                                                                                           tccctccggtttcggcacccGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACGCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGcc                                                                                                                                       <  2:297305/232‑1 (MQ=255)
                                                                                                              cgacggttccgcccACCGCCAGCGTATGTTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTGCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCATCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑AGCACT‑ACCGTGG‑GGTTGAATTTACGCGGCAGGTAGGTCTGTCCGAgg                                                                                                                                   <  2:549445/239‑1 (MQ=255)
                                                                                                                cgGC‑TTCGGCCCCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑GCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATc                                                                                                                                <  2:165527/251‑1 (MQ=255)
                                                                                                                  gC‑TTCGGCAGCCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑CCCGTGG‑TTTTGAATTTCCGCGGCAGGTAGGTCTGGCCGAGGATCgg                                                                                                                              <  1:41356/251‑1 (MQ=255)
                                                                                                                   ggtgcGGCCACCGCCAGCGTATGCTCAAGCGGCAGATAGCCAACCTCAC‑GCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGTG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCgg                                                                                                                              <  2:276765/247‑1 (MQ=255)
                                                                                                                     ttCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCAAT‑ACCGTGG‑TTTTGAATTTAAGCGGGCGGTTGGGCTGGCCGGGGAtggggg                                                                                                                            >  1:577264/1‑246 (MQ=255)
                                                                                                                       cGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑AACGGGG‑GTTTGGATTTACGCGGCAGGGAGGGCTGGGCGAGGGTCgggttt                                                                                                                          >  1:870869/1‑247 (MQ=255)
                                                                                                                         gacgcgcgcggctttgtcTCCAACGGCTGATAGCCCAACGCAC‑TTGCGGTGCTGGCGTAGGGTTTCTTGTTGCCGTGCGCGATGGAGAGCCCACCGCCCGCAAACAGGTTAAAGCCCACCCGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCATGTAGGTCTGGCCGAGGATCGGttcttc                                                                                                                        <  2:595214/234‑1 (MQ=255)
                                                                                                                          cccccGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTtca                                                                                                                       <  2:734652/248‑1 (MQ=255)
                                                                                                                           ccccgacAGCGTATGCTCCAGCGGCCGATAGCCAAACTCCC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCCCCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTtcat                                                                                                                      <  1:446851/245‑1 (MQ=255)
                                                                                                                              cGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGGGG‑TTTTGAATTTACGCGGGCGGGAGGGATGGCCGGGGAGCGGTTCTtcttca                                                                                                                    >  1:591177/1‑250 (MQ=255)
                                                                                                                                 cAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑GTTACT‑ACCGTGG‑GTTTTAATTTTTGCGGCAGGGAGGGCTGGCCgggggtagggtcttctttcgt                                                                                                                  >  1:443166/1‑229 (MQ=255)
                                                                                                                                         gCTCCAGCGGCAGATATCCAAACGCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCag                                                                                                                   <  2:599951/240‑1 (MQ=255)
                                                                                                                                         gCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑GTTTTAATTTACGCGGCAGGGAGGGCTGGCCGCGGATCGGTTCTTTATCag                                                                                                                   >  1:599951/1‑240 (MQ=255)
                                                                                                                                            ccAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGGGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACttt                                                                                                      <  2:121693/250‑1 (MQ=255)
                                                                                                                                               gCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTcc                                                                                                  <  1:228784/251‑1 (MQ=255)
                                                                                                                                               gCGGCAGAGAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTcc                                                                                                  <  1:353713/251‑1 (MQ=255)
                                                                                                                                                       aaGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAg                                                                                          <  1:724296/250‑1 (MQ=255)
                                                                                                                                                         ggccaaacccaaccgcgGTGCGGGCGGTGGTTTGCTTGTTGCCGTGTTCGATGGAAAGGCCTCCGCCCACCATCAGGTTAAAGCCCACCAGTGTGCCGTTTTCGTCGATCGCCACTA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGc                                                                                         <  2:267499/238‑1 (MQ=255)
                                                                                                                                                         gCCAAACTCAC‑GCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGcc                                                                                        <  1:35057/251‑1 (MQ=255)
                                                                                                                                                           caacctcac‑tcGCGGTGGGGGCGTAGGTTTTATTGTTGACGTGTTCGATGGAAAGCACGCCGACCGCCAGCAGGTTGAAGCTCACCAGTTTGCCGTTATCGGCGAGCGCCGCGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑CCCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTCCCTGGTCGAGCCAg                                                                                      <  1:145042/247‑1 (MQ=255)
                                                                                                                                                              actcac‑tcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGAtc                                                                                   <  1:593214/251‑1 (MQ=255)
                                                                                                                                                                 cac‑tcGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTTAATTTTCGCGGCAGGTAGGGCTGGCCGAGGGTCGGTTCTTCATCTGTAGTGGCGACTTTTTTCTTGGCGAGCCAGATCTcc                                                                                >  1:433492/1‑251 (MQ=255)
                                                                                                                                                                  cc‑tcGCGGTGCGGGCGTAGGTTTGCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGTCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCg                                                                               <  1:557371/250‑1 (MQ=255)
                                                                                                                                                                      cgcgGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑GTCCAT‑ACCGTGG‑GTTTGAATTTATGCGGCCGGGAGGGCTGGGCGAGGGTAGGGTTTTCCACCGAAGTGGCGGCTTTTTTCTGGGCGGGCCACATCTCCGCac                                                                            >  1:255625/1‑250 (MQ=255)
                                                                                                                                                                      cgcgGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑CTCACT‑ACCGTGG‑TTTTTTATTTACGCGGCAGGGAGGGCTGGCCGCGGGTCGGGTCTTCATTAGTAGGGGGGGCATTTTCCTTGTTGGGCCAGATCTCCGCAt                                                                            >  1:72015/1‑251 (MQ=255)
                                                                                                                                                                                gCGTAGGTTTTCTTGTTTCCGTGTTCGATGCAAACCCCACCGCCCACCAACAGGTTACAGCCCACCAGTTTTCCGTTTCCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCTCATACGCGCGGGt                                                                  <  2:657461/251‑1 (MQ=255)
                                                                                                                                                                                   tatgTTTTCTTGTTGCCGTGTTCGATGGACAGCCCACCGCCCACCAACAGGTTAAAGCCCTCCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGGTCTGTGGCGGG‑AGCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGGTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTAc                                                                <  2:784314/247‑1 (MQ=255)
                                                                                                                                                                                                  ccGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGGGAT‑ACA‑CCCGGGG‑TTTTGAATTTTCGCGGGAGGGAGGGCTGGGCGGGGATCGGGTTTTCATCTGAAGGGGGGGCTTTTTTCTGGTCGAGCCCGATATCCGCCTACGCCCGGGGAGGAGGCcacaaaagcgc                                                >  2:265041/1‑244 (MQ=255)
                                                                                                                                                                                                      gtTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCCGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑GTCACT‑AACGTGG‑GTTTGGATTTACTCGGGAGGGAGGGCTGGCCGGGGGGCGGGTATTTATTAGGAGAGGGGGCTTTTTCCTGGTCGGGCCAGATCTCCGCATAAGGGCGGGGACGGGGACACAGATggccaga                                             >  1:718718/1‑244 (MQ=255)
                                                                                                                                                                                                         cGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGTTTGCCGTGTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGTCAGGTAGGTCTGGCCGAGTATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAt                                          <  1:406523/250‑1 (MQ=255)
                                                                                                                                                                                                             ggAAAGCCCCCCGCACACCTACAGGTTAAAGCCCACCAGTTTGCCGTGTGCGGCGATCGACACGA‑AGTTCTGGT‑CGTTGGCG‑TGCAGATCGATATCGGTCTGGGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGTCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCATCAGATGCTCAGAAAtcttc                                      <  2:100880/250‑1 (MQ=255)
                                                                                                                                                                                                             ggAAAGCCCACCGCCCACCAACAGGTTAAAGCCCCCCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCTCT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTTGTGGCGACTTTTTCCTGTTCGAGCCAGATTTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAAtcttc                                      <  2:309822/250‑1 (MQ=255)
                                                                                                                                                                                                                    ccaccgcccaccAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTTAATTTTCGCGGCAGGTAGGTCTGGCCGGGGATCGGGTCTTTATCAGGAGGGGCGACTTTTTTCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCCGAAATCTTCTTCGcccc                              >  1:348374/1‑250 (MQ=255)
                                                                                                                                                                                                                         gcccaccAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGt                         <  1:237734/251‑1 (MQ=255)
                                                                                                                                                                                                                         gcccaccAACAGGTTAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGt                         <  1:217702/251‑1 (MQ=255)
                                                                                                                                                                                                                                       tAAAGCCCACCAGTTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTg            <  1:814668/250‑1 (MQ=255)
                                                                                                                                                                                                                                           gCCCACCAGTTTGCCTGTTGCGGCGATCCCCACGA‑TGTTCTTGT‑CGATGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAAGTTTCGCGGCAGGTAGTTCGGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCATCAGATGCGCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCt       <  2:333631/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ccaGTTTGCCGTTGTCTGCGATCTCCACTA‑TGTTCATGT‑CGTTGGCG‑TGCCGCCCGGTATCGCTCTGTGTCGGG‑ATCACA‑ACCGTGG‑TTTTGAATTTGCGCGGCAGGTAGGTCTGGCCGAGGATCGTTTCTTCGTCAGTAGAGGCGACTTTTTCCTGGGCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCGCAGAAATCTTCCTCGCCCACTCGTACGCTTCCGCGTGCAgttgtgac  <  2:206561/251‑8 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                            
GACGCTCCACTTCCGCTTTAAACGTCTCAACCCCCACGCGCTCCAGCGTGTATTTGGTTTTGGCATTTTTACGATCGGTTCGGTTACCCCAGTCACGCTGAGTTGTCACGA‑CGGC‑TTCGGCCACCGCCAGCGTATGCTCCAGCGGCAGATAGCCAAACTCAC‑TCGCGGTGCGGGCGTAGGTTTTCTTGTTGCCGTGTTCGATGGAAAGCCCACCGCCCACCAACAGGTTAAAGCCCACCAGCTTGCCGTTTTCGGCGATCGCCACGA‑AGTTCATGT‑CGTTGGCG‑TGCAGATCGATATCGTTCTGTGGCGGG‑ATCACT‑ACCGTGG‑TTTTGAATTTACGCGGCAGGTAGGTCTGGCCGAGGATCGGTTCTTCATCAGTAGTGGCGACTTTTTCCTGGTCGAGCCAGATCTCCGCATACGCGCGGGTACGAGGCAACAGATGCTCAGAAATCTTCTTCGCCCACTCGTACGCTTCCGCGTGCAGCTGCGAC  >  NC_000913/2889124‑2889611

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