Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 3,878,107 T→C 15.6% G42G (GGA→GGG)  Y75_RS19225 ← RHS element protein RhsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077793,878,1070TC15.6% 383.5 / 8.5 167G42G (GGA→GGG) Y75_RS19225RHS element protein RhsA
Reads supporting (aligned to +/- strand):  ref base T (81/52);  new base C (2/29);  total (83/84)
Fisher's exact test for biased strand distribution p-value = 1.18e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 167 total aligned reads displayed.

CAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCAGAGTGAAATAATTAGGAAA  >  NC_007779/3877915‑3878339
                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
cAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATg                                                                                                                                                                                             >  2:555560/1‑238 (MQ=11)
cAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACCGCTGTAGGGGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGGCACTCCGCCGGGGCCCACCGAACAGGCCACGCCGGCGAGGGCACCAATg                                                                                                                                                                                             >  2:603216/1‑238 (MQ=2)
            ggtatCCGGAGGCTTTTTCCAGCCGGGGCCGAGGCTCCCCCCCGGGGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCt                                                                                                                                                                                 <  2:1432511/236‑1 (MQ=12)
            gATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCt                                                                                                                                                                                 <  2:89241/238‑1 (MQ=32)
            gATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCt                                                                                                                                                                                 <  1:106943/238‑1 (MQ=32)
            gATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCt                                                                                                                                                                                 <  2:475531/238‑1 (MQ=32)
            catttCCGCGGGCATTTTCCAGCCGGGGCGGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAGCGGCAGCGGGCCGGGCAGGGCGGTGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GGGCTGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCGCCAATGCGCACCCCGGCt                                                                                                                                                                                 <  2:1128882/234‑1 (MQ=11)
             atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGGACCAATGCGGAGCCCGGCTGa                                                                                                                                                                               >  2:25136/1‑239 (MQ=14)
             atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTg                                                                                                                                                                                <  1:555560/238‑1 (MQ=32)
             atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGa                                                                                                                                                                               >  1:1168649/1‑239 (MQ=32)
             atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGa                                                                                                                                                                               >  1:150977/1‑239 (MQ=32)
             atatCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATGGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGa                                                                                                                                                                               >  1:1699076/1‑239 (MQ=17)
                tccgagggCTTTTCCCAGCCGGGGCCGAGGCTCCCCACCGGGGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAAcc                                                                                                                                                                            <  2:1035019/233‑1 (MQ=12)
                tcccccggtaTTTTCCAGCCGGGGCCGGGGCTCCCCACCGGGGCGGGCGGTTTTGGCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAAcc                                                                                                                                                                            <  2:667704/230‑1 (MQ=255)
                tCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAAcc                                                                                                                                                                            <  1:1799757/239‑1 (MQ=32)
                tCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGGCGGTTTCACCGGGCAGGACCTTT‑GCACC‑GAGCAGGGGGTTGACC‑GGATGGCCGGACGGCACCCCGCCGGGGCACACCGAACCGGCCACGCCGGTGGGGGCCACACTGCGCACCCGCGCTGAAc                                                                                                                                                                             >  2:1497813/1‑238 (MQ=11)
                     aGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGa                                                                                                                                                                        <  1:1761835/238‑1 (MQ=32)
                       gCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACa                                                                                                                                                                     >  1:667704/1‑239 (MQ=32)
                                 gCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTg                                                                                                                                                                         >  1:1375105/1‑225 (MQ=32)
                                     gggCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGAACCCCGAACAGGCCCCGCCGGGGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAAAAAtcccccccccat                                                                                                                                                        >  2:1573724/1‑227 (MQ=11)
                                      ggCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATa                                                                                                                                                       >  1:1614439/1‑238 (MQ=32)
                                       gCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCAt                                                                                                                                                        >  1:1343200/1‑236 (MQ=32)
                                               cccccACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCCCCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGAGGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATg                                                                                                                                             <  1:1600947/238‑1 (MQ=17)
                                                 cccACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGt                                                                                                                                            >  1:90548/1‑238 (MQ=32)
                                                 cccACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTc                                                                                                                                           >  1:262236/1‑239 (MQ=32)
                                                    aCCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGc                                                                                                                                         >  1:622363/1‑238 (MQ=33)
                                                        gcgcGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGACCAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATCCCACCGCCATACTGCCGCATGCCGCCCTg                                                                                                                                     >  2:1442136/1‑238 (MQ=255)
                                                             ggCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCAGACCGAACAGGCCACGCCGGTAGGGGCACCAATCCGCACCCCGGCTGAACCCGGAACAATCCAACCGCAATACTGCGGCATGTCGCCCTGAcccg                                                                                                                                >  2:179309/1‑238 (MQ=11)
                                                             ggCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATa                                                                                                                                                       >  1:1076458/1‑215 (MQ=32)
                                                              gCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATa                                                                                                                                                       >  1:782386/1‑214 (MQ=32)
                                                                            tAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCt                                                                                                                                                                                 <  1:66422/174‑1 (MQ=33)
                                                                            tAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCt                                                                                                                                                                                 >  2:66422/1‑174 (MQ=33)
                                                                            tAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGc                                                                                                                 >  1:1029594/1‑238 (MQ=35)
                                                                             aaCTGCTGTGGGTGGGGGAGAGGATGAACGGCGGCGGGCGGGGCAGGGCGATGTCGGTTTCACCGGGAAGGCCCTTT‑GCCCC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCGGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCt                                                                                                                <  1:59013/238‑1 (MQ=255)
                                                                             aaCTGCTGGGGGGGGGGGAGAGGATGAACGGGGGCGGGCCGGGCGGGGCGGTGTGGGTTTCGCGGGGAGGGCCCTTT‑GCGCC‑GGGGGGGGGTTTGCCG‑GGATGGCGGGACGCCCCCCCGCCGGGGCACCCCAAACGGCCCACGCCGGTGGGGGCACCATTGCGCCCCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCt                                                                                                                <  2:1629339/238‑1 (MQ=2)
                                                                                 gctgTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATa                                                                                                                                                       >  1:1490661/1‑195 (MQ=33)
                                                                                     tAGGTGCGGGAGAGGATGAACGGCAGCGGGACGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACATTT‑GCACC‑GCGCAGGGGGTTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGTACCGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGCACCCTGGACAATGCTACCGCCATACTACGTCATGTCGACCTGCCGCGCTGCCGTGTTTCCCCTCATTtccc                                                                                                        >  2:1098450/1‑236 (MQ=11)
                                                                                       ggTGCGGGGGAGGAGGACGGGCGGCGGGCCGGGCGGGGGGTTGTCGTTTTCCCCGGGAGGGCCCTTC‑GCCCC‑GGGCGGGGGATTGCCC‑GGATGGCGGGACGTCCCCCCGCGGGGGCCCCCCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACATTGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCAt                                                                                                      <  2:338192/238‑1 (MQ=14)
                                                                                          gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGaa                                                                                                                                                                              >  1:1505425/1‑163 (MQ=33)
                                                                                          gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCtt                                                                                                  >  1:1406818/1‑239 (MQ=40)
                                                                                          gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCtt                                                                                                  >  1:1222404/1‑239 (MQ=40)
                                                                                          gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCtt                                                                                                  >  1:1363279/1‑239 (MQ=40)
                                                                                          gCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCtt                                                                                                  >  2:1745720/1‑239 (MQ=40)
                                                                                                  ggATGAAGGGCGGGGGGCCGGGCAGGGGGAGGTCGTTTCCCCCGGGAGGGCCCTTC‑GCCCC‑GGGGAGGGGATTGCCC‑GGAGGCCGGGACGTCCCCCCGCCGGGGCACACCGAACGGGCCACGCCGGTGGGGGCACAAAGGCGCACCCCGGCTGACCCCGGAACATTGTTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTAt                                                                                           <  2:846807/238‑1 (MQ=14)
                                                                                                  ggATGAACGGCAGCGGGCCGGGCGGGGCGATGTCGGTTTCACCGGGAAGGCCCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCGGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTAt                                                                                           <  1:354944/238‑1 (MQ=255)
                                                                                                   gATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGCTGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGG‑TTTCCGCTCATTTCTCATCCTTGTTTTATCgg                                                                                        >  2:361523/1‑239 (MQ=255)
                                                                                                        aCGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTa                                                                                    >  1:870961/1‑239 (MQ=255)
                                                                                                        aCGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTa                                                                                    >  1:303119/1‑239 (MQ=255)
                                                                                                          ggcgggggcgcggggAGGGGGGTGTGGGTTCCCCCGGGGAGCCCCTTG‑GCCCG‑GAGGAGGGGATGGCCG‑GGAGGGGCGGGCGCCCCCCCCCCGGGGACCCCCGAACAGGCCCCGCGGGGGGGGCCACCAATGCGCCCCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTAtt                                                                                  <  2:34535/229‑1 (MQ=9)
                                                                                                               cgggccgggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATa                                                                                                                                                       >  1:844187/1‑165 (MQ=33)
                                                                                                                 ggccgggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTct                                                                             >  1:1656522/1‑237 (MQ=255)
                                                                                                                 ggccgggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTct                                                                             >  1:940081/1‑237 (MQ=255)
                                                                                                                      ggcAGGGGGATGTGGGTTTCCCCGGGAAGGCCCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGAGGGCCGGACGTCCCTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGAc                                                                       <  1:1315722/238‑1 (MQ=255)
                                                                                                                      ggcAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACa                                                                      >  1:1266732/1‑239 (MQ=255)
                                                                                                                          gggggggtggggTTTCCCCCGGGAGGACCTTT‑GCCCCGGACGGGGGGATGGACG‑GGTGGGCGGGCCTGCCCCCCGCCGGGGCCCCCCAACCGGGCACCGCGGGGGGGGGCCCAAAGGGGCCCCCGGGCTGACCCCTGAACAAGTCTCCCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAgg                                                                    <  2:920572/228‑1 (MQ=9)
                                                                                                                           ggCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGa                                                                 >  1:1643791/1‑239 (MQ=255)
                                                                                                                            gggATGTCGGTTCCCCCGGGAAGGCCTTTG‑GCACC‑GAGCAGGGGATTGCCG‑GGAGGGCCGGACGTCCCCCCGCGGGGGCACCCCGAACAGGCCCCGCCGGTGGGGGCACCAATGCGCCCCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGa                                                                 <  1:223536/236‑1 (MQ=14)
                                                                                                                            gCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCCCC‑GAGCAGGGGATTGCCG‑GGAGGGCCGGACGTCACCCCGCCGGGGCACCCCGAACAGGCCACGCCGGTGGGGGCACCAAGGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGa                                                                 <  1:1151404/238‑1 (MQ=255)
                                                                                                                            gCGATGGCGGTTTCCCCGGGAAGGCCTTTT‑GCCCC‑GAGCGGGGGATTGCCC‑GGAGGGCCGGACGTCACTCCGCGGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCCCCCCGGCTGAACCCGGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGa                                                                 <  1:853709/238‑1 (MQ=14)
                                                                                                                              gtgggcGGTTACA‑CGGGAGGGCGCTTT‑GCCCCGGGGCGGGGGAGTGCCG‑GGGGGGCCGGCCGCCCCCCCGCGGGGCCCCCCCGACCGGGCCACGCGGGTGGGGGCACCAATCCGCCCCCCGGCTGACCCCTGAACATTGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGaaa                                                               <  2:1330375/233‑1 (MQ=2)
                                                                                                                                   cgtttcccccGGGAAGCCCCTTG‑CCCCG‑GGGCGGGGGATGGCCC‑GGGGGGCGGGACGCCCCCCCGCGGGGCCACACCGAACGGGCCCCGCGGGGGGGGGCACCAAGGCGCCCCCCGGCTGACCCCTGAACAATGCTTCCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTAc                                                         <  2:126235/231‑1 (MQ=14)
                                                                                                                                   cGGTTTCACGGGGAAGGACCTTT‑GCCCC‑GAGCAGGGGTTTGACC‑GGATGGCCGGACGCCACCCCGCGGGGGCACCCCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCCCCCCGGCTGACCCCTGAACAATGCTCCCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTAc                                                         <  1:603216/239‑1 (MQ=14)
                                                                                                                                   cGGTTTCACCGGGAAGGACCTTT‑GCGCC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACCCCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTa                                                          <  1:179217/238‑1 (MQ=255)
                                                                                                                                      tttCACCGGGAGGGCCCTTT‑GCACC‑GAGCGGGGGATTGACC‑GGAGGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAg                                                       <  1:413019/238‑1 (MQ=255)
                                                                                                                                       ttccacGGGGAGGACCTTTGGCACG‑GGGGAGGGGATTGACA‑GGAGGGGCGGCGCTCCCCCCCCGGGGGCCCCCCGAACAGGCCCCGCCGGTGGGGCCACCAATGCGCCCCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAg                                                       <  2:1407289/233‑1 (MQ=14)
                                                                                                                                          aCCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATa                                                                                                                                                       <  1:1329910/138‑1 (MQ=34)
                                                                                                                                              ggAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAg                                              >  1:731097/1‑239 (MQ=255)
                                                                                                                                               gAAGGCCCTTG‑GCCCC‑GAGCAGGGAATTGACC‑GGAGGGCGGGGCGTCCCCCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAg                                              <  1:1447426/238‑1 (MQ=14)
                                                                                                                                                         tggcgcc‑GAGCGGGGGTTTACCG‑GGATGGCCGGGGGTCACCCCGCGGGGGAACACCGAACGGGCCACGCCGGTGGGGGCCCCAATGCGCCCCCGGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGaa                                    <  1:216420/234‑1 (MQ=14)
                                                                                                                                                         t‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCAccc                                                                                                                                                                                      >  1:1556524/1‑92 (MQ=37)
                                                                                                                                                            cccc‑GAGCGGGGGATTGCCC‑GGATGGCGGGACGTCCCTCCGCCGGGGCCCACCGAACGGGCCACGCCGGTGGGGGCCCCAATGCGCCCCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTg                                 <  1:126657/237‑1 (MQ=14)
                                                                                                                                                             ccc‑GAGCAGGGGATTGCCC‑GGAGGGCCGGACGCCACTCCGCGGGGGCCCACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTg                                 <  1:353105/237‑1 (MQ=255)
                                                                                                                                                             acccgggagggggTTTGCCC‑GGCGGGCCGGCCTGCCCCCCCCCGGGCCCCACCGAACAGGCCCCGCCGTGGGGGCCACCATGCCGCCCCCCGGCTGAACCCTGAACAAGGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGACGAAACCCTACAGATTTATAAGAAAGCCTGAAGTg                                 <  2:303336/230‑1 (MQ=14)
                                                                                                                                                             aCC‑GGGCGGGGGATGGCCG‑GGATGGCGGGACGTCGCTCCGGCGGGGCCCACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCCCCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCCGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTg                                 <  1:494825/238‑1 (MQ=14)
                                                                                                                                                             aCC‑GAGCGGGGGATTGACC‑GGATGGCCGGCCGTCCCTCCGCCGGGGCACACCAAACAGGCCACGCCGGTGGGGGCACCAATGCGCCCCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTg                                 <  1:102113/238‑1 (MQ=255)
                                                                                                                                                              cc‑GGGCAGGGGGTTAGCA‑GGATGGCCGAATGCCCCACCGGCGAGAAACACCGAACAGGCCACGCCGGTGGGGGCCCAAATGCGCCCCCCGGATTAACCCTGAACATTGTTTCCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTgag                               <  2:1363279/239‑1 (MQ=14)
                                                                                                                                                              cc‑GGGCAGGGGATGGCCC‑GGAGGGCGGGAGCTCCCCCCGCCGGGGCCCACCGAACGGGCCAGGCGGGTGGGGGCACCAAGGCGCCCCCCGGCTGACCCCTGAACAAGGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTgag                               <  2:579777/239‑1 (MQ=14)
                                                                                                                                                               cggAGCGGGGGATTGACC‑GGGGGGCGGGGCGCCCCCCCCCCGGGGCCCCCCGAACAGCCCACGCCGGGGGGGCCACCAATGCCCCCCCCGGCTGACCCCGGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTgag                               <  2:940081/237‑1 (MQ=14)
                                                                                                                                                                 gAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCa                                                                                                              >  2:891389/1‑158 (MQ=38)
                                                                                                                                                                  agggggggATTGACC‑GGATGGCGGGAGGTCACCCCGCGGGGGCACACCGAACGGGCCACGCCGGTGGGGGCACAAATGCGCCCCCCGGCTGAACCCTGAAAAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGat                             <  1:47472/234‑1 (MQ=14)
                                                                                                                                                                  aGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTc                                                                              >  1:597673/1‑189 (MQ=255)
                                                                                                                                                                  aGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTgag                               >  1:878206/1‑236 (MQ=255)
                                                                                                                                                                  aGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGata                            >  1:522535/1‑239 (MQ=255)
                                                                                                                                                                  aGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGata                            >  2:267928/1‑239 (MQ=255)
                                                                                                                                                                  aGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGata                            >  1:562990/1‑239 (MQ=255)
                                                                                                                                                                  aGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGata                            >  1:790169/1‑239 (MQ=255)
                                                                                                                                                                   ggaggtgcatGACCGGGATGGCCGGGCGCCCCTCCGCCGGGCCACCCCGAACCGGCCCCGCCGGTGGGGGCACCAAGGCGCACCCCGGCTGAACCCTGAACAATGATACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGat                             <  2:1367631/229‑1 (MQ=14)
                                                                                                                                                                        ggATGGACG‑GGAGGGCCGGACGCCCCCCCGCGGGGGCACCCCAAACAGGCCACGCGGGTGGGGCCCCCAATGCGCCCCCCGGCTGACCCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCAt                       <  1:1442136/238‑1 (MQ=14)
                                                                                                                                                                          ttggccc‑GGAGGGCCGGCCCCCCCCCGGCCGGGGCACACCGACCGGCCCACGCCGGGGGGGGCACCAATGCGCACCCCGGCTGAACCCGGAACAAGGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCa                     <  2:262236/233‑1 (MQ=14)
                                                                                                                                                                            tGACC‑GGATGGCCGGACGTCATTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCagag                  >  2:1269201/1‑239 (MQ=255)
                                                                                                                                                                            tGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGtttt                                                                                                                     >  1:1041863/1‑140 (MQ=40)
                                                                                                                                                                            tGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCagag                  >  1:128040/1‑239 (MQ=255)
                                                                                                                                                                            tGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCagag                  >  1:314693/1‑239 (MQ=255)
                                                                                                                                                                                   gATGACCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGCTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCCCATTTCTCCCCCTTGTTTTATCgg                                                                                        >  2:385182/1‑163 (MQ=255)
                                                                                                                                                                                    aTGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGaa                                                                <  1:1104490/186‑1 (MQ=255)
                                                                                                                                                                                       gCCGGACCTCACCCCGCGGGGGCACACCGAACGGGCCACGCGGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCAGAGTGAaataat         <  2:1656522/238‑1 (MQ=255)
                                                                                                                                                                                         cGGACGTCACTCCGCCGGGGCACACCGAACAGGCCTCGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGtt                                                                                               >  2:327777/1‑150 (MQ=255)
                                                                                                                                                                                         cGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGtt                                                                                               <  1:327777/150‑1 (MQ=255)
                                                                                                                                                                                           gccgccACTCCGCCGGGCCACCCCGACTGCACCCCGCTGGTGGAGGCCCAAAGCCGCCCCCCGCCTGAACCCTGAGCAATTTTTCCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCAGAGTGAAATAATTAgg     <  2:1029594/233‑1 (MQ=17)
                                                                                                                                                                                             cGTCCCTCCGCCGGGGCCCCCCAAACAGGCCACGCCGGTGGGGGCACCAATGCGCCCCCCGGCTAAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCAGAGTGAAATAATTAGGaa   <  1:179309/238‑1 (MQ=255)
                                                                                                                                                                                                   |                                                                                                                                                                                                                                        
CAGCTGTAAGCGGATATCCGCAGGCATTTTCCAGCCGGGGCCGAGGCTCCCCACCGGCGCGGGCGTTTTTGTCCGGTAACTGCTGTAGGTGCGGGAGAGGATGAACGGCAGCGGGCCGGGCAGGGCGATGTCGGTTTCACCGGGAAGGACCTTT‑GCACC‑GAGCAGGGGATTGACC‑GGATGGCCGGACGTCACTCCGCCGGGGCACACCGAACAGGCCACGCCGGTGGGGGCACCAATGCGCACCCCGGCTGAACCCTGAACAATGCTACCGCCATACTGCGTCATGTCGCCCTGACGCGCTGCCGGTTTTCCGCTCATTTCTCATCCTTGTTTTATCGGGTTATTTGTCTGCTGACAGGCGAAACCCTACAGATTTATAAGAAAGCCTGAAGTGAGATATCCATCAGAGTGAAATAATTAGGAAA  >  NC_007779/3877915‑3878339

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