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.

TAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTTT  >  NC_000913/1159595‑1160085
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                             
tAACTCCACCGCCGACCGACAACTCGGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCTGCTGCCAGGTGTCGGTAATGGCGCCGTTGAGTTCAAATTCCTCCCCTTTACTGACTGTCCCATCCTCCGCTTTATAGGCGGTTGCGCCGTTGCTGCCGGGGATAGgtg‑t                                                                                                                                                                                                                                                       <  2:275327/250‑1 (MQ=255)
      caccGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTGCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTGATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTg                                                                                                                                                                                                                                                <  2:450995/251‑1 (MQ=255)
       cccgccgACCGTCAACTCGGGCATGACAGGCATCCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGGCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTTGCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTgg                                                                                                                                                                                                                                               <  2:288745/250‑1 (MQ=255)
         cgccgACCGTCAACTCGTGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTCGCAGATTAGGAATAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAATTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCGCCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGAc                                                                                                                                                                                                                                             <  2:159914/251‑1 (MQ=255)
          gccgACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACt                                                                                                                                                                                                                                            <  1:407133/251‑1 (MQ=255)
                                aCAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AAt                                                                                                                                                                                                                     <  1:458316/251‑1 (MQ=255)
                                      aaCCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTCCTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGaa                                                                                                                                                                                                              <  1:42644/251‑1 (MQ=255)
                                       aCCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTTCTGGGGGTAAGTG‑TAACGGAGGAGTGAGCGAACTTATCCTGGTC‑AATAC‑GGa                                                                                                                                                                                                               >  1:360205/1‑249 (MQ=255)
                                             aGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGc                                                                                                                                                                                                        <  2:4332/250‑1 (MQ=255)
                                                   ggAACATTTTAACCGTGGTGCGTGGCAGAGTAGGATTAACGGCGTTTCCGTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACgtg                                                                                                                                                                                                 <  2:260710/250‑1 (MQ=255)
                                                        aTTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTGCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTggtggt                                                                                                                                                                                            <  1:78599/251‑1 (MQ=255)
                                                                 gtggtgCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCt                                                                                                                                                                                   <  1:6215/251‑1 (MQ=255)
                                                                          ggCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATcc                                                                                                                                                                           <  1:249447/250‑1 (MQ=255)
                                                                               aTTAGGATTAACGGCGTTTCCTTCGTTATCCTCGGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCTTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCATTCCg                                                                                                                                                                     <  1:311971/251‑1 (MQ=255)
                                                                                                  ccTTCGTTAACCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGt                                                                                                                                                   <  1:265906/250‑1 (MQ=255)
                                                                                                            ccTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGGTGCCGGGGGGAGGGG‑GAAAGGGGGGCTGAGCGACATTTTACTGCTC‑AATAA‑CGAAGATGGCGAACGTGGGGGGCAGAAAGAGATTCATCCAGTACGCGTTCAGACCCAGCGCGTATGT‑GTTAccc                                                                                                                                       >  1:122129/1‑250 (MQ=255)
                                                                                                                   aaTATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGATTTGCCGGTGGTGACGGGGGGAGGGG‑TACCGGGGGACTGAGCGGCATTTTCCTGCTC‑CATAC‑GGAAGATGGGCAACGGGGGGGGCAGACGGCCATTAATCAAGTCCGATTTCCAAACCAACTCGTAGTT‑GTTTCCgggggatg                                                                                                                                >  2:469499/1‑245 (MQ=255)
                                                                                                                   aaTATAGCGCGGTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTg                                                                                                                                <  1:9481/251‑1 (MQ=255)
                                                                                                                             gTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTCTGCCCGGGGTAGGTG‑TACCGGGGGGCTGGGCGGCATTTTCCTGCTC‑CATAA‑GGAAGAAGGGTAACGTGGTGGGCGAACGGCCATTCATCCCGTCCGATTTCCGACCCCGCCCGGAGTT‑GTTACCGGGGGTTGGGGCCAgct                                                                                                                       >  1:125191/1‑248 (MQ=255)
                                                                                                                             gTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGGTAGGGG‑TACCGGGGGACTGAGCGACAATATCCTGCTC‑AATAC‑GGAAGAAGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGGGATTGGAGCCCGat                                                                                                                       >  1:279114/1‑250 (MQ=255)
                                                                                                                                     aaaTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTg                                                                                                              <  2:474797/251‑1 (MQ=255)
                                                                                                                                       aTGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTCCGCCGGGGGTAGGGGGAACCGG‑GGGCTGAGCGGCATTAACCTGCTC‑CATACAGGAGGA‑GGGTAAAGTGGGGGTCAGACGGCTATTTATCCCGGCCGAATTCAGACACAGCTCGGTGTT‑TTTACCgggggtttgaggccggtattttcctga                                                                                                             >  1:362527/1‑225 (MQ=255)
                                                                                                                                                 ccAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGAt                                                                                                  <  1:145172/251‑1 (MQ=255)
                                                                                                                                                   aGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTTCTGCCGGGGGTAGGGG‑TAACGGGGGACTGGGCGACATTATCCTGCTC‑AATAC‑GGGAGATGGCTAAAGTGGGGGGCAGACGGCTATTCATCCAGTCCGATTTCAGAACCAGATCGTAGTTTGTT‑CCGGGGGTTGGGGCCAGAAATTTGCCTGAACCGTCACCATct                                                                                                >  1:27844/1‑250 (MQ=255)
                                                                                                                                                       gTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCAttt                                                                                             <  1:438826/250‑1 (MQ=255)
                                                                                                                                                                  cgcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGGTCCAGCCCCCGCTGTATAGGCGGTTTCGCCGTTGCTGCCGGGGATATGTG‑TACCGGTGGACTGATCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCGAACGTGGTGGTCTGACGGCTATTCATCCAGTCCGATTTCAGAACCAGCTCTTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGGCACGATCATTTTGCGGCTGG‑aaa                                                                                <  2:317548/251‑1 (MQ=255)
                                                                                                                                                                   gcCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCTCCGTTTCTGGCGGGGATAGGTG‑TACCCGGGGACTGAGCGGCATTTTCCTTCTC‑AATAA‑GGGAGATTGCTTAAGCGGTGGGCCGACCGCCATTTCTTCAGTCCGATTTCAGACCCAGCTCGTTGTT‑GGTTCCGGTGGTTGGGGGCCGGTATTTTCCTGTACTGTCCCGTTCATTTTTCGGCTgggaa                                                                                 >  2:123213/1‑247 (MQ=255)
                                                                                                                                                                               aaaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGGCGGGGGTAAGGG‑GACCGGGGGACTGAGCGACAATATCACGCTC‑AATAC‑GGCAGATGGCGCAAGCGGGGGGAAGACGGCGATTACTCCGGTCCGATTTTACACCCAGCTCGCAGAT‑TTTACCGGTGGTTGGAACCCGACATTTGTCTGAAACGTCCAGATCCTTTTTCGGCTGG‑AGAAAAGAggggtagc                                                                   >  1:20401/1‑245 (MQ=255)
                                                                                                                                                                               aaaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCCGGGGTAGGGG‑TACCGGGGGACTGAGCGACATTATCCTGCTC‑CATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGGCCGATTTCAGACCCCACCCGGAGGT‑GTTACCGGGGGTTGGGGGCCGAAATTTTTCTGAACTGTCACGATCATTTTTCGGCTGGGAAAATaagggatcgc                                                                    >  1:220342/1‑242 (MQ=255)
                                                                                                                                                                                aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGATGCCGGGGGTAGGGG‑TACCGGTGGACTGGGCGGCATTATACTGCTCCAAAAG‑GGAA‑ATGGCTAACGGGGGGGTCAGACGGCAATTTCTCCAAACCGAATTCAGAACCCGCCCGGAGTT‑GTTAAAGGTGATTGGAGACAAATATTTTCCTGAACTGTCACGATCATTTTTCGGGTGG‑GAAAGAGATGtatatat                                                                  >  1:69129/1‑248 (MQ=255)
                                                                                                                                                                                aaTTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCt                                                                  <  1:58420/251‑1 (MQ=255)
                                                                                                                                                                                     cACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGt                                                             <  1:123213/251‑1 (MQ=255)
                                                                                                                                                                                            ttACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTc                                                        <  1:425670/249‑1 (MQ=255)
                                                                                                                                                                                            ttACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGGTGCTGCCGGGGGTAGGTG‑TACCGGTGGAATGAGCGGCATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGGGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGGGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑GAAATAGAGGGATAGCTGGCGGAGGTCgg                                                      >  1:273508/1‑250 (MQ=255)
                                                                                                                                                                                                 gactgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGGAAGGGG‑TACCCGGGGGCTGAGCGGCATTATCCCGCTC‑AATAC‑GGAAGATGGCTAACGCGGTGGTCAGGCGGCTATTCCTCCCGTCCGATTTCAGACCCCGCTCGTTGTT‑GTTACCGGTGATTGGAGCCAGGTATTTGCCCGAACTGTCACCATCATTTTGCGGCTGG‑AAAAAAGAGGGATAGCTGGGGTAgggtgacccgg                                                 >  1:343994/1‑242 (MQ=255)
                                                                                                                                                                                                    tgTCCCATCCACCGCTTTATAGGCGGTTTTGCCGTTGCTGCCGGGGATAGGGG‑GACCGGTGGACCGAGCGACATTATCCTGCTC‑AATAA‑GGGAGATGGGTCACGGGGGGGGCCGACGGCTATTCATCCAGTCCGGTTTCAGACCCAGCTCGTAGGT‑GTTACCGGTGATTGGAGCCAGATATTTTCCTGAAATGTc                                                                                                       >  1:184156/1‑195 (MQ=255)
                                                                                                                                                                                                    tgTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTc                                                                                                       <  2:184156/195‑1 (MQ=255)
                                                                                                                                                                                                     gTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGt                                              <  2:290853/250‑1 (MQ=255)
                                                                                                                                                                                                       cccATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCa                                            <  1:469499/250‑1 (MQ=255)
                                                                                                                                                                                                         cATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTg                                         <  1:147687/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGTCCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTTG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGc              <  2:37416/251‑1 (MQ=255)
                                                                                                                                                                                                                                           gcCGGGGATAGGTG‑TACCGGTGGACTGGGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGTCGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtggt         <  2:318935/249‑1 (MQ=255)
                                                                                                                                                                                                                                            gcGGGGATAGGTG‑TACCGGTGGGCTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGTTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCgtggtgtg      <  2:322877/250‑1 (MQ=255)
                                                                                                                                                                                                                                                ggATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGttt  <  1:127442/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                             
TAACACCACCGCCGACCGTCAACTCTGGCATGACAGGCAACCGATAGCTGGTGAACATTTTAACCGTGGTGCGTGGCAGATTAGGATTAACGGCGTTTCCTTCGTTATCCTCTGCAATATAGCGCGTTGCGCCAAATGTCAGCTGCCAGTTGTCGGTAATTGCGCCGTTGAGTTCAAATTCCACCCCTTTACTGACTGTCCCATCCACCGCTTTATAGGCGGTTTCGCCGTTGCTGCCGGGGATAGGTG‑TACCGGTGGACTGAGCGACATTATCCTGCTC‑AATAC‑GGAAGATGGCTAACGTGGTGGTCAGACGGCTATTCATCCAGTCCGATTTCAGACCCAGCTCGTAGTT‑GTTACCGGTGATTGGAGCCAGATATTTGCCTGAACTGTCACGATCATTTTGCGGCTGG‑AAAATAGAGGTATAGCTGGCGTAGGTCGACCAGTTGTCATTGATGTCAAACACCAGACCAGCGTAAGGCGTGGTGTGGTTT  >  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: ‑