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

GCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTT  >  NC_000913/1159589‑1160084
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
gCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATa                                                                                                                                                                                                                                                        <  1:247466/251‑1 (MQ=255)
    gTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGgtg                                                                                                                                                                                                                                                    >  1:56650/1‑251 (MQ=255)
         caccacCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCg                                                                                                                                                                                                                                               <  1:113817/251‑1 (MQ=255)
             accGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTCGGTGTACCGGTgg                                                                                                                                                                                                                                           >  1:384460/1‑251 (MQ=255)
              ccgccgACCGTCAACTCGGTCAGGACTGGCAACTGGTAGCTCGTGGAAGTTTTAGTCGTGGTGCGTGGCAGATTTGGGTTAACGGTGTTTCCTGCGTTATCCTCGGCAATATAGCGCGTTGCTCCAAATGTCAGCTGCCAGTTTTCGGTTAGTGCTCTTTTTAGTTGAAATTCCACCCCTTTACAGACTGTCGCATCCTCCTCGTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGa                                                                                                                                                                                                                                          <  2:240954/251‑1 (MQ=255)
                gccgACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACt                                                                                                                                                                                                                                        <  2:26305/251‑1 (MQ=255)
                           aCTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTTCTGTCGGGGATAgggggcccgggggaccgcgcgccattc                                                                                                                                                                                                                             >  1:233824/1‑226 (MQ=255)
                                ggCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATc                                                                                                                                                                                                                           <  1:363457/248‑1 (MQ=255)
                                ggCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTg                                                                                                                                                                                                                        <  2:86485/251‑1 (MQ=255)
                                   aTGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCCGTTGTGCCGGGGATAGGGGTACCGGTGGACGGGGTGGCATTATCCTggcc                                                                                                                                                                                                                     >  2:362515/1‑248 (MQ=255)
                                         ggCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTTCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGGCTGGGGGGCATTATCCTTCTCCATACg                                                                                                                                                                                                               >  2:120843/1‑251 (MQ=255)
                                             aCCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTTACCTGCTCAATACggg                                                                                                                                                                                                             >  1:70329/1‑248 (MQ=255)
                                              ccGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGa                                                                                                                                                                                                          <  1:120843/251‑1 (MQ=255)
                                                 aTAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATgg                                                                                                                                                                                                       <  1:235221/251‑1 (MQ=255)
                                                              aTTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGGGGGCTGGGCGACATTATCCCGCTCAAAACGGAAGGTGGCCAACgggggggg                                                                                                                                                                                          >  1:392194/1‑244 (MQ=255)
                                                                      cGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGCTTTCCCCTTTCTGCCGGGGGTAGGTGTACCGGTGGACTGAGCGAAATTATCCTTCGCAATACGGGAGATGGTCAACGTGGTGGTGAGAccgc                                                                                                                                                                                  >  1:26597/1‑248 (MQ=255)
                                                                                  cAGAGTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGGCAGCTGCCAGTTGTAGGTAATTGCGCCGTGGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGt                                                                                                                                                                      <  1:63420/251‑1 (MQ=255)
                                                                                     aTTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCg                                                                                                                                                                   <  2:233824/251‑1 (MQ=255)
                                                                                        aGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGAtt                                                                                                                                                                <  2:59125/251‑1 (MQ=255)
                                                                                             tAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAg                                                                                                                                                            <  1:212495/250‑1 (MQ=255)
                                                                                                 ggCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCGGTGGGCTGAGCGACATTTTCCTGCTCAATACGGGAGATGGCTCACGTGGGGGTCCGACGGCTTTTCATCCAGTCCCATTTCAGAcccc                                                                                                                                                       >  1:69339/1‑250 (MQ=255)
                                                                                                 ggCGTTTCCTTCGGTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAACTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCa                                                                                                                                                       <  2:16775/251‑1 (MQ=255)
                                                                                                   cGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGGACCGGTGGGCTGGGCGGCATTATCCTGCTCCATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTTATCCAGTCCGGtttaaggaccagc                                                                                                                                                     >  1:145248/1‑241 (MQ=255)
                                                                                                      ttCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTc                                                                                                                                                   <  2:347402/250‑1 (MQ=255)
                                                                                                            cgttagcttctGCAATATCGCGCGTTGCGCCACATTTCAGCTGTCAGTTGGCGGTAAGTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTGAGAGGCGCGTTCGCTGTTGCTGCTGTGTATAGGTGTACAGGTGGACTGAGCGACCTTATCCTGCTCAATACTGAAGATGGCTAACGTGGTGTTTTGACGGTTATTCATCCAGTCCGATTTCCGACCCAgtttttagtt                                                                                                                                             <  2:362223/242‑10 (MQ=255)
                                                                                                            cGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgtt                                                                                                                                            <  2:269359/250‑1 (MQ=255)
                                                                                                             gTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGGGGGGGTCAGACGGCTATTCATCCAGTCCGATTTTAGACCCAGCTCGTAgttgt                                                                                                                                          >  1:229598/1‑251 (MQ=255)
                                                                                                                            ataGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAg                                                                                                                           <  1:345170/251‑1 (MQ=255)
                                                                                                                                       cgcCAAATGTCACCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCCTCCCCCGCTTTATAGGCGGTTTCTCCCCTGCTGCCGGGGGTAGGGGCACCGGGGGGCTGTGCCACATTTTTCTGCTCCATACGGGGGGTGGGTAAAGGGGGGGTCAGGCGGCTATTTATCCTACCCGATTTCCGACCCAGCCCGTTGTTGGTACCGGGGTTTGGGACCAGataattt                                                                                                                >  2:136481/1‑247 (MQ=255)
                                                                                                                                          cAAATGTCAGCTGCCAGTTGTCGGGAATTGCGCCGTTGTGATCAAATTACACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCCCCGCTGCTCCCGGGGGAAGGGGGAACGGTGGACAGAGCGACAATATCTCGCTCCAAACCGAAGAAGGCTAAAGAGGAGGGCAGACGGGTATTCATCCACACCGATTTCCGAACCAGACCGGAGGTGGTACCGGGGATTGGAGCCAGAGAttttccc                                                                                                             >  2:53494/1‑247 (MQ=255)
                                                                                                                                            aaTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGGACCGGGGGACTGAGCGAGATTTTTCTTCTCAATACGGGAAAGTGCTAACGGGGTGGGGCGACGGCCATTCATCCCGTCCGGTTTTAGACCCAGCTCGGAATTGTTACCGGTGATTTGAGGCAGAAATTTGCCTg                                                                                                            >  1:419944/1‑250 (MQ=255)
                                                                                                                                               gTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACt                                                                                                        <  1:84113/251‑1 (MQ=255)
                                                                                                                                                   gCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGCTGCTGCCGGGGGTAGGTGGACCGGTGGACTGAGCGACATTATCCTGCACAATACGGGAGATGGGTAACGTGGGGGGCAGAAGGCTCTTCATCCAGTCCGATTTCAGACCCAGCCCGTAGTTGTTACCGGTGGTTGGCGGCAGATATTTTACTGAACTGTc                                                                                                     >  1:348115/1‑250 (MQ=255)
                                                                                                                                                      gCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGa                                                                                                 <  1:74878/251‑1 (MQ=255)
                                                                                                                                                        cAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATc                                                                                               <  2:56650/251‑1 (MQ=255)
                                                                                                                                                        cAGTTGTAGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGGCCTTGCTTCCGGGGGAAGGTGTAACGGGGGACGGAGCGAAATTATACTGCTCAATACGGGAGATGGCTAAAGTGGGGGTCAAACGGCTATTCATTCAGAACGATTTCAAACCCCGCTCGTAGGTGGTAACCGGGATTAGGGCCAGATTTTTGCCTGAACTGACACCATc                                                                                               >  1:19826/1‑251 (MQ=255)
                                                                                                                                                         aGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGGAGAAGGCTAACGTGGTGGGCAGACGGCTATTCATCCAGTCCGCTTTTAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATc                                                                                               >  1:218013/1‑250 (MQ=255)
                                                                                                                                                                ggTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTTCTGCCGGGGATAGGGGTACCGGTGGGCTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAAAGTGGGGGTCAGACGGCTATTCATCCAGTCCGATTTTAGACCCAGCTCGTAGTTGTTTCCGGGGATTTGGGGCAGATATTTTCCTGAACTGTCACGCTCAtttttcg                                                                                       >  1:330222/1‑248 (MQ=255)
                                                                                                                                                                   aaTTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCAGCCACCGCGTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                                                                                     <  2:368089/250‑1 (MQ=255)
                                                                                                                                                                          ccGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGATCTGTCACGATCATTTTGCGGCTGCAAAAt                                                                             <  1:315331/251‑1 (MQ=255)
                                                                                                                                                                              tGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATaga                                                                          <  1:334908/250‑1 (MQ=255)
                                                                                                                                                                                aGTTCAATTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGGCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGt                                                                       <  1:404844/251‑1 (MQ=255)
                                                                                                                                                                                aGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGTTGCCCGGGCTAGGGGTACCGGTGGACTGGGGGACATTTTTCTGCTCAATAAGGAGAATGGGAAACGGGGTGGTCAGACGGGCATTCATACCGTCCGTTTTTAGACCCAGCTCGTAGTTGTTACCGGGGATTGGAGCCAGATATTTGCCCGTAATGTCACGAGCATTTTGGGGCTGGGAAAttagggg                                                                       >  2:125584/1‑245 (MQ=255)
                                                                                                                                                                                  ttCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGGCGGGGGTAGGGGTAACGGGGGAATGAGCGGCATTTTTCTGGTCAATAAGGGAGAAGGCTAACGGGGGGGGCGGACGGGCATTTAACCCGTCCGATTTTAGACCCAGCACGTAGGTGGTAACGGGGGTTGGAGCCAGATATTTGCCTGAACTGGTACGAACAtttttggggtgggaaagagagggat                                                                     >  2:133368/1‑230 (MQ=255)
                                                                                                                                                                                    cAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGCTTCTGGCGGGGGGAGGTGGACCGGGGGACTGAGCGACCTTATCCCGCTCAATATGGAAGATGGGCAACGGGGGGGTCAGACGGCTATTTATCCCGTCCGACTTTAGGACCAACTCGTATTTGTTTCCGGGGATTGGAGCCCGATATTTTCCTTAACTGTCACCGTTATTTTTCGGCTGGGAAATAGAggggta                                                                    >  1:189456/1‑246 (MQ=255)
                                                                                                                                                                                          ccACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCg                                                             <  2:69339/251‑1 (MQ=255)
                                                                                                                                                                                                  ttACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGGCGGTTCTGCCGGGGGTAAGGGTAACGGGGGGATGAGCGGCATTATCCTTCTCAATAAGGAAGATGGCTAACGGGGGGGGCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTAACGGTGATTGGAGCCAGATATTTGCCTGAAATGTCAAGGACCTTTTTCGGCTGGGAAATAGAGGTATAGCTGGCGttgggtgg                                                     >  1:225890/1‑244 (MQ=255)
                                                                                                                                                                                                    actgactgTCCCATCCGCCGCTTTATAGGCGGTTTCGTCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGGCGGCTATTCAGCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGAc                                                    <  2:12856/250‑1 (MQ=255)
                                                                                                                                                                                                    actgactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTTCTTCCGGGGGTAGGTGTACCGGTGGGCTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGGGGGGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCCGCTCGTAGTTGTTACCGGTGATTGGGGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGGAAATAGAGGGATAGATGGCGTTGGGCGAcc                                                   >  1:8685/1‑251 (MQ=255)
                                                                                                                                                                                                    actgactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGAc                                                    <  1:218221/250‑1 (MQ=255)
                                                                                                                                                                                                       gactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAg                                                 <  1:258812/250‑1 (MQ=255)
                                                                                                                                                                                                                   cACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGt                                    <  1:33244/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGAGATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTa                 <  1:171723/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTa                 >  1:44477/1‑249 (MQ=255)
                                                                                                                                                                                                                                             ggtgctgcCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCAGTGATGTCAAACACCAGACCAGCGTAAGGCg            <  2:303167/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   gcCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtgtt        <  1:362515/249‑3 (MQ=255)
                                                                                                                                                                                                                                                    ccGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtggtgtg     >  1:444356/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     cGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACTCCAGACCAGCGTAAGGCgtggtgtgg    <  1:4239/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       gggATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACCCCAGACCAGCGTAAGGCGTGGTGTGGtt  >  2:335267/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
GCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTT  >  NC_000913/1159589‑1160084

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: