Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,541,0370AG78.0% 82.1 / 21.0 41I518I (ATT→ATC) narZnitrate reductase Z subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (4/5);  new base G (19/13);  total (23/18)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.21e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAAT  >  NC_000913/1540790‑1541282
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
gtCGTTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGata                                                                                                                                                                                                                                                    <  2:418528/251‑1 (MQ=255)
   ggTTCCAGGCGGGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATATTGTGCCCCGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATAc                                                                                                                                                                                                                                                 <  1:293973/251‑1 (MQ=255)
                                    aGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACAGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGAGAATCATCGAGCGCCCATGCGCTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCCCGCCCGTCATTTGCTCAcccc                                                                                                                                                                                                                  >  2:417884/1‑249 (MQ=255)
                                             gtgGGCGCCGTTTTCCCTGGCCGGCATAGCTTGCCCAGCCGCCCCCCCTTTGCCCGACACAACCGCAGAAAATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTAATACCAGTGGTTAACCCCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTGTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCATGCTggagt                                                                                                                                                                                                        <  2:140598/250‑4 (MQ=255)
                                             gtgGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGgtg                                                                                                                                                                                                       >  1:378251/1‑251 (MQ=255)
                                                    cAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGtt                                                                                                                                                                                                >  1:305721/1‑251 (MQ=255)
                                                                     cATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCt                                                                                                                                                                               <  1:189158/251‑1 (MQ=255)
                                                                        aGTGTTCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATCCCAGTGGTTAACACCTGCGCCGAGGATAATCAGCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTGGATTTCAGCGTAATCTTTc                                                                                                                                                                            <  1:326686/251‑1 (MQ=255)
                                                                             gCCCAGCCTCGACCGCTTTCCCCGACACAACCGCAGAAGATCAGCAGATTGATCAGCCCGCGGTAGTTCATGTCCATGTGATACCAGTTGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGTAGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATATTGCGCACt                                                                                                                                                                       <  2:42075/251‑1 (MQ=255)
                                                                                    cgccACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCa                                                                                                                                                                <  1:269995/251‑1 (MQ=255)
                                                                                       cACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCt                                                                                                                                                             >  2:409916/1‑251 (MQ=255)
                                                                                         ccGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc                                                                                                                                                           >  2:203309/1‑251 (MQ=255)
                                                                                         ccGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTc                                                                                                                                                           <  1:409916/251‑1 (MQ=255)
                                                                                                    gACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGAt                                                                                                                                                >  1:150995/1‑251 (MQ=255)
                                                                                                        caACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCg                                                                                                                                              >  2:136829/1‑249 (MQ=255)
                                                                                                          aCCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGc                                                                                                                                          <  1:233756/251‑1 (MQ=255)
                                                                                                           ccGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGcc                                                                                                                                         >  2:143911/1‑251 (MQ=255)
                                                                                                                gaagaTCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAAt                                                                                                                                    >  1:172095/1‑251 (MQ=255)
                                                                                                                    aTCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGc                                                                                                                                >  1:234857/1‑251 (MQ=255)
                                                                                                                              tGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGAAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCaa                                                                                                                       >  2:223942/1‑250 (MQ=255)
                                                                                                                                 tcatcCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATc                                                                                                                    >  2:224036/1‑250 (MQ=255)
                                                                                                                                             aGTTCATGTCCATGTGATACCAGTGGTTACCACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGacc                                                                                                       >  1:384207/1‑251 (MQ=255)
                                                                                                                                             aGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGac                                                                                                        <  1:450694/250‑1 (MQ=255)
                                                                                                                                                tCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGaccac                                                                                                     <  1:481371/250‑1 (MQ=255)
                                                                                                                                                tCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGaccacc                                                                                                    >  1:140516/1‑251 (MQ=255)
                                                                                                                                                          gtgATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGgtgt                                                                                          >  2:473723/1‑251 (MQ=255)
                                                                                                                                                                            acCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGAGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGGGTTGCCATCGACTACCGTc                                                                         >  1:396762/1‑250 (MQ=255)
                                                                                                                                                                                              catcGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAAtt                                                      >  2:70679/1‑251 (MQ=255)
                                                                                                                                                                                                    gcgcCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGgcgc                                                 >  1:178815/1‑250 (MQ=255)
                                                                                                                                                                                                                     gggCGGTATCGGCAAATTCTCTGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTATCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTTCCACCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTgtgt                               <  1:273205/251‑1 (MQ=255)
                                                                                                                                                                                                                     gggCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTgtgt                               <  1:448797/251‑1 (MQ=255)
                                                                                                                                                                                                                      ggCGGTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGtt                              <  2:305721/251‑1 (MQ=255)
                                                                                                                                                                                                                          gTATCGGCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTAc                          <  1:136829/251‑1 (MQ=255)
                                                                                                                                                                                                                               ggCAAATTCACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGc                     >  1:155484/1‑251 (MQ=255)
                                                                                                                                                                                                                                       cACGGGCGATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGc             >  2:232037/1‑251 (MQ=255)
                                                                                                                                                                                                                                              gATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGac                                                                                                        <  1:488340/153‑1 (MQ=255)
                                                                                                                                                                                                                                              gATGGTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGac                                                                                                        >  2:488340/1‑153 (MQ=255)
                                                                                                                                                                                                                                               aTGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTc     <  1:392868/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gttcaaaTATACTGGCGCGGCACGTCGGTAATTTGCTCACCCCAGTCTGGGGTGTGCGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGATCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACTGGCCATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAAt  <  2:140516/246‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAAt  >  1:235372/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  gTTTCGATATACTGGCGCGGCACGCCAGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAAt  <  2:378251/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  gTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGAGTTTTACGCTGCGAAAATGCGGATTTTCAAt  >  2:33938/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
GTCGGTTCCAGTCGAGCGCAAAGGCCAGTGGCAACCAGCCGGTTTGTGGGCGCAGTTTTTCCTGGCCGACATAGTGTGCCCAGCCGCCACCGCTTTGCCCGACACAACCGCAGAAGATCAGCATATTGATCATCCCGCGGTAGTTCATGTCCATGTGATACCAGTGGTTAACACCTGCGCCGAGGATAATCATCGAGCGCCCATGCGTTTTATGGGCGGTATCGGCAAATTCACGGGCGATGGTTTCAATATACTGGCGCGGCACGCCGGTAATTTGCTCACCCCAGGCTGGGGTGTACGGTTTGATTTCAGCGTAATCTTTCGCACTGTTTTCATCTTCCAGCCCGCGATCGAGGCCGTAATTCGCCAGTACCAAATCATAAACGCTGACCACCGGACAGGTGTTGCCATCGACTAACGTCAGGTTTTTAACGGGCAATTGGCGCACCAGCACCGGATTGTGTTTTACGCTGCGAAAATGCGGATTTTCAAT  >  NC_000913/1540790‑1541282

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑