Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,346,6450GA27.3% 79.4 / 11.6 44S33F (TCT→TTT) pdeRcyclic di‑GMP phosphodiesterase PdeR
Reads supporting (aligned to +/- strand):  ref base G (23/9);  new base A (1/11);  total (24/20)
Fisher's exact test for biased strand distribution p-value = 3.25e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAA  >  NC_000913/1346400‑1346885
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
gCTCAACAAAGAAAACCCGGTGTGATAAATACCGTGCAACCGCCGGTGTGTCAGGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAATATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTTGTTTCGGTGG‑TAGaaa                                                                                                                                                                                                                                                  <  2:432827/249‑1 (MQ=255)
     gCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCa                                                                                                                                                                                                                                           >  1:257312/1‑250 (MQ=255)
         aGGAAAGCCCGGGTGATAACTCACGAGAAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACa                                                                                                                                                                                                                                      <  1:192618/251‑1 (MQ=255)
                   gtgtGATAAATCACGAGCAACCGCCGGTGTGTCATGACATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCATACCAGA‑TGCACAGAAACATC‑GCTTTCCTCGACGGTCC‑GACTCAT‑CATCAGGCTGGAGGTGATGGAC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGAACACTGCAATGA‑ACGATCAGGGTCTGTGGTTTCGGTGG‑TAGAAAA‑GCAAAAAACATCCCTGCttt                                                                                                                                                                                                                            >  2:694046/1‑251 (MQ=255)
                    tgtgATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTc                                                                                                                                                                                                                           >  1:499177/1‑251 (MQ=255)
                     gtgAAAAACCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTc                                                                                                                                                                                                                           <  1:545813/250‑1 (MQ=255)
                        aTAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTc                                                                                                                                                                                                                       <  1:83887/251‑1 (MQ=255)
                                gaaGCACCCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCa                                                                                                                                                                                                             <  1:81526/249‑1 (MQ=255)
                                 ggcccaCCGCCGGTGTGTCATGCCATTCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGTTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAg                                                                                                                                                                                                            <  2:493242/246‑1 (MQ=255)
                                 gAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCAATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑Gcc                                                                                                                                                                                                              >  2:264331/1‑249 (MQ=255)
                                            ggTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGAAATAT‑AAGATCAGGTTCTGTGGTTTCGGTGG‑TAACCAA‑ACGAAAAACATCGATTCTTTCAGTC‑AATC‑GCCCGTAGGGAGtgtg                                                                                                                                                                                                 >  2:112197/1‑251 (MQ=255)
                                                 gtCAGGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGCCCTACCAGA‑TGCAAAAAAAGAGC‑GCTTTCTCCGAGGGTCA‑AACT‑ATACCTAGGGCGGGAGGTGATGACC‑GTCTTTCCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACGGAAATTA‑ACGATAAGGTTTTGGGTTTTGGGGGG‑TAGAAAA‑GCAAAAAACATCGCTTTTTCCCGCC‑AACC‑GCCAGTATGGATTGTGCGAt                                                                                                                                                                                             <  1:417014/250‑1 (MQ=255)
                                                                 gcgTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACa                                                                                                                                                                             >  2:751109/1‑250 (MQ=255)
                                                                      cccacccaaTTCCGGTTTTTTGTTTTTTC‑TCCTACCAAAAGGCACAGAAA‑ATCCGCTTTCATCGACG‑TCAAGACTCAT‑CATCGGGTTGGGGGTTA‑GACCCTTCATTTCCCTG‑AGGCGAGCAGCTTTTTCGG‑CAGAAAACTGAAATTG‑ACGTTCTGTTTTTGGGTTTTGGGTGG‑TAAAAAA‑CCGAAAAACACCGCTGTTTCCCATC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGtt                                                                                                                                                                        <  1:722993/242‑1 (MQ=255)
                                                                                 cccGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCt                                                                                                                                                             >  1:575805/1‑250 (MQ=255)
                                                                                       tATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGtt                                                                                                                                                      >  2:182708/1‑251 (MQ=255)
                                                                                        aTTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTACTCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACA‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGAATGGATTGTGCGATCCGAGCACGTTGTACAACGTTGTGGACTACCTAACGGttt                                                                                                                                                     >  2:326631/1‑251 (MQ=255)
                                                                                                   gTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGAACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCat                                                                                                                                          >  1:306902/1‑251 (MQ=255)
                                                                                                           aacatgaaaaaaaaGGTC‑GCTTTCGTCGACGTGCA‑GACTCAT‑CACCGGGCGTGAGTGGATGCCC‑TTC‑TTTCCCTGTATGGGAGCGGCCT‑TTCGGGCAGATAAC‑GCAATTTTAGGATGAGGTTGTGTGGTTGGGGGGG‑TAGAAAG‑GCGAAAAACCTCGCTGCTTTCCTTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATAtcctcc                                                                                                                                   <  2:269068/240‑1 (MQ=255)
                                                                                                             a‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGCGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGt                                                                                                                                >  1:140756/1‑251 (MQ=255)
                                                                                                                               ttttCATCGACGTTCA‑GCCTAAC‑CATCGGGGGGGGGGTGAGGCCC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTGCGG‑AAGATAACGGAAATGT‑CGGAACGGGTTTTGGGGTTT‑GGTGGGAAAAAAA‑GCGAAAAACACCGCTGCTTCCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTcc                                                                                                                  <  1:590688/248‑1 (MQ=255)
                                                                                                                                 ttCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCgg                                                                                                              >  1:660372/1‑251 (MQ=255)
                                                                                                                                       gACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaact                                                                                                        >  2:167234/1‑249 (MQ=255)
                                                                                                                                       gACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGAGGTtaaat                                                                                                        >  1:509205/1‑251 (MQ=255)
                                                                                                                                          ccggccagcctcatcagcgggggggAGGTGATGACC‑GTCATTCCCGTG‑ATGCGAGCAGCCTTTTCGG‑CAGATACCTGAAATGT‑ACGTTCAGGTTCTGTGTTTTGGGTGG‑AAAAAAA‑GCGAAAACCACCCCTGTTTCCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaata                                                                                                       <  1:339316/228‑1 (MQ=255)
                                                                                                                                           gTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTtaaataaa                                                                                                     >  1:465890/1‑250 (MQ=255)
                                                                                                                                                     at‑catcaGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGaacaac                                                                                            >  2:615991/1‑250 (MQ=255)
                                                                                                                                                            ggggggggggTGATGACC‑GTCTTTTCCCGG‑AGGCGAGCTGCCTTTTCGG‑AGGAAAACTGCAATGT‑ACGCTCGGGTTCTGTGTTTTGGGTGG‑TAAAAAA‑GCGCAAAACTCCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACtt                                                                                     <  1:615991/243‑1 (MQ=255)
                                                                                                                                                                    ggTGATGACC‑GTCTTTTCCTTG‑AGGCGACCAGCCTGTTCGG‑AAGATAACTGAATTGT‑CGGATCGGGTTTTGGGGTTTCGGTGG‑TAAAAAA‑GCGCAAAACACCGCTGCTTCCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGaaa                                                                             <  1:750916/251‑1 (MQ=255)
                                                                                                                                                                                   cccgatcggaggcgggcagcCTGTTCGG‑CAGATAACTGAAATGT‑ACGATCAGGTTCTGTGGTTTGGGTGG‑TAAAAAA‑GCGAAAAACATCGCTGCTTCCCGTA‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTa                                                                <  1:748340/236‑1 (MQ=255)
                                                                                                                                                                                      cccTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGCAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCtcaa                                                              >  2:272560/1‑246 (MQ=255)
                                                                                                                                                                                          g‑ATGCGAGCAGCTTTTCCGG‑CAAATAACGGCAATGT‑CCGATCAGGTTCTGTGTTTTCGGTGG‑TAAAAAA‑GCGCAAACCATCGTTTTTTCCCGTC‑ATCC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCa                                                          <  1:514655/249‑1 (MQ=255)
                                                                                                                                                                                          g‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTa                                                                >  2:748340/1‑243 (MQ=255)
                                                                                                                                                                                            aTGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTg                                                       >  1:574844/1‑251 (MQ=255)
                                                                                                                                                                                                    agcCTGTTGGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGGGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCACTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAATCATAGtttttttt                                               >  2:513196/1‑251 (MQ=255)
                                                                                                                                                                                                     gcCTGTTGGG‑AAGAAACCTGAAATGT‑ACGATCGGGTTCTGTGGTTTGGGGGG‑AAAAAAA‑GCGAAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTta                                              <  1:182708/251‑1 (MQ=255)
                                                                                                                                                                                                         tgtggggcagaaacgttcAATGT‑ACGATCGGTTTTTGTGGTGCGGGTGG‑TAAAAAA‑GCGAAAAACATCGCTTCTTTCCG‑C‑AACCCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTtata                                            <  1:751109/233‑1 (MQ=255)
                                                                                                                                                                                                                    aaacctgaattatccgttaaGGTTCTGTGGTTTGGGTGG‑AAAAAAA‑GCGAAAAACAGCGCTGCTTCCCGAC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGATTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGt                                 <  1:264331/232‑1 (MQ=255)
                                                                                                                                                                                                                     gaacatcaatttccaaTCAGTTTTTGTGGTTTCGGTGT‑AAGAAAA‑GAGAAAAACTTCGCT‑CCTGCCCTCTACCC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGt                                 <  1:530907/235‑1 (MQ=255)
                                                                                                                                                                                                                          ccaatttacgctcggtttCTGTGTTTTCGG‑GG‑TAAAAAAAGCAAAAAACATCGCTG‑TTTCCCTT‑AACCCGCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCa                           <  1:272560/235‑1 (MQ=255)
                                                                                                                                                                                                                                       aGGTTCTGTGGTTTCGGTGG‑TAGAAAC‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCCGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGGTCAACGACTTGTTTTTTTTATTTATTTTTTTGTAATCAATTTTTGGtaaaa               >  2:483421/1‑250 (MQ=255)
                                                                                                                                                                                                                                              gtgGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATAGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAAAGCATTGTTTTTTTTATATATTTTTTTGTAATCCAGTTTTGGAAAAACGCCAGtt       >  1:598535/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  tttCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTTTATTTTTATTTTTTATCCAGTTTTGGAAAAATGCCCGTTTTc    >  2:463732/1‑250 (MQ=255)
                                                                                                                                                                                                                                                    tCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTTTATATTTATTTGTAATCCAGTTTTTGaaaaaagcacgtttttta  >  2:247201/1‑237 (MQ=255)
                                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
GCTCAGCAAAGGAAAGCCCGTGTGATAAATCACGAGCAACCGCCGGTGTGTCATGCCATGCAGACGCGTTGCCAGCCCATTCCCGTTTATTGATTTTTCGTCCTACCAGA‑TGCACAGAAAGATC‑GCTTTCATCGACGGTCA‑GACTCAT‑CATCAGGCTGGAGGTGATGACC‑GTCATTTCCCTG‑ATGCGAGCAGCCTGTTCGG‑CAGATAACTGCAATGT‑ACGATCAGGTTCTGTGGTTTCGGTGG‑TAGAAAA‑GCGCAAAACATCGCTGCTTTCCGTC‑AACC‑GCCAGTATGGATTGTGCGATCCGAGAAAGTTGTACAACGTTGTGGACTCCCTAACGGTTTTCATGATGCATATCCTCCGGTTAACAGCGGAGTCCGCGGTTAAATAAAAGGAACAACATGACTTACATGAAATTAACGGCGGCTAAACGCATTGTTTTTTTTATATATTTATTTGTAATCCAGTTTTGGAAAAACGCCAGTTTTCAA  >  NC_000913/1346400‑1346885

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