Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,130,3420AC38.2% 90.4 / 11.6 55G174G (GGT→GGG) flgAflagellar basal body P‑ring formation protein FlgA
Reads supporting (aligned to +/- strand):  ref base A (23/11);  new base C (1/20);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 3.21e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATCGCGGCAACCACATAATTTCCTGTGGCCTG  >  NC_000913/1130093‑1130590
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
tACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCa                                                                                                                                                                                                                                                          <  1:4849/250‑1 (MQ=255)
            agGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATcc                                                                                                                                                                                                                                             >  2:427917/1‑250 (MQ=255)
                gAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACAtt                                                                                                                                                                                                                                         <  2:215605/251‑1 (MQ=255)
                     gCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGAGCACATTGAcgc                                                                                                                                                                                                                                    <  1:206334/251‑1 (MQ=255)
                              aCTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGGCCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCg                                                                                                                                                                                                                           <  2:375746/251‑1 (MQ=255)
                                        tATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTAccc                                                                                                                                                                                                                  <  2:10806/250‑1 (MQ=255)
                                           tCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGcc                                                                                                                                                                                                               >  2:522047/1‑250 (MQ=255)
                                                                                aagcggaaaGGGAGGGGTAATCTTAACCAGCTAACAGTTTTAAAAAAATATCCCCTG‑CTGCATCA‑ACAACGCGGCTGACTACCTGCCCCGATACCATGCCCCCCCGCGCATTCTTTGCGAGGGCTGCATTTTTCGGCGCCTGCCCTTCTGCGTGGG‑CGCTAACCCCACCCCCGCTGGCGATCAAATGGACGTGTGGTCCCGTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATc                                                                                                                                                                          <  2:200360/243‑1 (MQ=255)
                                                                                     gaaaGGGACGGGTAATCTTTAACAGCTAACAGGTTTATAAAAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATCCCATGCGCCCCCGCGCATTCTGGTCGACGGCGTCATTGTTCAGCGCCTGCCCTTCGTCGTTGG‑CGCTAACCCCATCCCCGCTGGCGATCACTGTGCCGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGa                                                                                                                                                                    <  2:361401/249‑1 (MQ=255)
                                                                                      cAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTCCATGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGACACTGGCTTAACTGGATAGGTTGATCggg                                                                                                                                                                       >  2:463819/1‑247 (MQ=255)
                                                                                       aaGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGAt                                                                                                                                                                   >  1:170573/1‑250 (MQ=255)
                                                                                         gggACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATCCCCAT‑CTGCATCA‑ACACCGCCGCTGACTCCCTGCCCCGATACCATGCGCACCCGCGCTTTCTGGGCGACGGCTGCATTGTTCAGCGCCTGACCTTCGGCGTTGG‑CGCAAAACCCATCCCCGCTGGCGATCACATTGACGCGTTGTCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGa                                                                                                                                                                <  1:454353/251‑1 (MQ=255)
                                                                                             aggttattCTTTAACAGCTTCCAGTTTTATAAGAATAT‑CCCATTCTGCACCA‑ACAACGCCGCTGACTCCCTGCCCCAATCCCAGGCGCACCCGCGCTTTCTGTGCGACGTCTGCATTGTTCACCGCCTGCCCTTCTGCGTGGC‑CTCTAACCCCATCACCGCGGGCGATCACATTGACGCGTTGCCCCGCTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATcgc                                                                                                                                                            <  1:462256/244‑1 (MQ=255)
                                                                                              gggaaaTCTTAAACAGCTTACAGTTTTAAAAAAATATCCCCAT‑TTCCATAC‑ACAACCCCGCGGACTCCCTCTCCCGATCCCAGGCGCCCCCGCGCTTTCTGTGCGACGGCTGTTTTGTTCAGCGCTTGACCTTCTGCGTTGG‑CCCTAAACCCCCCCCCGCGGGCGATCACTTTGCCGCTTTGTCCCGCTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATcgcg                                                                                                                                                           <  2:509556/247‑1 (MQ=255)
                                                                                               ggTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCgcgc                                                                                                                                                          >  1:87581/1‑251 (MQ=255)
                                                                                                  aatttttaccatcTTACAGGTTAATAAAACTATTCCCAT‑CTGCA‑CACAAAACTCCTGTGACTTCCTGCCCCAATCCCATGCGCACCCGCGCATTCTGTGCGACGGTTGCTTTGTTCATCGCCTGCCCTTGGGTGTTGG‑CGCTAAACTCAGCCCCTCTGGCGTCCTCTTGGACGTGTTTTCCGTTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAg                                                                                                                                                        <  1:578531/238‑1 (MQ=255)
                                                                                                  aaTCTTTAACATCTTACAGGTTAATAAGAATATCCCCAT‑CTGCATCA‑ACAACGCCGTGGCCTCCCTGTCCCGATCCCTTGCGCCCCCGCGCATTCTGTGCGACGGCTGCATTTTCCAGCGCTTGCCCTTCTGCGTTGG‑CGCTAACCCCATCACCGCGGGCGATCACATTGACGCGTTGTCCCGCTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAg                                                                                                                                                        <  1:97243/250‑1 (MQ=255)
                                                                                                     ttttAACGGCTTACAGGTGTAAAAGAATTTCCCCAT‑CGGCACAA‑ACAACGCCGCTGCCTACCTGTCCCGATCCCATGCGCACCCGCTCATTCTGTGCGAGGGCTGCATTGTTCAGCGCGTGACTTTCGGCGTTGG‑CGCTAAACACATCACCGCGGGCGTTCACATTGACGCGTTGTCCCGCTTTTGCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCt                                                                                                                                                     <  2:125695/249‑1 (MQ=255)
                                                                                                       ttAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCACCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCCCCCGCGCATTCTGTGCGACGGCTGCATTTTTCAGCGCCTGCCCTTCGGCGTTGG‑CGCTAAACCCATCACCGCTGGGGATCACATTGACGCGTTGTCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCt                                                                                                                                                     <  1:100514/248‑1 (MQ=255)
                                                                                                                tgcAGGTTTAAAAAAATCTTCCCAT‑CTACACCA‑ACAACGCCGCTGTCTCCCTGTCCAGCTCCCAGGCGCCCCCGCTCATTCTGTGCGACGGCTGCATTTTTCGGCGCCTGCCGTTCGGCGTTGG‑CGCTAAACCCACCCCCGCTGGCGATCATATGGACGCGTTTTCCCGCTTTGACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaa                                                                                                                                          <  2:104921/248‑1 (MQ=255)
                                                                                                                 actgttttATAAGATTATCCCCAC‑CTGCATCA‑ACAACGCCGCTGCCTCCCTGTCCCGATACCATGCGCCCCCGCGCTTTTTGTGCGACGGCTGCATTGTTCAGCGCCTGCCCTTCTGCTTTGG‑CGCTAAACCCACCCCCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATcaac                                                                                                                                         <  1:60376/245‑1 (MQ=255)
                                                                                                                  cAGGTTTATAAAAATATCCCCTT‑CTACAACA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGAGACGGCTGAATTGTTCAGCGCCTGCCCTTCTGCGTTGG‑CGCTAACCCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGCACCTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCCGGCTAATGGCATcacca                                                                                                                                        >  2:579507/1‑247 (MQ=255)
                                                                                                                      tttATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttg                                                                                                                                   >  1:285051/1‑251 (MQ=255)
                                                                                                                        tataAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGACATCAACAAGttgat                                                                                                                                 >  1:387707/1‑251 (MQ=255)
                                                                                                                        aaaaaaaaTATCCCCTT‑CTGCATCA‑ACACCGCCGCGGCCTCCCTGCCCGGATACCATGCGCCCCCGCGCTTTCGGTGGGACGGTTGATTTGTTCGGCTCGTGACCTGCTGCGTTGC‑CGCTAACCCCATCCCCGCTGGGGATCACATTGAGGCTTTGCCCCGCTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgat                                                                                                                                 <  2:482910/245‑1 (MQ=255)
                                                                                                                         ataAGAATATTCCCAT‑CTGCATCA‑ACACCGCCGCTGACTCCCTGCCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGCCCTTCTGCGTTGG‑CGCTAAACCCATCCCCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGttgatt                                                                                                                                <  1:224973/251‑1 (MQ=255)
                                                                                                                            aaaaTATTCCCCT‑CTCCATCA‑TCACCGCCGCTTACTCCCTGCCCCGATCCCATGCCCCCCCCCGTATTCTGTGCGACGGCTGCATTGTTCAGCGCTTTCCTTTCTGGGTTGG‑CTTTAACCCCATCCCCGCGGGCTATCACATTGACGCTTGGTCCCGCTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTtgattgat                                                                                                                             <  1:262822/249‑1 (MQ=255)
                                                                                                                              aaTATCCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTTTCCCGTAACCATGCCCACCCGCGCATTCTGTGCGCGGGCTGCTTTGTTCAGCGCCTGCCCTTCTG‑GTTGGGCGCAAACCCCATCCCCGCTGGGGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGatat                                                                                                                           <  2:261803/251‑1 (MQ=255)
                                                                                                                                 atccccAT‑CTGCATCA‑ACAACGCCGCGGACTACCTGCCCCGATACCATGCCCCCCCGCGCATTCTGTGCGACGGCTGCTTTGTTCAGCGCTGGCCCTTCTGCGTTGG‑CTCTAAACCCATCCCCGCGGGCGATCAAATTGACGCGTTGCCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCa                                                                                                                        <  1:266100/248‑1 (MQ=255)
                                                                                                                                   ccccAT‑CTGCATCA‑ACAACGCCGCTGCCTCCCTGCCCCGATCCCGTGCCCCCCCGCGCATTCTGTGCGACGTCTGCATTGTTCAGCGCCTGCCCTTCTGCGTTGG‑CGATAACCCCATCACCGCTGGCGATCACATTGCCGCGTTGTCCCTTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATcc                                                                                                                         <  2:437590/247‑1 (MQ=255)
                                                                                                                                      cAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCAcc                                                                                                                   >  1:339873/1‑251 (MQ=255)
                                                                                                                                         tcttcATCA‑ACAACCCCGCTGACTCCCTGCCCCGATACCATGCCCCCCCGCGCATTCTGTGCGACGGCTGCTTTGTTCAGCGCCTGACCTTCGGCGTTGG‑CGCTAACCCCACCCCCGCGGGCTATCACATTGACGCTTTGCCCCTCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATTGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCg                                                                                                                  <  2:41090/246‑1 (MQ=255)
                                                                                                                                             gcttcaccaacGCCTCGGACTCCCTGCCCCGATCCCATGCGCCCCCGCGCATTCGGTGCGACGGCGGCATTGTTCAGCGCCGGACTTTCGGCTTTGG‑CGCTAACCCCACCCCCGCTGGCGATCACATTGACGCGTTTTCCCTTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTAc                                                                                                               <  2:339526/243‑1 (MQ=255)
                                                                                                                                                          cgcTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGt                                                                                                  >  1:532807/1‑250 (MQ=255)
                                                                                                                                                           gcTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAACCTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGAAGGGTAt                                                                                                >  1:562262/1‑251 (MQ=255)
                                                                                                                                                                                 aTGCGCACCCGCACATTCTGTGCGACGACTGCATAGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCCCCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCACGCCTGGCGCAACTGGGTTAACTGGCTAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGACAGGGTATCCAGCCGTCCGCGTTTCCGttt                                                                          >  2:359841/1‑251 (MQ=255)
                                                                                                                                                                                        cccGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTg                                                                   >  1:64969/1‑251 (MQ=255)
                                                                                                                                                                                         ccGCGTATTCTGGGCGAGGGCGGCATTGTCCAGCGCCTGCCCTTCTGCGTTGG‑CGCTAACCCCATCCCCGCGGGGGATCACATTGAGGCGTGGCCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGc                                                                  <  1:414754/251‑1 (MQ=255)
                                                                                                                                                                                              cATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGccagct                                                             >  2:159077/1‑251 (MQ=255)
                                                                                                                                                                                                            gCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTtcctcc                                               >  1:408983/1‑251 (MQ=255)
                                                                                                                                                                                                                cATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCcgc                                           >  2:375973/1‑251 (MQ=255)
                                                                                                                                                                                                                  ttGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCgcgc                                         >  2:522068/1‑251 (MQ=255)
                                                                                                                                                                                                                        ggcgctggcccTTCTGCTTGGG‑CGCAAACCCCACCACCGCTGGCGATCACATTGCCGCTTTGCCCCTTTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATggg                                   <  2:170573/242‑1 (MQ=255)
                                                                                                                                                                                                                        aGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATggg                                   >  1:585204/1‑251 (MQ=255)
                                                                                                                                                                                                                        aGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATggg                                   >  2:1204/1‑251 (MQ=255)
                                                                                                                                                                                                                            catGACCTTCTGCGTGGG‑CGAAAAACCCAGCCCCGCTGGCGATAACATTGACGCTTTGCCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATc                               <  1:114706/249‑1 (MQ=255)
                                                                                                                                                                                                                             cTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATcg                              >  1:132563/1‑251 (MQ=255)
                                                                                                                                                                                                                               aacctcctgctgggg‑CGCAAACCCCATCCCCCCTGGCGATCACATTGACGCGTTTCCCCTCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATcgcg                            <  1:2374/238‑1 (MQ=255)
                                                                                                                                                                                                                                     cTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTtcc                                                  >  2:538687/1‑223 (MQ=255)
                                                                                                                                                                                                                                     cTGCGTGGG‑CGCAAACCCCACCCCCGCGGGCGATCACATGGACGCGTGGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTtcc                                                  <  1:538687/223‑1 (MQ=255)
                                                                                                                                                                                                                                           tGG‑CGCAAACCCCATCCCCGCGGGCGATCACATTGCCGCGTTGTCCCGCTTTTCCCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATCGCGGCAACCACATa                 <  1:107739/250‑1 (MQ=255)
                                                                                                                                                                                                                                             g‑cgcAAAACCAAGCACCGGTGGAGATCACATGGAGGCGTTGTCCCGTTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCTGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATCGCGGCAACCACATAAtt              <  2:285051/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gcTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATCGCGGCAACCACATAATTTc            >  2:223272/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        cATCCCCGTTGGCGATCACATTGAGGGGTTGCCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATCGCGGCAACCACATAATTTCCTGTGGcc    <  2:550181/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          tCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATCGCGGCAACCACATAATTTCCTGTGGCCTg  >  2:157489/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          tCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATCGCGGCAACCACATAATTACCTGTGGCCTg  >  1:414564/1‑251 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
TACAGGCTTCAGAGGCGAAGTGCGATCAATACTCATGGTTTATTCCTCATTGAGGGCGCTTTTATCATGTGTTGCTTATTTATCGGCAAGGGACGGGTAATCTTTAACAGCTTACAGGTTTATAAGAATATTCCCAT‑CTGCATCA‑ACAACGCCGCTGACTACCTGTCCCGATACCATGCGCACCCGCGCATTCTGTGCGACGGCTGCATTGTTCAGCGCCTGACCTTCTGCGTTGG‑CGCTAAACCCATCACCGCTGGCGATCACATTGACGCGTTGTCCCGCTTTTACCCGCCATGCCTGGCGAAACTGGGTTAACTGGATAGGTTGATCGGGTGATAGATCGCGCAGGCTAATGGCATCAACAAGTTGATTGATATCCAGCACCGTACGCGGTGGCAGGGTATCCAGCCGTCCGCGTTTCAGTTTGACATTGCCAGCTTCCAGCTTTCCTCCCCGCGCAATGGGCATCGCGGCAACCACATAATTTCCTGTGGCCTG  >  NC_000913/1130093‑1130590

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: ‑