Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,8370AG30.1% 160.7 / 16.5 73I572T (ATC→ACC) fhuEferric coprogen/ferric rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (4/47);  new base G (22/0);  total (26/47)
Fisher's exact test for biased strand distribution p-value = 5.83e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTTT  >  NC_000913/1159589‑1160085
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
gCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAAGGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGAt                                                                                                                                                                                                                                                            <  1:704125/250‑1 (MQ=255)
 ccAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAg                                                                                                                                                                                                                                                          <  1:63479/251‑1 (MQ=255)
      tAACCCCACCGCCGACCTGCAACTCTGGCGTGACAGGCTACCGATAGCGGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCTAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGGCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGttgt                                                                                                                                                                                                                                                      <  2:722343/250‑1 (MQ=255)
      tAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGTCCTTGCTTCCGGGGGTAggggg                                                                                                                                                                                                                                                      >  2:649035/1‑247 (MQ=255)
      tAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGgtgt                                                                                                                                                                                                                                                      <  1:807038/250‑1 (MQ=255)
        acacCACCGCCGACCGTCAACTCTGGCATTACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCGTTGCTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTAcc                                                                                                                                                                                                                                                   <  2:382551/251‑1 (MQ=255)
               cgccgACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGAc                                                                                                                                                                                                                                           <  2:852364/251‑1 (MQ=255)
                   gACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCTCCGCTTCTGTCGGGGGGAGGGGTACCGGTGGGCTggg                                                                                                                                                                                                                                       >  2:42538/1‑249 (MQ=255)
                         cAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTTGACTGAGCGACa                                                                                                                                                                                                                                 <  1:690470/250‑1 (MQ=255)
                           aCTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTa                                                                                                                                                                                                                              <  1:430304/251‑1 (MQ=255)
                                 gCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGc                                                                                                                                                                                                                        <  2:581115/251‑1 (MQ=255)
                                   aTGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCGGTGGGCTGAGCGGCATTATCCTGCTc                                                                                                                                                                                                                      >  1:286850/1‑251 (MQ=255)
                                     gACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCaa                                                                                                                                                                                                                    <  2:453676/251‑1 (MQ=255)
                                      aCAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTTCCCTTTCTTTCGGGGGTAGGTGTACCggaggatcgagcgcccgcttccggctctcc                                                                                                                                                                                                                    >  2:388334/1‑227 (MQ=255)
                                      aCAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAAt                                                                                                                                                                                                                   <  2:445556/251‑1 (MQ=255)
                                       cAGGAAACCGATAGCTGGTTAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATa                                                                                                                                                                                                                  <  2:1272/251‑1 (MQ=255)
                                                    gCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGGACCGGTGGGCTGGGCGGCATTTTcccaggtttgaaaagcggcgggtt                                                                                                                                                                                                     >  2:242614/1‑229 (MQ=255)
                                                    gCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCt                                                                                                                                                                                                                          <  1:242614/230‑1 (MQ=255)
                                                    gCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTa                                                                                                                                                                                                    <  1:813671/251‑1 (MQ=255)
                                                       ggTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCATTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACAGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAAc                                                                                                                                                                                                  <  1:488553/250‑1 (MQ=255)
                                                        gTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGGCGCCGGGGGTAGGGGTACCCGGGGGCGGGGCGGCAATATCCTGCTCATTACGGGAGATggggcacgg                                                                                                                                                                                                 >  2:183694/1‑244 (MQ=255)
                                                                  tAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCTTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCat                                                                                                                                                                                       <  2:50357/250‑2 (MQ=255)
                                                                                ggCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCa                                                                                                                                                                        <  2:286850/251‑1 (MQ=255)
                                                                                         ggATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGAttt                                                                                                                                                               <  1:540365/251‑1 (MQ=255)
                                                                                             tAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCCGTGGTGCCGGGGGTAGGGGTACCGGGGGGCTGAGCGGCATTTTCCCGCTCAATACGGGAGATGGGCTACGTGGGGGTCAGACGGCTATTCATCCAGTCCGAtttt                                                                                                                                                              >  1:762176/1‑247 (MQ=255)
                                                                                              aaCGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTTCCGGGGGTAGGGGTAACGGGGGAATGAGCGGCATTATTCTGCTCCATACGGAAGATGGCGTACGTGGAGGGCAGACGGCTATTTATCCAGTCCGATTTCAGAc                                                                                                                                                          >  1:20786/1‑251 (MQ=255)
                                                                                               aCGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGAcc                                                                                                                                                         <  1:238136/251‑1 (MQ=255)
                                                                                                  gCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGTGTACCGGTGGAATGAGCGGCATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGGCAGACGGCTATTCATCCCGTCCGAttttcgaacc                                                                                                                                                        >  1:528101/1‑242 (MQ=255)
                                                                                                     tttCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGTGTACCGGTGGGCTGAGCGACATTATCCTGCTCCATACGGAAGATGGGTAACGTGGGGGTCAGACGGCTATTCATCCCGTCCGCtttttgaacccgctc                                                                                                                                                   >  1:212699/1‑239 (MQ=255)
                                                                                                      ttCCTTCTTTATCCGCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTc                                                                                                                                                   <  2:670741/250‑1 (MQ=255)
                                                                                                       tCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCg                                                                                                                                                  <  1:532736/250‑1 (MQ=255)
                                                                                                           tCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgt                                                                                                                                              <  1:840415/250‑1 (MQ=255)
                                                                                                            cGTTATCCTCTGCAATATATCGCGGTGCGCCACATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgttg                                                                                                                                            <  2:327702/251‑1 (MQ=255)
                                                                                                                 tCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGCTGCTGCCGGGGGTAGGTGTACCGGGGGGCTGAGGGACACTATCCTGCTCCATACGGAAGATGGGTAACGTGGGGGGCAGAAGGCTACTCATCCACTCCGATTTTAGACCCAGCTCGTAGTTGTTAcc                                                                                                                                       >  2:98564/1‑251 (MQ=255)
                                                                                                                        cAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCCTTGCTGGCGGGGGCAGGTGGACCGGTGGACTGAGCGAAATTATCCTTCTCAATACGGAAGGTGGCTAAAGTGGGGGGCAGACGGCTATTCCTCCAGGCCGCTTTTAGACCCCGCTCGGAGTTGTTACCGGGGAtt                                                                                                                                >  1:568868/1‑251 (MQ=255)
                                                                                                                        agatatTGCGGGTTGCGCCAAAGGGCAGCTGCCAGTTGGCGGTAATTGCGCCGTGGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCGCCGCTTTATAGGCGGTGTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGAtt                                                                                                                                <  2:298894/249‑1 (MQ=255)
                                                                                                                          atatAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTgg                                                                                                                              <  1:601731/251‑1 (MQ=255)
                                                                                                                            ataGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGAAGGGGTACCGGGGGACTGAGCGACATTATTCTGCTCAATACGGGAGATGGCTAACGGGGGGGGCAGAAGGCTTATCATCCAGTCCGATTTCCGACCCAACTCGGAGATGTTACCGGTGATTgg                                                                                                                              >  1:775641/1‑249 (MQ=255)
                                                                                                                            ataGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAg                                                                                                                            <  1:649035/251‑1 (MQ=255)
                                                                                                                            ataGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGCTTCTGCCGGGGATAGGGGGGCCCGGGGGGTGGGCGGACTTATACTGCGCAATACAGAAGGAGGCTAAAGTGGGGGGCAGACGGCGTTTTATCCAGTCCGATTTCAAACCCAGCCCGCAGGTGTTAcccggggttgggg                                                                                                                            >  2:243257/1‑240 (MQ=255)
                                                                                                                                            aaTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGa                                                                                                            <  1:570665/251‑1 (MQ=255)
                                                                                                                                                  aGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGGTTCTGCCGGGGGTAGGGGTACCGGTGGACGGGGGGGCATTATCCTTCTCAATACGGAAGGTGGGCAACGGGGGGGGCAGACGGCTATTCATCCAGGCCGAGTTCAGACCCAGCTCGTAGGTGTTTCCCGGGATTTGAGCCCGATATTTGCCTGGACttgc                                                                                                      >  2:154658/1‑248 (MQ=255)
                                                                                                                                                      gCCAGTTGTCGGGAACGGCGCCGCTGACTGCAAATTCCACCACTTTACTGACTGTCCCAGCCTGCGTTTGATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCAGCATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGAGCGCTATTCATCCAGTCCGAGTTCAGACCCAGCACGTAGTTGTGACCGGTGATTGGAGCCAGATATTTACCTGAACTGTcgc                                                                                                    <  2:581386/249‑3 (MQ=255)
                                                                                                                                                       ccAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACGGAAGAAGGCTAAAGGGGGGGGCAGACGGCTATTCATCCAAGCCGATTTCAGAACCAACTCGTAGTTGTTAACGGGGATTTGGGCCAGATATTTTCCTggaccttctcgat                                                                                                 >  1:70830/1‑239 (MQ=255)
                                                                                                                                                       ccAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTCCCGTTGCTGCCGGGGGTAGGTGTACCCGGGGGATGGGGGACATTATCCTGCTCCATACGGAAGATGGGTAACGTGGGGGGCAGAAGGCTATTCATCCAGTCCGATTTTAGACCCAGCTCGTAGTTGTTACCGGGGGTTGGAGCCAGATATTTGCCTGCACTGTCACGAt                                                                                                 >  1:443340/1‑251 (MQ=255)
                                                                                                                                                             gTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCAttt                                                                                            <  1:304051/250‑1 (MQ=255)
                                                                                                                                                                 gTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCCCCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                                                                                      <  2:16505/251‑1 (MQ=255)
                                                                                                                                                                   aaTTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGCTTCTGCCGGGGGTAGGGGTACCGGGGGACCTAGCCACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGGGGTCAGACGGCTATTCATCCACTCCGATTTCAGACCCCGCTCGTTGTTGTTACCGGGGGTTGGGGCCCGATATTTGCCTGAACTGTCACGATCATTTTGCGGCt                                                                                     >  1:324783/1‑251 (MQ=255)
                                                                                                                                                                      tGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGTGGGGGGCAGACGGCTATTTATCCAGTCCGATTTCAGACCCAGCTCGTAGATGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGCTCCTTTTGCGGCTggg                                                                                  >  1:799450/1‑250 (MQ=255)
                                                                                                                                                                        cgcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATAtt                                                                                                                   <  1:484740/216‑1 (MQ=255)
                                                                                                                                                                         gcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAAGTGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGTGGGGGGCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTTCCCGTGATTTGAGCCCGATATTTGCCTGAAATGTCACGATCATTTTTTGGCTgggaa                                                                                >  1:821892/1‑247 (MQ=255)
                                                                                                                                                                         gcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCCGGGGTAGGTGTACCCGTGGACTGAGCGGCATTATCCTGCTTAATACGGGAGATGGCTAACGTGGTGGTCAGACGGCTATTTATCCAGGCCGATTTCCGACCCCGCTCGTAGTTGTTTCCGGTGGTTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTTCGGCTgggaa                                                                                >  1:277370/1‑247 (MQ=255)
                                                                                                                                                                          ccGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTGCCGGGGGTAGGGGTACCGGTGGGCTGAGGGACATTATCCTGCTCAATAAGGAAGATGGCTAACGTGGGGGGCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCGCGTAGTTGTTACCGGTGATTTGGGCCAGATATTTGCCTGGACTTTCACGATCATTTTGCGGCTGGGAAAt                                                                              >  1:79126/1‑251 (MQ=255)
                                                                                                                                                                               gAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAgg                                                                         <  1:342825/251‑1 (MQ=255)
                                                                                                                                                                                agtgcaaattcatcccttgtACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGCCATTAGCCTGCTCAAGACGGAAGATGGCTAACGTGGGGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGt                                                                        <  2:343673/232‑1 (MQ=255)
                                                                                                                                                                                   tCAAATGCCACCCCTTTACTGACTGGCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGtat                                                                      <  2:196079/250‑1 (MQ=255)
                                                                                                                                                                                    cAAATTCCGCCCCTTTACTGACCGTCCCATCCTCCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATCTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGACTGAACTGTCACGATCATTTTGCGGCTGTAAAATAGAGGTATAg                                                                    <  2:179270/251‑1 (MQ=255)
                                                                                                                                                                                    cAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGtata                                                                     <  1:42538/250‑1 (MQ=255)
                                                                                                                                                                                    cAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGGTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAg                                                                    <  2:832510/251‑1 (MQ=255)
                                                                                                                                                                                     aaaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGc                                                                   <  1:8172/251‑1 (MQ=255)
                                                                                                                                                                                             ccccTTTACTGACTGTACCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAg                                                           <  2:568868/251‑1 (MQ=255)
                                                                                                                                                                                              cccTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATGTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAgg                                                          <  1:747285/251‑1 (MQ=255)
                                                                                                                                                                                                                aTCCACCGCTTTATAGGCGGTTTCGCCGGTGCTGCCGGGGGTAGGGGGACCGGTGGGCTGGGCGACATTATTCTGGTCAATAAGGAAGATGGCTAACGTGGTGGTGAGACGGCTATTCATCCAGGCCGATTTTAGACCCAGCTCGTAGGTGTTACCGGGGATTGGAGGCAGATATTTGCCTGAACTGTCACGATCATTTTGTGGTGGGGAAATTGAGGGATAGATGGCGGAGGGCGGACAATTGTCATTg                                         >  2:288954/1‑250 (MQ=255)
                                                                                                                                                                                                                     ccGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCa                                   <  1:183694/251‑1 (MQ=255)
                                                                                                                                                                                                                           tataGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATTTCAAacac                              <  1:403262/250‑1 (MQ=255)
                                                                                                                                                                                                                              aGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACaccag                           <  1:480368/250‑1 (MQ=255)
                                                                                                                                                                                                                                gCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGGGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACaccagacc                        >  2:133623/1‑251 (MQ=255)
                                                                                                                                                                                                                                gCGGATTCGCCGTTGCTGCCGGGGATCGGTGTATCGGTGTACTGAGCGTCATTAGCCTGTTCTATACGGAAGATCGCTAACGTTGTGTTCAGGCGGCTATTCATCCAGTCCGCTTGCAGACCCAGCTGGTAGTTGTTACCGGTGATTGGAGCCACATATTTGCATCAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTAAAACaccaaacc                        <  2:449669/251‑1 (MQ=255)
                                                                                                                                                                                                                                 cGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACaccagacca                       <  1:357899/251‑1 (MQ=255)
                                                                                                                                                                                                                                  ggTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGGAGGTCGACCAGTTGTCATTGATGTCAAACaccagacc                        >  2:687613/1‑249 (MQ=255)
                                                                                                                                                                                                                                         ccGTGGCTGCCGGGGATAGGTGTACCGGTGGAAGGAGAGACATTATCCTGCTCAATACGGAAGAGGGCTAACGTGGTGGTCAGACGGCTTTTCACCCAGTACTATTTCAGACCCAGCTCGTAGTTGTTACCGGGGATTTGTGCCAGATATTTGTCTGAACTGTCACGATCATTTTTCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGTCCAGTTGTCATTGATGTCAAACTCCAGACCAGCGTa                  <  2:537487/248‑1 (MQ=255)
                                                                                                                                                                                                                                          cGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATAGTTGCCTGAACTGGCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCCACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGc              <  1:683541/251‑1 (MQ=255)
                                                                                                                                                                                                                                                      ggATAGGTGTACCGTTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGTCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGGAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGGGGTGTGGttt  >  2:416392/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                           
GCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTTT  >  NC_000913/1159589‑1160085

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: