Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,917,0880GC22.9% 171.7 / 11.8 71A678P (GCT→CCT) barAsensor histidine kinase BarA
Reads supporting (aligned to +/- strand):  ref base G (8/46);  new base C (16/0);  total (24/47)
Fisher's exact test for biased strand distribution p-value = 2.97e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGATTAGCGAAAGCGGTATCGATGACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATGCAATGGCCGGGCAAAAAGAGAAGCTGCTTGG  >  NC_000913/2916848‑2917338
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                  
cGATTAGCGAAAGCGGTATCGATGACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑GTGACAAACCCCCTCACCTGa                                                                                                                                                                                                                                                           >  1:760613/1‑250 (MQ=255)
   ttAGCGAAAGCGGTATCGATGACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGGTG‑TTTACAACCCCCCTCAACTGAAACt                                                                                                                                                                                                                                                       >  1:259056/1‑251 (MQ=255)
     aGCGAAAGCGGTACCGATGACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGGAAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTTCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTa                                                                                                                                                                                                                                                     <  2:679379/251‑1 (MQ=255)
        gAAAGCGGTATCGATGACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCGACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCg                                                                                                                                                                                                                                                  <  1:433210/251‑1 (MQ=255)
             cGGTATCGATGACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTT‑ATTAAAACCCCCCCAACCTGAAAATTATCAgcgcc                                                                                                                                                                                                                                             >  1:144861/1‑250 (MQ=255)
               gTATCGATGACCGATTTCCTGATGGTGGCACTTTCTTGCTATGCACAAGTCAATGCTGAAAAACTCAAGCAAGTTGGTAGCGGCGCGTGTCTGCTGAAACCATTAACTCCTACACGCCTGTGGCCTGCCCTGACGGAAGTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGGAAGCTGGCAATGACAGTTATGGCGGTTG‑ATGACAACCACGCTAACCTGAAACTTATCGGCGCAtt                                                                                                                                                                                                                                           <  1:779734/251‑1 (MQ=255)
                   cGATGACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑t                                                                                                                                                                                                                                       <  1:727257/250‑1 (MQ=255)
                       gACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGAGGTTG‑ATGACAAACCCGCCAAACCGAAAACTATCGGGGCATTTC‑Tggg                                                                                                                                                                                                                                    >  1:186772/1‑248 (MQ=255)
                       gACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCAACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGACAGTAAGCTGGCAATGCCAGTCATGGCGGTTG‑TTGAAAACCCCCCCCACCTGaaaattttgcggccgtttgtggg                                                                                                                                                                                                                                     >  2:76715/1‑229 (MQ=255)
                        aCCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGa                                                                                                                                                                                                                                 <  1:767077/251‑1 (MQ=255)
                           gATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATTACAACCCCCCTCAACTGAAACTTATTGGGGCAATGC‑TGGGAGatat                                                                                                                                                                                                                              >  1:88664/1‑251 (MQ=255)
                             tttCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATg                                                                                                                                                                                                                             <  2:88664/250‑1 (MQ=255)
                                   gactcTGGCACTTCCTTGCCATGCACAAGTCAATTCATAAAAACGCAATCCGGATGGTAATGGAGCGTGACTGCGGAAACCAATAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTTCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCTCGCTAACCTGAAACTTATCGGCGCATTGA‑TGGAAGATATGGTGC‑aa                                                                                                                                                                                                                      <  2:757599/246‑1 (MQ=255)
                                   gATGCTGTCACTTCCTTGCCATGCACAAGACAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AAc                                                                                                                                                                                                                     <  2:149460/251‑1 (MQ=255)
                                    aTGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGAAAACCCCGCTAACCTGAAACTTATCGGCGCATTGCATGGaaaaatggggcaac                                                                                                                                                                                                                       >  1:831878/1‑238 (MQ=255)
                                      gcgggCACTTCCTTGCCATGCACAAGTCAATGCGGAAAAACTCAAGCAAGAGGGTATCGGCGCGTGTCTGCTGAAACCAGTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTGATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAt‑g                                                                                                                                                                                                                 <  2:103857/248‑1 (MQ=255)
                                         ggCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑Tgg                                                                                                                                                                                                             <  1:788545/251‑1 (MQ=255)
                                                       aagaacaAGTCAATGCTGAAAAACTCAATCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACTCCTGTTGCCTGCCCTTACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATTAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGc                                                                                                                                                                                             <  2:252842/247‑1 (MQ=255)
                                                             aaGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCAt                                                                                                                                                                                       <  2:73519/251‑1 (MQ=255)
                                                              aGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTAATCACAAAAAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑GTGACAACCCCCCCAACCCGAAAATTATCGGCGCAATGC‑TGGGAAATTAGGTGG‑AACAT‑G‑TGG‑AAACTTGC‑GAGAGAGGGCATc                                                                                                                                                                                      >  1:127268/1‑251 (MQ=255)
                                                               gTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCCCTAACCTGAAACTTATTGGCGCATTTC‑TTGAAGGAATGTGGC‑AAAAT‑G‑TGG‑GACttttccataacggggatc                                                                                                                                                                                       >  1:11684/1‑235 (MQ=255)
                                                                 cAATGCTGAAAAACTCAAGCAAGAGGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAgg                                                                                                                                                                                   <  2:216949/251‑1 (MQ=255)
                                                                        gAAAAACTCCACCAAGATGGTATCGGCGCGTGTCTGCAGAAACCATTAACCCCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTg                                                                                                                                                                             <  2:162121/250‑1 (MQ=255)
                                                                         aaaaaCTCAAGCAAGAGGGTATCGGCGCGTGTCTGCGGAAACCATTAACACCTACACGCCTGGTGCCGGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAGCGCCGCTAAGCTGAAACGTAGCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGaa                                                                                                                                                                           <  1:90557/251‑1 (MQ=255)
                                                                                caatcaagATGGTATCGGCCCGTGGCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTTCCCTGACGGAATTTGGTCATCACAAACAAAACACGCTTGTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCAGTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGc                                                                                                                                                                      <  2:63162/245‑1 (MQ=255)
                                                                                caagcaagATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGc                                                                                                                                                                      <  1:521291/249‑1 (MQ=255)
                                                                                caagcaagATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCCCTAAACCGAAACTTATCGGCGCATTGC‑TGGGAGATATGGGGC‑AAACT‑G‑TGG‑GACTTTGC‑GATAGCGGGCATCAGGCGGGTGGACGGGCgg                                                                                                                                                                    >  1:120422/1‑250 (MQ=255)
                                                                                     aagATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCa                                                                                                                                                                                        <  2:330371/226‑1 (MQ=255)
                                                                                          ggTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑GTGACAACCCCCCTAACCTGAAAATTATCGGGGCCTTGC‑TGGGAGATATGGTGC‑AAAATTG‑GGG‑AA‑TTTGC‑GATAGCGGGCATAAGGGGGGTGTACggggggaaaacattgc                                                                                                                                                          >  1:138478/1‑238 (MQ=255)
                                                                                             aTCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGt                                                                                                                                                        <  1:468499/250‑1 (MQ=255)
                                                                                                ggCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATc                                                                                                                                                                                                                                                   <  1:445128/162‑1 (MQ=255)
                                                                                                  cgccggtCTGCTGAACCCATAAACACCTACACGCCTGTTGCCTGCCCGGACGGAATTTTTTCAACACAAACAAAACACGCCTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGATAACCCCGCTAGCCTGAAACTTATCGGCGCATGGC‑TGGAAGATAGGGTGC‑AACAT‑G‑TGG‑AACGGTGC‑GATAGCGGGCATCAGGCGGTTGAACGGTCGAAACTGATGCCGTTCGAtt                                                                                                                                                  <  2:339586/246‑1 (MQ=255)
                                                                                                    cgttcctgctgAAACCAGTAACTCCCTCACGCCTGTCTCCTGCCCTGTCGGAATTTTGTCATCACAAACCAAACACGCTTTTGCCTGTATCCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACATCCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGAGATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGAGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGAtt                                                                                                                                                  <  2:482839/244‑1 (MQ=255)
                                                                                                       gtCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGAt                                                                                                                                              <  1:368630/250‑1 (MQ=255)
                                                                                                          tgctgAAACCATTAACACCTACACGCCTGTTGCCTACCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGa                                                                                                                                                                                                                                    <  1:103282/166‑1 (MQ=255)
                                                                                                           gcttaaaCCAATAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTGGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCTAAACAGATGCCGTTCGATTTGATCTTa                                                                                                                                          <  2:738094/246‑1 (MQ=255)
                                                                                                           gctgAAACCATTAACACCTACACGCCGGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTa                                                                                                                                          <  2:73295/250‑1 (MQ=255)
                                                                                                              gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑GTGGCAACCCCCCTAACCTGAAACTTATCGGCGCCTTTG‑TGGAAGGTATTGTGA‑AACAA‑G‑TGGGAACTTTG‑‑GGTTGCGGGGATCAGGGGGGTGTAAGGGCGCAACAGaagacctttggttttattttaatgg                                                                                                                                      >  1:861415/1‑227 (MQ=255)
                                                                                                              gAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAAt                                                                                                                                        <  1:62070/249‑1 (MQ=255)
                                                                                                                     ttAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCCCCAACCTGAAAATTATCGGCGCATTGC‑TGGAAGATATGGGGC‑AACAT‑GGTGG‑A‑CTTTGC‑GATAGCGGGCATCAGGCGGGTTAACGGGCGGAACAGATGCCGTTCGATTTTGTCTTTATGGATAtta                                                                                                                                >  1:738918/1‑249 (MQ=255)
                                                                                                                          acCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGGCAACCCCCCCAACCTGAAACTTATTGGCGCATTTC‑TGGAAGATATGGTGCCAAAAT‑‑‑TGG‑GACTTTTC‑GTTAGCGGGCATCCAGCGGTTGAACGGGCGGAACAGATGCCCTTTGATTTTATCTTAAATGATAttttaattc                                                                                                                          >  1:717047/1‑244 (MQ=255)
                                                                                                                          acCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACg                                                                                                                                                                         <  2:823600/204‑1 (MQ=255)
                                                                                                                          acCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAAACCCCCTAACCCTAAACTTATCGGCGCATTTC‑TGGGAGATATTGGGC‑AACAT‑G‑TGG‑GACTTTTC‑GATAGCGGGGCTCAAGCGGTTGAACg                                                                                                                                                                         >  1:823600/1‑204 (MQ=255)
                                                                                                                             cacacGCCTGTTGCTTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCgg                                                                                                                       <  2:711084/250‑1 (MQ=255)
                                                                                                                                    cTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTTTAT‑ACAACCCCCATAACCCTAAAATTATCGGCGCACTGC‑TTGAATATATGGGGC‑AACAT‑A‑TGG‑AACATTTC‑GCTTTCGGGCATAAGGCGGTTTAAAGGGCGAAAGAGATGCCGGTCGATATGTTCTTTATTGGTATTCacatcacggtcagagg                                                                                                                >  2:288087/1‑236 (MQ=255)
                                                                                                                                    cTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGATG‑ATGAAAAACCCGCTAACCTGAAACTTATCGGCGGATTGC‑TTGAAGAAATGGGGC‑AAAAT‑G‑TGG‑AACTTTTC‑GATAGGGGGGATCAGGCGGGTGAAAGGGCGAAAAAGATGACGTTTGGTTTTATCTTAATGGGTATTCAAATGACGGACa‑tggg                                                                                                               >  1:601296/1‑250 (MQ=255)
                                                                                                                                     tGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑TTGACAACCCCCCCAACCTGCAACTTATTGGCGCAATGC‑TGGAAGGAATGGTGC‑AACAT‑T‑TGG‑GACTTTGC‑GCTAACCGGCATCCGGCGGGTGGAAGGGCGGAACAAATGCCGGTCGGTTTGCTCTTAATAGATATTTAAATGCCGGACa‑aggat                                                                                                              >  1:518208/1‑251 (MQ=255)
                                                                                                                                             gccCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTAACGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCg                                                                                                      <  1:332217/251‑1 (MQ=255)
                                                                                                                                                 tGACGGAATTTTTTCATCACAAACAAAACTCGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAATCTTATCGGCGCATTGC‑TGTAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TTTATGGCATTCggg                                                                                                    <  2:747442/249‑1 (MQ=255)
                                                                                                                                                  gACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCt                                                                                                 <  1:458455/251‑1 (MQ=255)
                                                                                                                                                      gAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGa                                                                                             <  1:302902/251‑1 (MQ=255)
                                                                                                                                                           ttGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTc                                                                                         <  1:5235/250‑1 (MQ=255)
                                                                                                                                                                atcaCAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATccac                                                                                   <  2:120422/251‑1 (MQ=255)
                                                                                                                                                                     aaacaaaacaCGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATGCGGGCCTGCGAACTCATCCACCAA‑Ct                                                                             <  2:760613/251‑1 (MQ=255)
                                                                                                                                                                        caaaacaCGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑GTTACAACCCCCCTCACCCGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AGCTTTGCCAATA‑CGGGCATCAGGCGGTTGAACAGGCGAAACAGAGGCCGTGCGATTTGTGCTTTATGGGAATTCAAAAGCCGGACA‑TTGATGGCATTCGTGCCCGCGCACTTAACTACCAA‑CTccc                                                                          >  1:201378/1‑251 (MQ=255)
                                                                                                                                                                                    gtGCCTGTAACCGAGGAAAGTACGCTTGCAATGACAGTCATGGTGGTTG‑ATGACATCCCCGCTACCCTTAAACTTATCGGCTCATTGC‑TGGAAGATATGGTGC‑AACTT‑G‑TGG‑AACTTTGC‑GATTGCGGGCATCAGGCGGTTTAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGTATATTCAAACGCCTGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCAGCAGcaac                                                               <  2:738176/249‑1 (MQ=255)
                                                                                                                                                                                       ccTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAAc                                                           <  1:239839/251‑1 (MQ=255)
                                                                                                                                                                                                 aTGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATACGCGCCTTAC‑TGGAAGATATGGTGC‑AACAA‑T‑TGG‑AACTTTGC‑GATAGAGGGCATCCAGCGGGTGAAAGGAGGAAAAAAAGGACGTTCTATTTGAGTTTTATGGATATTAAAATTCACGACAATGG‑TGGGATTCGGGCCTGCGGACTCATCCACCAAACTccccccaagcacaaaaagccgggtttc                                                  >  1:802033/1‑227 (MQ=255)
                                                                                                                                                                                                   gAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATcgc                                               <  1:165076/251‑1 (MQ=255)
                                                                                                                                                                                                                cAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCAt                                  <  1:76715/251‑1 (MQ=255)
                                                                                                                                                                                                                 aaTGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATTGATATTCAAATGCCGGACA‑TGGATGGCAGTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATg                                 <  1:849243/251‑1 (MQ=255)
                                                                                                                                                                                                                  aTGACAGTCATGGCGGCTG‑ATGACAACCCCGCTAACCTGAAACTTAGCGTCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATGc                                <  2:738918/251‑1 (MQ=255)
                                                                                                                                                                                                                           aTGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGc                                                                                               <  1:292319/180‑1 (MQ=255)
                                                                                                                                                                                                                               cGGTTG‑ATGACAACCCCTCTAACCGGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATGCAATGGCCGGGCaa                   <  2:186772/251‑1 (MQ=255)
                                                                                                                                                                                                                                ggTTG‑ATGACATCCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GGTAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCAATCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGc                                                                  <  2:534783/203‑1 (MQ=255)
                                                                                                                                                                                                                                ggTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGc                                                                  >  1:534783/1‑203 (MQ=255)
                                                                                                                                                                                                                                   tg‑atgaCAACCCCGCTAACCTGAAACCTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑GGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATGCAATGGCCGGGCaaaaa                <  2:259056/250‑1 (MQ=255)
                                                                                                                                                                                                                                       tgaCAAACCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACTCCGGTTATCGCGGTAACGGCGCATGCAATGGCCGGGCAAAAagcg             >  1:417618/1‑248 (MQ=255)
                                                                                                                                                                                                                                           aaCCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCCTGCAATGGCCGGGCAAAAAGAGAAgct        >  2:854023/1‑251 (MQ=255)
                                                                                                                                                                                                                                               ccGCTCACCTGAATCTTGTCGGCGCATTGC‑TGGAAGATATGGTGC‑AAAAT‑G‑TGG‑AAATTTGC‑GATAGCGGGCAGCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCAGTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATTCAATGGTCGGGCAAAAAGAGAAGCTGCtt    <  2:717047/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATGCAATGGCCGGGCAAAAAGAGAAGCTGCTTg   >  1:266515/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 gCTAACCTGAAACTTAACGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATGCAATGGCCGGGCAAAAAGAGAAGCTGCTTgg  <  2:619289/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                                  
CGATTAGCGAAAGCGGTATCGATGACCGATTTCCTGATGCTGGCACTTCCTTGCCATGCACAAGTCAATGCTGAAAAACTCAAGCAAGATGGTATCGGCGCGTGTCTGCTGAAACCATTAACACCTACACGCCTGTTGCCTGCCCTGACGGAATTTTGTCATCACAAACAAAACACGCTTTTGCCTGTAACCGATGAAAGTAAGCTGGCAATGACAGTCATGGCGGTTG‑ATGACAACCCCGCTAACCTGAAACTTATCGGCGCATTGC‑TGGAAGATATGGTGC‑AACAT‑G‑TGG‑AACTTTGC‑GATAGCGGGCATCAGGCGGTTGAACGGGCGAAACAGATGCCGTTCGATTTGATCTTAATGGATATTCAAATGCCGGACA‑TGGATGGCATTCGGGCCTGCGAACTCATCCACCAA‑CTCCCGCATCAGCAACAAACGCCGGTTATCGCGGTAACGGCGCATGCAATGGCCGGGCAAAAAGAGAAGCTGCTTGG  >  NC_000913/2916848‑2917338

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