Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,159,8370AG25.7% 157.3 / 10.6 70I572T (ATC→ACC) fhuEferric coprogen/ferric rhodotorulic acid outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base A (11/41);  new base G (18/0);  total (29/41)
Fisher's exact test for biased strand distribution p-value = 1.49e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTG  >  NC_000913/1159597‑1160079
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           
acacCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGggggtaggggtaac                                                                                                                                                                                                                                                  >  1:478908/1‑249 (MQ=255)
   ccaccGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGt                                                                                                                                                                                                                                               <  2:231178/251‑1 (MQ=255)
      ccgccgACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTGTGCCGGGGGTAGGGGTACCgggggg                                                                                                                                                                                                                                            >  2:143724/1‑247 (MQ=255)
                   aCTCTGGCATGACAGGCAACCGATAGCTGGAGAACATTTTAAATGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACACCTTTACTGACTGTCCCATCAACCGCTTTATAGGCGGTTTTGGCGTTGGCGCCGGGGGTAGGGGTAACGGGGGACTGAgagaaatct                                                                                                                                                                                                                              >  1:154168/1‑243 (MQ=255)
                          cATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCt                                                                                                                                                                                                                       <  1:23354/251‑1 (MQ=255)
                           aTGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTTCTGCCGGGGGTAGGTGTAACGGGGGACTGAGGGACATTTTCCTGCTc                                                                                                                                                                                                                      >  1:474661/1‑251 (MQ=255)
                              aCAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGGCATTATCCTGCTCAAt                                                                                                                                                                                                                   >  1:137191/1‑251 (MQ=255)
                                       cGATAGCTGGTGACCATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTATCTGACTGTCCCATCCCCCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCTACATTATCCTGCTCAATACGGACGAt                                                                                                                                                                                                        <  2:368110/251‑1 (MQ=255)
                                        gATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGGTTCCGGGGGTAGGGGGACCGGGGGGCTGGGCGACATTTTTCTGCTCAATAAGGAAGATg                                                                                                                                                                                                       >  2:86467/1‑251 (MQ=255)
                                            gCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGGTTTGCCGCTGCTGCCCGGGAGAAGAGGACCGGGGGGATGAGGGACATTTTCCTGCTCAATACGGGAGATgggc                                                                                                                                                                                                    >  1:127685/1‑248 (MQ=255)
                                               ggTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAAc                                                                                                                                                                                                 <  1:456739/250‑1 (MQ=255)
                                                gTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGTGTACCGGGGGGCTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACgg                                                                                                                                                                                               >  1:13293/1‑250 (MQ=255)
                                                   aaCATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACgtggt                                                                                                                                                                                            <  2:530221/251‑1 (MQ=255)
                                                        tttAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTc                                                                                                                                                                                        <  2:194008/250‑1 (MQ=255)
                                                                 ggtgCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTAt                                                                                                                                                                             <  1:143724/251‑1 (MQ=255)
                                                                   tgCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCCGGGGGCTGAGCGACATTATCCTGCTCCATACGGGAGATGGCTAACGGGGGGGGCAGAAGGCTAttt                                                                                                                                                                           >  1:135328/1‑250 (MQ=255)
                                                                        ggCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCTCCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCa                                                                                                                                                                     <  2:14474/251‑1 (MQ=255)
                                                                            gATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGATAGGGGTAACGGTGGGCTGGGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGGGGGGGGCAGACGGCTAt                                                                                                                                                                             >  1:175601/1‑240 (MQ=255)
                                                                            gATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTAt                                                                                                                                                                             <  2:175601/240‑1 (MQ=255)
                                                                                 ggATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCCGTTCTGCCGGGGGAAGGTATACCGGTGGGCTGAGCGGCATTTTTCTGCTCAATAAGGAAGGTGGCTAAAGTGGTGGTCAGACGGCCATTCCTCCAGTCCGAttt                                                                                                                                                            >  1:366565/1‑251 (MQ=255)
                                                                                          gCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAg                                                                                                                                                  <  1:86467/251‑1 (MQ=255)
                                                                                               tCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCCTTGCTGCCGGGGGTAGGTGTTACCGTGGACTGGGCGAACTTATCCTGCCCAATAAGGGAGAAGGGTAACGTGGGGGGCAGAGGGGTTATCATCTAGTACGCTTTCAGACCCCGCTctg                                                                                                                                             >  2:145730/1‑249 (MQ=255)
                                                                                                  gaatttctgctctGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGATGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgt                                                                                                                                          <  2:466115/242‑1 (MQ=255)
                                                                                                    cggtATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCTCCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgtt                                                                                                                                         <  2:131021/247‑1 (MQ=255)
                                                                                                     gTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgttgt                                                                                                                                       <  1:330459/251‑1 (MQ=255)
                                                                                                     gTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGCTTCTCCGTTTCTGCCGGGGGTACGGGTTCCGGGGGACTGTGCGGCACTCTCCTTCTCAATACGGAAAATAGCTAACGTTGTGGGCAACGGGCTATTCATCCAGTcccgttttcgaccccgaactgttgtgt                                                                                                                                        >  2:96620/1‑226 (MQ=255)
                                                                                                       tATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAgttgtt                                                                                                                                      <  1:383885/250‑1 (MQ=255)
                                                                                                         tCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTGCCGGGGGTAGGGGTACCGGTGGGCTGAGCGACATTTTCCTGCTCCATACGGGAGATGGCTAACGTGGTGGGCAGACGGCTATTCATCCCGTCCGATTTCAGACCCCGCTCGCAGTTGTTAcc                                                                                                                                   >  1:332616/1‑251 (MQ=255)
                                                                                                               gCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGAt                                                                                                                             <  1:431905/251‑1 (MQ=255)
                                                                                                                        cgcgTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAg                                                                                                                    <  2:78646/251‑1 (MQ=255)
                                                                                                                        cgcgTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAg                                                                                                                    <  1:462317/251‑1 (MQ=255)
                                                                                                                         gcgTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCCCCCGTGGTGCCGGGGGTAGGTGTACCGGTGGGCTGAGGGGCATTATCCTGCCCAATAACGGAGAGGGGTAACGTGGGGGGCAGACGGCTATTTATCCCAGCCGAGTTCAGAACCCGCCCGTAGTTGTTTCCGGTGGTTGGAggccg                                                                                                                    >  2:71738/1‑246 (MQ=255)
                                                                                                                          cgTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGat                                                                                                                  <  2:135328/251‑1 (MQ=255)
                                                                                                                           gTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCCGGGAGAGGGGTACCGGTGGACTTGACGGCATTATCCTGCTCAATACGGAAGATGGGTAAAGTGGGGGGGAGACGGCTATTCATCCAGTCCGATTTCCGACCCAGCTCGTAGGTGTTACCCGTGATTGGAGCCAGat                                                                                                                  >  1:493730/1‑250 (MQ=255)
                                                                                                                            ttGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTCCCGTTGCTGGCGGGGGTAGGGGTACCCGGGGGCTGAGGGACATTATTCTGCGCAAAAAGGGAGATGGCTAAAGGGGGGGTCAGAAGGCCATTTATCCAGGCCGATTTCAGACCCAACTCGTAGGTGTTACCGGGGATTTGAggccggatt                                                                                                                >  1:260516/1‑243 (MQ=255)
                                                                                                                               cgcCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTg                                                                                                            <  1:277067/251‑1 (MQ=255)
                                                                                                                                 ccAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGcc                                                                                                          <  1:371821/251‑1 (MQ=255)
                                                                                                                                     aTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGATAGGTGTACCGGGGGACTGGGCGACCTTATCCTGCTCAATACGGAAGAAGGGTAAAGGGGTGGGCATACGGCTAATCATCCAAGCCGATTTTAGATCCCACTCGTAGGTGTTAACGGGGATTGGAGCCAGATATTTGCCTgt                                                                                                       >  1:588607/1‑249 (MQ=255)
                                                                                                                                     aTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCTCCGGTGCTGCCCGGGATAGGTGTACCGGTGGGCTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGTGGGGGGCAGACGGCTATTCATCCCGTCCGATTTTAGACCCAACTCGTAGTTGTTACCGGGGATTGGAGCCAGATATTTTCCTGa                                                                                                       >  1:163003/1‑250 (MQ=255)
                                                                                                                                       ggcAGCTGCCAGTTGTCTGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTAGAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACTTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAAc                                                                                                     <  2:235846/248‑1 (MQ=255)
                                                                                                                                       gTCAGCTGCCTGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAAc                                                                                                     <  2:127685/250‑1 (MQ=255)
                                                                                                                                       gTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACt                                                                                                    <  2:137191/251‑1 (MQ=255)
                                                                                                                                       gTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACt                                                                                                    <  2:474661/251‑1 (MQ=255)
                                                                                                                                                            aTTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTg                                                                               <  2:498589/251‑1 (MQ=255)
                                                                                                                                                            aTTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTGCCGGGGGAAGGTGTAACCGTGGGCAGAGCGACATTATTCTGGTCAATAAGCAAGATGGCTAACGTGGTGGTCAGAAGGCTATTCATCCAGTTCGAGTTAATAACCAGGACGGCGTTGTTAGCGGTGATTGGAGGCCGATAATTTGCCGGaatggcccgcatattttgcggctg                                                                                >  1:219031/1‑229 (MQ=255)
                                                                                                                                                                  ccGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAAt                                                                         <  1:230522/251‑1 (MQ=255)
                                                                                                                                                                   cGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATa                                                                        <  1:460894/251‑1 (MQ=255)
                                                                                                                                                                        aGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGt                                                                  <  2:163003/251‑1 (MQ=255)
                                                                                                                                                                              aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCt                                                            <  2:13293/251‑1 (MQ=255)
                                                                                                                                                                              aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCt                                                            <  1:61996/251‑1 (MQ=255)
                                                                                                                                                                              aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCCCCGTTGCTGCCCGGGGTAGGAGTACCCGCGGAGTGAGCGACATTATCCTGCTCAATACGGGAGATGGCCAAAGTGGGGGGCAGACGGGCATTCATCCAGCCCGGTTTTAGAACCAGCTTGTAGTTGTTTCCGGTGGTTTGAGCCAGAAATTTGTCTGAACTTTTCCGATCATTTTTGGGCTGGGAAAAAAAGGTATAGCt                                                            >  1:455113/1‑251 (MQ=255)
                                                                                                                                                                               aTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGCTTCGGCCGGGGAGAGGTGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGGGGGGGGCAGACGGCTATTCATCCCGTCCGATTTTAGACCCAGCTCGGAGGTGTTTCCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGGTGGAAAAAAGAGGAATAGCTg                                                           >  1:335058/1‑251 (MQ=255)
                                                                                                                                                                                    nccccTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCTTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCCATACGGAAGATGGCTAAAGTGGTGGTCAGACGGCTATTCAACCAGTCCGATTTCAGACCCAGCTCGCAGGTGTTACCGGTGAATGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGGAAATAGAGGGATAGCTGGCgg                                                       >  2:266955/2‑249 (MQ=255)
                                                                                                                                                                                       ccTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGt                                                   <  1:286475/251‑1 (MQ=255)
                                                                                                                                                                                          ttACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCg                                                 <  2:356873/250‑1 (MQ=255)
                                                                                                                                                                                              tgactgTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGTCGGGGGTATGTGTACCGGGGGGCTGAGGGGCATTATCCTGCTAAATAAGGGAGAAGGGCTCCGTGGGGGGCAGACGGCTATTCATCCCGTCCGATTTCAAAACCAGCCCGTAGTTTTTAACGGGGATTGGAGGCAGAATTTTGCCTGAAATGTCACGATCATTTTGCGGCTGGGGAATAGATGGTTAGCTgggcgaagtctgcacc                                            >  2:130761/1‑237 (MQ=255)
                                                                                                                                                                                              tgactgTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGGCGGGGGTAGGTGTCCCCGTGGGCTGAGTGGCATTTTCCTGCTCCATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTCGTTTTTTCCGGTGGTTGGAGCCCGGTATTTGCCCGAACTGTCACGATCATTTTTCGGGTGGGAAATAGAGGGATAGGTgggggagggcggcccc                                            >  2:219627/1‑237 (MQ=255)
                                                                                                                                                                                                 ctgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGGGTACCGGGGGACTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGGGGTTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGGAAATAGAGGTATAGCTGGCGTAGGTCGACCAGttt                                         >  1:221432/1‑250 (MQ=255)
                                                                                                                                                                                                          ccACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGGACCGGGGGGCTGAGCGGAATTATTCTGCTCAATACGGAAGATGGCTAACGTGGGGGTCAAGCGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGGGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTTCGGGTGGGAAATTGAGGTGTAACTGGCGTAGGGCGGCCAGTTGTCATTGAtt                                >  1:195511/1‑250 (MQ=255)
                                                                                                                                                                                                            aCCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGGCGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACGGTCACGATCATTTTGCGGCTGGAAAATAGAGGTGTAGCTGGCGTATGTCGACCAGTTGTCATTGATGTc                              <  2:157547/251‑1 (MQ=255)
                                                                                                                                                                                                            aCCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTc                              <  1:219627/251‑1 (MQ=255)
                                                                                                                                                                                                               tcTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCaa                            <  1:372423/249‑1 (MQ=255)
                                                                                                                                                                                                               gCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCaaa                           <  1:315602/251‑1 (MQ=255)
                                                                                                                                                                                                                  ttATAGGCGGTTTCGCCGTTGCTGCCGGGGGTAGGGGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGGAGATGGCTAAAGTGGTGGTCAGACGGCTATTCATCCAGGCCGATTTCAGACCCAGCTCGTGGTTGTTAACGGGGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGGAAATAGAGGGATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCaaaaac                        >  1:211256/1‑248 (MQ=255)
                                                                                                                                                                                                                   tataGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGATCTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCa                             <  1:560299/245‑1 (MQ=255)
                                                                                                                                                                                                                               cGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTaa           <  1:548999/251‑1 (MQ=255)
                                                                                                                                                                                                                                gCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGGGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACAACAGACCAGCGTa            >  2:423262/1‑249 (MQ=255)
                                                                                                                                                                                                                                  cGTTGCTGCAGGGGATAGGTGTACCGGTTGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAgtc        <  2:366565/251‑3 (MQ=255)
                                                                                                                                                                                                                                          gcCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtggtg  <  2:478908/250‑1 (MQ=255)
                                                                                                                                                                                                                                            cGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCg              <  1:381718/236‑1 (MQ=255)
                                                                                                                                                                                                                                            cGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCg              >  2:381718/1‑236 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                           
ACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTG  >  NC_000913/1159597‑1160079

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: