Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,040,6230GA64.9% 48.8 / 92.0 94noncoding (107/120 nt)rrfA5S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (15/18);  new base A (37/24);  total (52/42)
Fisher's exact test for biased strand distribution p-value = 1.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-02
Rejected as consensus: Frequency below/above cutoff threshold.

CGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATG  >  NC_000913/4040377‑4040864
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
cGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACt                                                                                                                                                                                                                                               >  2:648992/1‑251 (MQ=14)
    gCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGcc                                                                                                                                                                                                                                            <  1:368863/250‑1 (MQ=14)
      gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAgg                                                                                                                                                                                                                                         >  1:198900/1‑251 (MQ=14)
      gTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAgg                                                                                                                                                                                                                                         >  1:320585/1‑251 (MQ=14)
        tGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGACGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGc                                                                                                                                                                                                                                        >  2:327982/1‑250 (MQ=14)
              aaCCGGTACTAATGAACCGTGCGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATc                                                                                                                                                                                                                                     >  2:513968/1‑247 (MQ=14)
              aaCCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATc                                                                                                                                                                                                                                     <  1:513968/247‑1 (MQ=14)
                   gTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAAc                                                                                                                                                                                                                            <  2:109877/251‑1 (MQ=18)
                        aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGaa                                                                                                                                                                                                                         >  2:496929/1‑249 (MQ=25)
                        aaTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGAGTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACGCAGAAGTGAAACGCCGGAGCGCCGATGGTAGTGTGGGG‑GCTCCGCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                                                                                                                       <  2:81705/251‑1 (MQ=255)
                         aTGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAgg                                                                                                                                                                                                                      <  2:579500/251‑1 (MQ=25)
                           gAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTGCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCt                                                                                                                                                                                                                    <  2:60195/251‑1 (MQ=25)
                            aaCCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAAGCATCAAATAAAACGAAAGGCTc                                                                                                                                                                                                                   >  1:711515/1‑251 (MQ=25)
                               cGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGt                                                                                                                                                                                                                <  2:425072/251‑1 (MQ=25)
                                gTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTc                                                                                                                                                                                                               >  1:41706/1‑251 (MQ=25)
                                gTGAGGATTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTc                                                                                                                                                                                                               >  2:218541/1‑251 (MQ=25)
                                     gCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAg                                                                                                                                                                                                          <  2:530613/251‑1 (MQ=38)
                                          aCCTTACAACGCCGACGGTGTTTTGGCGGCTTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGCAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGACCTCAGATGTGAAACGCCGTAACGCCGCTGGCAGTGTGGGG‑TCTCCTCATGCGGGAGTAGGAAACTGCCACGCATCACATACAACGAAAGGCTCAGTCGGAAGACTgg                                                                                                                                                                                                     >  1:263958/1‑251 (MQ=11)
                                               aCAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCtt                                                                                                                                                                                                >  2:415158/1‑251 (MQ=38)
                                                cAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCttt                                                                                                                                                                                               >  1:57038/1‑251 (MQ=38)
                                                 aaCGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTc                                                                                                                                                                                              <  2:143473/251‑1 (MQ=38)
                                                             gtTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgtt                                                                                                                                                                                  <  2:198900/251‑1 (MQ=25)
                                                             gtTTTGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCtgatt                                                                                                                                                                                 >  2:276493/1‑248 (MQ=14)
                                                                 tGGCGGAT‑GAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgttt                                                                                                                                                                              <  2:135206/250‑1 (MQ=14)
                                                                   gCGGCTTGAGAGAAGATTTTCCGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttgt                                                                                                                                                                            >  2:47134/1‑251 (MQ=18)
                                                                    cGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCACATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTc                                                                                                                                                                           >  2:540321/1‑251 (MQ=25)
                                                                     cggatGAGAGAAGATTTACAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCTCCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGTTAGTGTGGGG‑TCGCCTCATGCGAGAGTAGGGAACTGCCAGTCATCAAATAATACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAACTGGTGTTTGTCg                                                                                                                                                                          <  2:419173/247‑1 (MQ=11)
                                                                      gattgcgagtaaactctCAGCCTAATACAGATTAAATAACAACGCAGAAGCGGCCTGATAAAAACGACTTTGCCTAGAGGCAGTAGCGCGGTGGTCCAAACCGACCCCATGCCGAACTCAGAAGTGAAACGCCGTACAGCCGACGATAGAGTGGGG‑TCTCCTCATGCGAGAGAAGGAAACTGCCAGGAATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCgg                                                                                                                                                                         >  2:367625/17‑251 (MQ=11)
                                                                       aTTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGt                                                                                                                                                                        >  1:297648/1‑251 (MQ=18)
                                                                        atGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGGCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTg                                                                                                                                                                       <  2:259896/250‑1 (MQ=17)
                                                                          gagagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGa                                                                                                                                                                      >  2:601616/1‑250 (MQ=33)
                                                                          gagagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGaa                                                                                                                                                                     <  1:500286/251‑1 (MQ=18)
                                                                            gagaAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                                                                                   <  2:206849/251‑1 (MQ=18)
                                                                            gagaAGAGTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTATCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGGAGGAAACTGCCAGGCAGCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACg                                                                                                                                                                   <  1:737059/251‑1 (MQ=255)
                                                                             agaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGc                                                                                                                                                                  >  2:41703/1‑251 (MQ=18)
                                                                              gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGct                                                                                                                                                                 >  1:154083/1‑251 (MQ=33)
                                                                              gaagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGct                                                                                                                                                                 <  1:130483/251‑1 (MQ=18)
                                                                              gaagaTTTTCAGACTGATACAGATTAAATCAGAATGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGAGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGATGTTTGTCGGTGAACGct                                                                                                                                                                 >  2:456913/1‑251 (MQ=255)
                                                                               aagaTTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAAATTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAACACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGctc                                                                                                                                                                >  2:604942/1‑251 (MQ=18)
                                                                                       cAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTCGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGt                                                                                                                                                        >  2:630736/1‑251 (MQ=18)
                                                                                       cAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAg                                                                                                                                                         <  2:209263/250‑1 (MQ=33)
                                                                                       cAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGt                                                                                                                                                        <  1:496929/251‑1 (MQ=18)
                                                                                         gCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAg                                                                                                                                                      <  1:81496/251‑1 (MQ=18)
                                                                                                 cAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATc                                                                                                                                              >  1:523096/1‑251 (MQ=33)
                                                                                                       aaaTCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATccgccg                                                                                                                                         <  1:244762/250‑1 (MQ=18)
                                                                                                          tCAGAACGCAGAAGAGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCggg                                                                                                                                       >  1:412827/1‑249 (MQ=18)
                                                                                                           cAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGc                                                                                                                                    <  2:309434/251‑1 (MQ=18)
                                                                                                              aaCGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCgg                                                                                                                                  >  2:39541/1‑250 (MQ=18)
                                                                                                                 gCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAttt                                                                                                                              <  1:41703/251‑1 (MQ=18)
                                                                                                                        cGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGtt                                                                                                                       <  1:215122/251‑1 (MQ=33)
                                                                                                                               aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGc                                                                                                                >  2:567802/1‑251 (MQ=33)
                                                                                                                               aTAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGc                                                                                                                <  1:415393/251‑1 (MQ=33)
                                                                                                                                 aaaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                                                                                                                       >  1:440382/1‑146 (MQ=32)
                                                                                                                                 aaaaCAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAg                                                                                                                                                                                                                       <  2:440382/146‑1 (MQ=32)
                                                                                                                                 aaaaCAGAAGTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGGCCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCGCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCaa                                                                                                              <  2:183626/251‑1 (MQ=14)
                                                                                                                                  aaaCAGAATTTGCCGGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTTAACGTTGCGAAGCAAc                                                                                                             <  1:631941/251‑1 (MQ=18)
                                                                                                                                     cAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAATTGAAACGCCGTAGCGCCGATGGTAGGGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGTCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGGTGTGTGTCGGTGAACGCGCTCTTGAGTAGGACTAATCCGCCGGGATCGGATTTGAACGTTGCGAAGCAACGGc                                                                                                          <  2:720104/251‑1 (MQ=11)
                                                                                                                                       gAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGCTGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGCGTAGGACAAATCACCCGTGAGCGGATTTGAACTTTGCGAAGCAACGGccc                                                                                                        >  1:445316/1‑251 (MQ=11)
                                                                                                                                         aTTTGCCTGGCGGCAGTAGCGAGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCgg                                                                                                      >  2:418785/1‑251 (MQ=18)
                                                                                                                                              ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGaa                                                                                                                                                                                                                                                 <  1:78237/107‑1 (MQ=33)
                                                                                                                                              ccTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGaa                                                                                                                                                                                                                                                 >  2:78237/1‑107 (MQ=33)
                                                                                                                                               cgggcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                <  2:153933/249‑1 (MQ=18)
                                                                                                                                               cTGGCGGCAGTCACGCGGGGGTCCCACCTGACCCCATGCCGAACTCAGCAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGCGAGTAGGAAACTGCCAGGCATCAACTAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACCCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                >  2:621712/1‑251 (MQ=11)
                                                                                                                                               cTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAggggg                                                                                                >  2:276496/1‑249 (MQ=18)
                                                                                                                                                tGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTgg                                                                                               >  2:503033/1‑251 (MQ=33)
                                                                                                                                                 ggcggcAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTg                                                                                                <  2:381725/249‑1 (MQ=33)
                                                                                                                                                 ggcggcAGTAGCGCGGCGGTCCCACCTGACCCCATGCCGCACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGc                                                                                              >  2:438038/1‑251 (MQ=18)
                                                                                                                                                       aGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAg                                                                                        <  2:523096/251‑1 (MQ=33)
                                                                                                                                                            cgcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttg                                                                                                                                                                             <  1:421469/161‑1 (MQ=35)
                                                                                                                                                            cgcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGttgtttg                                                                                                                                                                             >  2:421469/1‑161 (MQ=35)
                                                                                                                                                            cgcgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCAGGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGTGCCTTTCGGTTGATCTGTTGTTTGTCGGTGATCGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGAGTTTAACGTTGCGAAGCAACGGCTCGGAGGGTGGCGGGCAGGAcgc                                                                                   <  2:320585/251‑1 (MQ=11)
                                                                                                                                                              cgGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgccc                                                                                 <  2:194620/251‑1 (MQ=33)
                                                                                                                                                               ggtggtCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATACGAGAGTAGGGCACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGAcgcccg                                                                                >  1:461188/1‑251 (MQ=18)
                                                                                                                                                                   gtccccccttccccagtGCCGAACTCAGAGGTGAAACTCCGTAGCGCCGAGGGTAGTGTGGGG‑TCTCCTCAGGCGAGAGTAGGGAACTGCAGGGCTTCAAAAAAAAGGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTAGCGGTTGTTTGTGGGTGAACGCTCCCCGGAGTAGGAAAAATCCGCCGGGAGGGGATTTGAACTTTGGGAAGCAACGGCCCGGAGGGTGGCGGTCAGGACGCCCGCCAt                                                                            <  2:428838/235‑1 (MQ=255)
                                                                                                                                                                      ccACCTGACCCCCTGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTCGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCCGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGACGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATaa                                                                          >  2:430058/1‑250 (MQ=255)
                                                                                                                                                                          cTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGc                                                                     >  1:510454/1‑251 (MQ=18)
                                                                                                                                                                               cccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc                                                                >  1:266442/1‑251 (MQ=18)
                                                                                                                                                                                ccATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGc                                                                >  1:38565/1‑250 (MQ=18)
                                                                                                                                                                                    gCCGAACCCAGAAGTAAACCGCCGTAGCGCGGATGGAAGTGTGGGG‑TCCCCTCAGGCGAGAGTAGGAAACTGCCAGGCATCAAAAAAAACAAAAGGCTCAGGCGGAAGACTGGCCCTTTCTTTTGATTTGTTGTTTGTGGGTGAACGCTCTCCTGAGAAGAAAAAATCCGCGGGGAGCTGATTTAAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCaa                                                           <  1:648753/251‑1 (MQ=11)
                                                                                                                                                                                     ccGAACTCAGAAATGAAACGCAGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAACTACAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGCCAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGAGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCAtaaa                                                           >  2:732342/1‑247 (MQ=11)
                                                                                                                                                                                       gAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGACGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAAtt                                                        >  1:678798/1‑251 (MQ=18)
                                                                                                                                                                                             aGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGTCGCCCGCCATAAACTGCCAGGCATCAAATTAAGc                                                    >  2:6200/1‑249 (MQ=18)
                                                                                                                                                                                               aaggaaaaCGCCGTAGCGCCGA‑GGTTGTGGGGGGGTCCCCTCAGGGGAGAGAAGGAACCCGCCAGACAAAAAAAAAAACGAGGGGCCCAGTGGGAAGCCTGCCCCTTGCTTTTTATCTGTTGTTTGTCGGAAAACGCCCCCCTGAGTAGGACAACGCCGCGGGGGGCGGATTTAAACGTGGAAAAGCACGGGCCCGGAGGGTGGCGGACAGTACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGaa                                                <  2:57038/246‑1 (MQ=2)
                                                                                                                                                                                                     aaCGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCAt                                          >  2:289523/1‑251 (MQ=33)
                                                                                                                                                                                                        gCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGCACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCt                                       >  1:301939/1‑251 (MQ=18)
                                                                                                                                                                                                           gTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGCGGGTGGCGGGCAGGACGCCCGACATAACCTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGAc                                    >  2:608517/1‑251 (MQ=18)
                                                                                                                                                                                                               cgcCGATGGTGGTGGGGGG‑TCTCCTCAGGGGAGAGTGGGGAACCGCCAGGCATCAAAAAAAACGAAGGGCTCAGTCGGAAGACGGGGCTTTTGGTTTTATCTGTTGTTTGTCGGAGAACGCTCTCCTGAGTAGAAAAAACCCGCCGGGAGCGTTTTTGAACGTTGCGAAGCAACGGCCGGGAGGGTGGGGGGCAGGCCGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGAt                                <  1:567802/251‑1 (MQ=11)
                                                                                                                                                                                                               cgcCGATGGTAGTGTGGGG‑TCTCCTCAGGCAAGAGTAGAAAACTGCCAGGCAACAAATAAAACGAAAGGCTTAGTCGGAAGACTGGGCTTTTCTTTTTATCTTTTGTTTGGGGGTGACCCCTCTCCGGAGTAGGACAAACCCGCCGGGGGCGGATTTGAACTTTGCGAAGCACCGCCCCGGAGGGTGGCGGGCGGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGAt                                <  1:604942/251‑1 (MQ=11)
                                                                                                                                                                                                                gcCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATg                               >  2:513191/1‑251 (MQ=18)
                                                                                                                                                                                                                 ccgctgggaGTGTGGGG‑TCTCCTCATGCGGGAGTAGGAAACTGCCAGGAATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGCGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGGGGCGGGCAGGCCGCCCGCAATAAACTGCCAGGCAGACAATTAAGCAGCAGGCAAGCCGGACGGATgg                              >  2:698054/9‑251 (MQ=11)
                                                                                                                                                                                                                    ggggtggtgggggg‑TCTCCTAATGCGAGAGTAGGAAACTGCCAGGCAAAAAAAAAAACGAAAGGCGCGGTGGGAAGACTGGTCTTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCCCATGAGTAGGAAAAACCCGCCGGGAGCGGATTTGAGGGTTGCGAAGCAACGGCCGGGAGGGGGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCt                           <  1:276496/241‑1 (MQ=11)
                                                                                                                                                                                                                    aTGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGACAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAAAGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGGGGCGGGCAGGACGCCCGCCAGAAAATGCCAGGCATCAAATTAAGCCGAAGGCCATCCCGCCGGATGGCCt                           >  2:441570/1‑251 (MQ=11)
                                                                                                                                                                                                                                             cGAGAGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCAGATTTGAACGTTGCGAAGCAACGGCACGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAAt   >  1:682042/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gagaGTAGGAAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATg  >  1:487345/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                 
CGATGCGTTGAGCTAACCGGTACTAATGAACCGTGAGGCTTAACCTTACAACGCCGAAGGTGTTTTGGCGGATTGAGAGAAGATTTTCAGCCTGATACAGATTAAATCAGAACGCAGAAGCGGTCTGATAAAACAGAATTTGCCTGGCGGCAGTAGCGCGGTGGTCCCACCTGACCCCATGCCGAACTCAGAAGTGAAACGCCGTAGCGCCGATGGTAGTGTGGGG‑TCTCCTCATGCGAGAGTAGGGAACTGCCAGGCATCAAATAAAACGAAAGGCTCAGTCGGAAGACTGGGCCTTTCGTTTTATCTGTTGTTTGTCGGTGAACGCTCTCCTGAGTAGGACAAATCCGCCGGGAGCGGATTTGAACGTTGCGAAGCAACGGCCCGGAGGGTGGCGGGCAGGACGCCCGCCATAAACTGCCAGGCATCAAATTAAGCAGAAGGCCATCCTGACGGATGGCCTTTTTGCATTGGCGCAGAAAAAAATG  >  NC_000913/4040377‑4040864

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