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.

GCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTT  >  NC_000913/1159589‑1160084
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
gCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATa                                                                                                                                                                                                                                                        <  1:247466/251‑1 (MQ=255)
    gTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGGTAGgtg                                                                                                                                                                                                                                                    >  1:56650/1‑251 (MQ=255)
         caccacCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCg                                                                                                                                                                                                                                               <  1:113817/251‑1 (MQ=255)
             accGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGGTCGGTGTACCGGTgg                                                                                                                                                                                                                                           >  1:384460/1‑251 (MQ=255)
              ccgccgACCGTCAACTCGGTCAGGACTGGCAACTGGTAGCTCGTGGAAGTTTTAGTCGTGGTGCGTGGCAGATTTGGGTTAACGGTGTTTCCTGCGTTATCCTCGGCAATATAGCGCGTTGCTCCAAATGTCAGCTGCCAGTTTTCGGTTAGTGCTCTTTTTAGTTGAAATTCCACCCCTTTACAGACTGTCGCATCCTCCTCGTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGa                                                                                                                                                                                                                                          <  2:240954/251‑1 (MQ=255)
                gccgACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACt                                                                                                                                                                                                                                        <  2:26305/251‑1 (MQ=255)
                           aCTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTTCTGTCGGGGATAgggggcccgggggaccgcgcgccattc                                                                                                                                                                                                                             >  1:233824/1‑226 (MQ=255)
                                ggCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATc                                                                                                                                                                                                                           <  1:363457/248‑1 (MQ=255)
                                ggCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTg                                                                                                                                                                                                                        <  2:86485/251‑1 (MQ=255)
                                   aTGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCCCGTTG‑TGCCGGGGATAGGGGTACCGGTGGACGGGGTGGCATTATCCTggcc                                                                                                                                                                                                                     >  2:362515/1‑248 (MQ=255)
                                         ggCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTTTCC‑GTTGCTGCCGGGGGTAGGGGTACCGGTGGGCTGGGGGGCATTATCCTTCTCCATACg                                                                                                                                                                                                               >  2:120843/1‑251 (MQ=255)
                                             aCCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTTACCTGCTCAATACggg                                                                                                                                                                                                             >  1:70329/1‑248 (MQ=255)
                                              ccGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGa                                                                                                                                                                                                          <  1:120843/251‑1 (MQ=255)
                                                 aTAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATgg                                                                                                                                                                                                       <  1:235221/251‑1 (MQ=255)
                                                              aTTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGGTAGGGGTACCGGGGGGCTGGGCGACATTATCCCGCTCAAAACGGAAGGTGGCCAACgggggggg                                                                                                                                                                                          >  1:392194/1‑244 (MQ=255)
                                                                      cGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGCTTTC‑CC‑CTTTCTGCCGGGGGTAGGTGTACCGGTGGACTGAGCGAAATTATCCTTCGCAATACGGGAGATGGTCAACGTGGTGGTGAGAccgc                                                                                                                                                                                  >  1:26597/1‑248 (MQ=255)
                                                                                  cAGAGTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGGCAGCTGCCAGTTGTAGGTAATTGCGCCGTGGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGt                                                                                                                                                                      <  1:63420/251‑1 (MQ=255)
                                                                                     aTTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCg                                                                                                                                                                   <  2:233824/251‑1 (MQ=255)
                                                                                        aGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGAtt                                                                                                                                                                <  2:59125/251‑1 (MQ=255)
                                                                                             tAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAg                                                                                                                                                            <  1:212495/250‑1 (MQ=255)
                                                                                                 ggCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGGTAGGTGTACCGGTGGGCTGAGCGACATTTTCCTGCTCAATACGGGAGATGGCTCACGTGGGGGTCCGACGGCTTTTCATCCAGTCCCATTTCAGAcccc                                                                                                                                                       >  1:69339/1‑250 (MQ=255)
                                                                                                 ggCGTTTCCTTCGGTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAACTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCa                                                                                                                                                       <  2:16775/251‑1 (MQ=255)
                                                                                                   cGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGGTAGGGGGACCGGTGGGCTGGGCGGCATTATCCTGCTCCATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTTATCCAGTCCGGtttaaggaccagc                                                                                                                                                     >  1:145248/1‑241 (MQ=255)
                                                                                                      ttCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTc                                                                                                                                                   <  2:347402/250‑1 (MQ=255)
                                                                                                            cgttagcttctGCAATATCGCGCGTTGCGCCACATTTCAGCTGTCAGTTGGCGGTAAGTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTGAGAGGCGC‑GTTCGCT‑GTTGCTGCTGTGTATAGGTGTACAGGTGGACTGAGCGACCTTATCCTGCTCAATACTGAAGATGGCTAACGTGGTGTTTTGACGGTTATTCATCCAGTCCGATTTCCGACCCAgtttttagtt                                                                                                                                             <  2:362223/242‑10 (MQ=255)
                                                                                                            cGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAg‑tt                                                                                                                                            <  2:269359/250‑1 (MQ=255)
                                                                                                             gTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTTGCC‑GTTGCTGCCGGGGGTAGGGGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGGGGGGGTCAGACGGCTATTCATCCAGTCCGATTTTAGACCCAGCTCGTAg‑ttgt                                                                                                                                          >  1:229598/1‑251 (MQ=255)
                                                                                                                            ataGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAg                                                                                                                           <  1:345170/251‑1 (MQ=255)
                                                                                                                                       cgcCAAATGTCACCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCCTCCCCCGCTTTATAGGCGG‑TTTCTCC‑CCTGCTGCCGGGGGTAGGGGCACCGGGGGGCTGTGCCACATTTTTCTGCTCCATACGGGGGGTGGGTAAAGGGGGGGTCAGGCGGCTATTTATCCTACCCGATTTCCGACCCAGCCCGTTG‑TTGGTACCGGGGTTTGGGACCAGataattt                                                                                                                >  2:136481/1‑247 (MQ=255)
                                                                                                                                          cAAATGTCAGCTGCCAGTTGTCGGGAATTGCGCCGTTGTGATCAAATTACACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCCCC‑GCTGCTCCCGGGGGAAGGGGGAACGGTGGACAGAGCGACAATATCTCGCTCCAAACCGAAGAAGGCTAAAGAGGAGGGCAGACGGGTATTCATCCACACCGATTTCCGAACCAGACCGGAG‑GTGGTACCGGGGATTGGAGCCAGAGAttttccc                                                                                                             >  2:53494/1‑247 (MQ=255)
                                                                                                                                            aaTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGGTAGGTGGACCGGGGGACTGAGCGAGATTTTTCTTCTCAATACGGGAAAGTGCTAACGGGGTGGGGCGACGGCCATTCATCCCGTCCGGTTTTAGACCCAGCTCGGAA‑TTGTTACCGGTGATTTGAGGCAGAAATTTGCCTg                                                                                                            >  1:419944/1‑250 (MQ=255)
                                                                                                                                               gTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACt                                                                                                        <  1:84113/251‑1 (MQ=255)
                                                                                                                                                   gCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GCTGCTGCCGGGGGTAGGTGGACCGGTGGACTGAGCGACATTATCCTGCACAATACGGGAGATGGGTAACGTGGGGGGCAGAAGGCTCTTCATCCAGTCCGATTTCAGACCCAGCCCGTAG‑TTGTTACCGGTGGTTGGCGGCAGATATTTTACTGAACTGTc                                                                                                     >  1:348115/1‑250 (MQ=255)
                                                                                                                                                      gCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGa                                                                                                 <  1:74878/251‑1 (MQ=255)
                                                                                                                                                        cAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATc                                                                                               <  2:56650/251‑1 (MQ=255)
                                                                                                                                                        cAGTTGTAGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTTGGC‑CTTGCTTCCGGGGGAAGGTGTAACGGGGGACGGAGCGAAATTATACTGCTCAATACGGGAGATGGCTAAAGTGGGGGTCAAACGGCTATTCATTCAGAACGATTTCAAACCCCGCTCGTAGGTGGTAACCGG‑GATTAGGGCCAGATTTTTGCCTGAACTGACACCATc                                                                                               >  1:19826/1‑251 (MQ=255)
                                                                                                                                                         aGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGGAGAAGGCTAACGTGGTGGGCAGACGGCTATTCATCCAGTCCGCTTTTAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATc                                                                                               >  1:218013/1‑250 (MQ=255)
                                                                                                                                                                ggTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTTCTGCCGGGGATAGGGGTACCGGTGGGCTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAAAGTGGGGGTCAGACGGCTATTCATCCAGTCCGATTTTAGACCCAGCTCGTAG‑TTGTTTCCGGGGATTTGGGGCAGATATTTTCCTGAACTGTCACGCTCAtttttcg                                                                                       >  1:330222/1‑248 (MQ=255)
                                                                                                                                                                   aaTTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCAGCCACCGCGTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGc                                                                                     <  2:368089/250‑1 (MQ=255)
                                                                                                                                                                          ccGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGATCTGTCACGATCATTTTGCGGCTGCAAAAt                                                                             <  1:315331/251‑1 (MQ=255)
                                                                                                                                                                              tGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATaga                                                                          <  1:334908/250‑1 (MQ=255)
                                                                                                                                                                                aGTTCAATTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGGCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGt                                                                       <  1:404844/251‑1 (MQ=255)
                                                                                                                                                                                aGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGTTGCCCGGGCTAGGGGTACCGGTGGACTGGGGGACATTTTTCTGCTCAATAAGGAGAATGGGAAACGGGGTGGTCAGACGGGCATTCATACCGTCCGTTTTTAGACCCAGCTCGTAG‑TTGTTACCGGGGATTGGAGCCAGATATTTGCCCGTAATGTCACGAGCATTTTGGGGCTGGGAAAttagggg                                                                       >  2:125584/1‑245 (MQ=255)
                                                                                                                                                                                  ttCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGGCGGGGGTAGGGGTAACGGGGGAATGAGCGGCATTTTTCTGGTCAATAAGGGAGAAGGCTAACGGGGGGGGCGGACGGGCATTTAACCCGTCCGATTTTAGACCCAGCACGTAG‑GTGGTAACGGGGGTTGGAGCCAGATATTTGCCTGAACTGGTACGAACAtttttggggtgggaaagagagggat                                                                     >  2:133368/1‑230 (MQ=255)
                                                                                                                                                                                    cAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GCTTCTGGCGGGGGGAGGTGGACCGGGGGACTGAGCGACCTTATCCCGCTCAATATGGAAGATGGGCAACGGGGGGGTCAGACGGCTATTTATCCCGTCCGACTTTAGGACCAACTCGTAT‑TTGTTTCCGGGGATTGGAGCCCGATATTTTCCTTAACTGTCACCGTTATTTTTCGGCTGGGAAATAGAggggta                                                                    >  1:189456/1‑246 (MQ=255)
                                                                                                                                                                                          ccACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCg                                                             <  2:69339/251‑1 (MQ=255)
                                                                                                                                                                                                  ttACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGGC‑GGTTCTGCCGGGGGTAAGGGTAACGGGGGGATGAGCGGCATTATCCTTCTCAATAAGGAAGATGGCTAACGGGGGGGGCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTAACGGTGATTGGAGCCAGATATTTGCCTGAAATGTCAAGGACCTTTTTCGGCTGGGAAATAGAGGTATAGCTGGCGttgggtgg                                                     >  1:225890/1‑244 (MQ=255)
                                                                                                                                                                                                    actgactgTCCCATCCGCCGCTTTATAGGCGG‑TTTCGTC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGGCGGCTATTCAGCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGAc                                                    <  2:12856/250‑1 (MQ=255)
                                                                                                                                                                                                    actgactgTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTTCTTCCGGGGGTAGGTGTACCGGTGGGCTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGGGGGGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCCGCTCGTAG‑TTGTTACCGGTGATTGGGGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGGAAATAGAGGGATAGATGGCGTTGGGCGAcc                                                   >  1:8685/1‑251 (MQ=255)
                                                                                                                                                                                                    actgactgTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGAc                                                    <  1:218221/250‑1 (MQ=255)
                                                                                                                                                                                                       gactgTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAg                                                 <  1:258812/250‑1 (MQ=255)
                                                                                                                                                                                                                   cACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGt                                    <  1:33244/251‑1 (MQ=255)
                                                                                                                                                                                                                                        cGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGAGATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTa                 <  1:171723/250‑1 (MQ=255)
                                                                                                                                                                                                                                         gCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTa                 >  1:44477/1‑249 (MQ=255)
                                                                                                                                                                                                                                             ggtgctgcCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCAGTGATGTCAAACACCAGACCAGCGTAAGGCg            <  2:303167/249‑1 (MQ=255)
                                                                                                                                                                                                                                                   gcCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtgtt        <  1:362515/249‑3 (MQ=255)
                                                                                                                                                                                                                                                    ccGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtggtgtg     >  1:444356/1‑251 (MQ=255)
                                                                                                                                                                                                                                                     cGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACTCCAGACCAGCGTAAGGCgtggtgtgg    <  1:4239/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       gggATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACCCCAGACCAGCGTAAGGCGTGGTGTGGtt  >  2:335267/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                        
GCCAGTTAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGG‑TTTCGCC‑GTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAG‑TTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTT  >  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: ‑