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

TAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTTT  >  NC_000913/1159595‑1160085
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                             
tAACTCCACCGCCGACCGACAACTCGGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCTGCTGCCAGGTGTCGGTAATGGCGCCGTTGAGTTCAAATTCCTCCCCTTTACTGACTGTCCCATCCTCCGCTTTATAGGCGGTTGCGCCGTTGCTGCCGGGGATAGgtgt                                                                                                                                                                                                                                                       <  2:275327/250‑1 (MQ=255)
      caccGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTGCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTGATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTg                                                                                                                                                                                                                                                <  2:450995/251‑1 (MQ=255)
       cccgccgACCGTCAACTCGGGCATGACAGGCATCCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGGCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTTGCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTgg                                                                                                                                                                                                                                               <  2:288745/250‑1 (MQ=255)
         cgccgACCGTCAACTCGTGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTCGCAGATTAGGAATAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAATTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCGCCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGAc                                                                                                                                                                                                                                             <  2:159914/251‑1 (MQ=255)
          gccgACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACt                                                                                                                                                                                                                                            <  1:407133/251‑1 (MQ=255)
                                aCAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAAt                                                                                                                                                                                                                     <  1:458316/251‑1 (MQ=255)
                                      aaCCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTCCTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGaa                                                                                                                                                                                                              <  1:42644/251‑1 (MQ=255)
                                       aCCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTTCTGGGGGTAAGTGTAACGGAGGAGTGAGCGAACTTATCCTGGTCAATACGGa                                                                                                                                                                                                               >  1:360205/1‑249 (MQ=255)
                                             aGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGc                                                                                                                                                                                                        <  2:4332/250‑1 (MQ=255)
                                                   ggAACATTTTAACCGTGGTGCGTGGCAGAGTAGGATTAACGGCGTTTCCGTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACgtg                                                                                                                                                                                                 <  2:260710/250‑1 (MQ=255)
                                                        aTTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTGCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTggtggt                                                                                                                                                                                            <  1:78599/251‑1 (MQ=255)
                                                                 gtggtgCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCt                                                                                                                                                                                   <  1:6215/251‑1 (MQ=255)
                                                                          ggCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATcc                                                                                                                                                                           <  1:249447/250‑1 (MQ=255)
                                                                               aTTAGGATTAACGGCGTTTCCTTCGTTATCCTCGGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCTTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCATTCCg                                                                                                                                                                     <  1:311971/251‑1 (MQ=255)
                                                                                                  ccTTCGTTAACCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGt                                                                                                                                                   <  1:265906/250‑1 (MQ=255)
                                                                                                            ccTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGGTGCCGGGGGGAGGGGGAAAGGGGGGCTGAGCGACATTTTACTGCTCAATAACGAAGATGGCGAACGTGGGGGGCAGAAAGAGATTCATCCAGTACGCGTTCAGACCCAGCGCGTATGTGTTAccc                                                                                                                                       >  1:122129/1‑250 (MQ=255)
                                                                                                                   aaTATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGATTTGCCGGTGGTGACGGGGGGAGGGGTACCGGGGGACTGAGCGGCATTTTCCTGCTCCATACGGAAGATGGGCAACGGGGGGGGCAGACGGCCATTAATCAAGTCCGATTTCCAAACCAACTCGTAGTTGTTTCCgggggatg                                                                                                                                >  2:469499/1‑245 (MQ=255)
                                                                                                                   aaTATAGCGCGGTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTg                                                                                                                                <  1:9481/251‑1 (MQ=255)
                                                                                                                             gTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTCTGCCCGGGGTAGGTGTACCGGGGGGCTGGGCGGCATTTTCCTGCTCCATAAGGAAGAAGGGTAACGTGGTGGGCGAACGGCCATTCATCCCGTCCGATTTCCGACCCCGCCCGGAGTTGTTACCGGGGGTTGGGGCCAgct                                                                                                                       >  1:125191/1‑248 (MQ=255)
                                                                                                                             gTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGGTACCGGGGGACTGAGCGACAATATCCTGCTCAATACGGAAGAAGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGGGATTGGAGCCCGat                                                                                                                       >  1:279114/1‑250 (MQ=255)
                                                                                                                                     aaaTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTg                                                                                                              <  2:474797/251‑1 (MQ=255)
                                                                                                                                       aTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTCCGCCGGGGGTAGGGGGAACCGGGGGCTGAGCGGCATTAACCTGCTCCATACAGGAGGAGGGTAAAGTGGGGGTCAGACGGCTATTTATCCCGGCCGAATTCAGACACAGCTCGGTGTTTTTACCgggggtttgaggccggtattttcctga                                                                                                             >  1:362527/1‑225 (MQ=255)
                                                                                                                                                 ccAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGAt                                                                                                  <  1:145172/251‑1 (MQ=255)
                                                                                                                                                   aGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTCTGCCGGGGGTAGGGGTAACGGGGGACTGGGCGACATTATCCTGCTCAATACGGGAGATGGCTAAAGTGGGGGGCAGACGGCTATTCATCCAGTCCGATTTCAGAACCAGATCGTAGTTTGTTCCGGGGGTTGGGGCCAGAAATTTGCCTGAACCGTCACCATct                                                                                                >  1:27844/1‑250 (MQ=255)
                                                                                                                                                       gTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCAttt                                                                                             <  1:438826/250‑1 (MQ=255)
                                                                                                                                                                  cgcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGGTCCAGCCCCCGCTGTATAGGCGGTTTCGCCGTTGCTGCCGGGGATATGTGTACCGGTGGACTGATCGACATTATCCTGCTCAATACGGAAGATGGCGAACGTGGTGGTCTGACGGCTATTCATCCAGTCCGATTTCAGAACCAGCTCTTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGGCACGATCATTTTGCGGCTGGaaa                                                                                <  2:317548/251‑1 (MQ=255)
                                                                                                                                                                   gcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCTCCGTTTCTGGCGGGGATAGGTGTACCCGGGGACTGAGCGGCATTTTCCTTCTCAATAAGGGAGATTGCTTAAGCGGTGGGCCGACCGCCATTTCTTCAGTCCGATTTCAGACCCAGCTCGTTGTTGGTTCCGGTGGTTGGGGGCCGGTATTTTCCTGTACTGTCCCGTTCATTTTTCGGCTgggaa                                                                                 >  2:123213/1‑247 (MQ=255)
                                                                                                                                                                               aaaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGGCGGGGGTAAGGGGACCGGGGGACTGAGCGACAATATCACGCTCAATACGGCAGATGGCGCAAGCGGGGGGAAGACGGCGATTACTCCGGTCCGATTTTACACCCAGCTCGCAGATTTTACCGGTGGTTGGAACCCGACATTTGTCTGAAACGTCCAGATCCTTTTTCGGCTGGAGAAAAGAggggtagc                                                                   >  1:20401/1‑245 (MQ=255)
                                                                                                                                                                               aaaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCCGGGGTAGGGGTACCGGGGGACTGAGCGACATTATCCTGCTCCATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGGCCGATTTCAGACCCCACCCGGAGGTGTTACCGGGGGTTGGGGGCCGAAATTTTTCTGAACTGTCACGATCATTTTTCGGCTGGGAAAATaagggatcgc                                                                    >  1:220342/1‑242 (MQ=255)
                                                                                                                                                                                aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGATGCCGGGGGTAGGGGTACCGGTGGACTGGGCGGCATTATACTGCTCCAAAAGGGAAATGGCTAACGGGGGGGTCAGACGGCAATTTCTCCAAACCGAATTCAGAACCCGCCCGGAGTTGTTAAAGGTGATTGGAGACAAATATTTTCCTGAACTGTCACGATCATTTTTCGGGTGGGAAAGAGATGtatatat                                                                  >  1:69129/1‑248 (MQ=255)
                                                                                                                                                                                aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCt                                                                  <  1:58420/251‑1 (MQ=255)
                                                                                                                                                                                     cACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGt                                                             <  1:123213/251‑1 (MQ=255)
                                                                                                                                                                                            ttACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTc                                                        <  1:425670/249‑1 (MQ=255)
                                                                                                                                                                                            ttACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTGCCGGGGGTAGGTGTACCGGTGGAATGAGCGGCATTATCCTGCTCAATACGGAAGATGGCTAACGTGGGGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGGGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGGAAATAGAGGGATAGCTGGCGGAGGTCgg                                                      >  1:273508/1‑250 (MQ=255)
                                                                                                                                                                                                 gactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGAAGGGGTACCCGGGGGCTGAGCGGCATTATCCCGCTCAATACGGAAGATGGCTAACGCGGTGGTCAGGCGGCTATTCCTCCCGTCCGATTTCAGACCCCGCTCGTTGTTGTTACCGGTGATTGGAGCCAGGTATTTGCCCGAACTGTCACCATCATTTTGCGGCTGGAAAAAAGAGGGATAGCTGGGGTAgggtgacccgg                                                 >  1:343994/1‑242 (MQ=255)
                                                                                                                                                                                                    tgTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGATAGGGGGACCGGTGGACCGAGCGACATTATCCTGCTCAATAAGGGAGATGGGTCACGGGGGGGGCCGACGGCTATTCATCCAGTCCGGTTTCAGACCCAGCTCGTAGGTGTTACCGGTGATTGGAGCCAGATATTTTCCTGAAATGTc                                                                                                       >  1:184156/1‑195 (MQ=255)
                                                                                                                                                                                                    tgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                                                                                       <  2:184156/195‑1 (MQ=255)
                                                                                                                                                                                                     gTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGt                                              <  2:290853/250‑1 (MQ=255)
                                                                                                                                                                                                       cccATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCa                                            <  1:469499/250‑1 (MQ=255)
                                                                                                                                                                                                         cATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTg                                         <  1:147687/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGTCCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTTGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGc              <  2:37416/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gcCGGGGATAGGTGTACCGGTGGACTGGGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGTCGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtggt         <  2:318935/249‑1 (MQ=255)
                                                                                                                                                                                                                                            gcGGGGATAGGTGTACCGGTGGGCTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGTTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtggtgtg      <  2:322877/250‑1 (MQ=255)
                                                                                                                                                                                                                                                ggATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGttt  <  1:127442/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                             
TAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTGTACCGGTGGACTGAGCGACATTATCCTGCTCAATACGGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTTGTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGGAAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTTT  >  NC_000913/1159595‑1160085

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