Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,336,9010TA35.6% 70.9 / 20.8 45Y44F (TAT→TTT) adiCarginine:agmatine antiporter
Reads supporting (aligned to +/- strand):  ref base T (25/4);  new base A (0/16);  total (25/20)
Fisher's exact test for biased strand distribution p-value = 7.49e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCT  >  NC_000913/4336651‑4337151
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          
gATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat                                                                                                                                                                                                                                                            >  1:144144/1‑251 (MQ=255)
gATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAAccatccat                                                                                                                                                                                                                                                            <  1:455444/251‑1 (MQ=255)
             aaTGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATccc                                                                                                                                                                                                                                               <  1:65983/251‑1 (MQ=255)
              aTGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCg                                                                                                                                                                                                                                              <  1:83120/251‑1 (MQ=255)
                                gaaaaaattAACTAAAATACCCTCCGCCAATGCCCCCCATGGGATTATTGCCGACCCAGCGGGCCAGCCAGTAGAGGACGTTGTTTTGATAACCGAGAAACGGGCAAAAGCAGCGGCGGGCGTAAGCGTAAGACCCCCCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCCCCGATAATCGTCACCACCCAACCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGttt                                                                                                                                                                                                                            <  1:468046/243‑1 (MQ=255)
                                  aagaagTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGc                                                                                                                                                                                                                          >  1:572875/1‑251 (MQ=255)
                                   agaagTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGc                                                                                                                                                                                                                          >  1:774893/1‑250 (MQ=255)
                                           cataaatctcccacgcaaagccccaccaTGGCGATATCGCAGATCAAAGGGCCCAGCCGGTAGAGGTCGTGGGTTTGATAACCGAGAAGCGGGCCAAGGCGGGGGCGGGCGTAGGCGTAAGAACCACCAGGATTTGGGTCAAGGAACGACTTTTTGGCGTAACCCATCGAGCGCCCCAGTGCACCAATAACCGTCCCCAACCATCAAAAAAGGGCACCCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAg                                                                                                                                                                                                                 <  2:29609/229‑1 (MQ=255)
                                               aaTATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCCCCAAACATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaa                                                                                                                                                                                                              >  1:741028/1‑250 (MQ=255)
                                               aaTATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGaaa                                                                                                                                                                                                              >  2:577193/1‑250 (MQ=255)
                                                  atCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGCGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGCACGACATTTTGGCGTATACCATCCAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATTGCAATCCCGCCAGTAGAGGCCAGGCTTGCAGGTAACAGAAAAac                                                                                                                                                                                                           >  1:750613/1‑250 (MQ=255)
                                                                                     gATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAggg                                                                                                                                                                        >  1:785506/1‑251 (MQ=255)
                                                                                              ngccagccaGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGGAACAAAAAAACACCTGACCCCAAAATATTTCCCGACACCCTCAAGGTTACgggg                                                                                                                                                                >  2:344473/2‑250 (MQ=255)
                                                                                                      agtaggggccgtggttttaaaaaCCGAGAAAGGGGCAAAAGCGGCGGGGGGGGTAAGCGAAAGACCCCCCAGGACTGGGGCCGAGGAACGACATTTTGGAGTATCCCAGCGAGAGCCCCATGCCACCGAAAAGCGTCCCCACCCACCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGccc                                                                                                                                                       <  2:137193/230‑1 (MQ=255)
                                                                                                                  tGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGCGACGGGGATTAAGCCCACTTTGTGagc                                                                                                                                            >  2:758694/1‑250 (MQ=255)
                                                                                                                              cGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGGGACGGGGATTAAGCCCACTTTGGGAGAATCAGCATCCGa                                                                                                                                >  2:42558/1‑250 (MQ=255)
                                                                                                                                       ccgatcggcggcggcgGGCGTAAGCGTAAAAACCACCAGGACTGGGGCCGAGGAACGAAATTTTGGCGTATCCCACCGAGACCCCCAGTGCACCGATAATCGCCACCAACCACCCAAAAAGGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAt                                                                                                                       <  1:760022/240‑1 (MQ=255)
                                                                                                                                        gcaaaggcggcggcgGGCGAAAGCGGAAAACCCCCCAGGACTTGGGTCGAGGAACGACATTTGGGCGTATACCATCGGGACCCCCGGTGCCCGGAAAATCGGCCCAACCCATCAAAAAAGGGCACTCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTa                                                                                                                     <  2:741028/242‑1 (MQ=255)
                                                                                                                                           aaaGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAAc                                                                                                                  >  1:242150/1‑251 (MQ=255)
                                                                                                                                            aaGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACt                                                                                                                 >  1:756475/1‑251 (MQ=255)
                                                                                                                                             agcagcGGCGGGCATAAGCGAAAGAACCACCAGGACGGGGGTCGAGGAACGACATTTTGGCGAATACCATCGAGACCCCAAGTGCCCCGAAAATCGCCACCAACCATCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTc                                                                                                                <  2:785506/251‑1 (MQ=255)
                                                                                                                                              gcagcGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAtagatcggg                                                                                                               >  2:760022/1‑243 (MQ=255)
                                                                                                                                                 ggggggggcgTAAGCGAAAGACCCACCAGGACTTGGGTCGAGGAACGACTTTTTGGCGTATCCCATCGAGAGCCCCAGTGCCCCGATAATCGTCACCAACCACCCAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGc                                                                                                            <  1:577193/246‑1 (MQ=255)
                                                                                                                                                  gggggggcgTAAGGGTAAGACCCACAAGAATTTGGGTCGAGGAACGACATTTTGGGGTATCCCATCGGGAGCCCCATACCCCCAATAATCTTCCCCACCCACCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCg                                                                                                           <  1:162125/247‑1 (MQ=255)
                                                                                                                                                  cggcggGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCg                                                                                                           >  1:394055/1‑251 (MQ=255)
                                                                                                                                                   ggcgggcgTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCg                                                                                                           >  2:460901/1‑250 (MQ=255)
                                                                                                                                                    gcgggcgAAAGCGAAAGACCCCCAAGGACTTGGGTCGAGGAACGACTTTTTGGCGTATACCACCGAGACCCCCAGTGCCCCGATAATCGTCACCACCCATCCATAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGaa                                                                                                         <  1:559045/251‑1 (MQ=255)
                                                                                                                                                      gggcgTAAGCGTAAGAACCCCAAGGTCTGGGGTCGAGGAACGACATTTGGGCTTAACCCAGCGAGAGCCCCAGTGCACCGATAATCGCCACCAACCAACCAAAAAGGACAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAgg                                                                                                       <  2:572875/251‑1 (MQ=255)
                                                                                                                                                            aagcataagACCCACCGGGACTGGGTGCGAGGAACGCCTTTTTGGCGTATACCATCGAAAGCCCCAGTGCACCGATAATCGCCACCACCCACCCAAAAATGGCAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCt                                                                                                 <  2:168558/251‑1 (MQ=255)
                                                                                                                                                                      accaccAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAAt                                                                                       >  2:485625/1‑251 (MQ=255)
                                                                                                                                                                             ggACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCt                                                                                                                                                                                                       >  2:320299/1‑132 (MQ=255)
                                                                                                                                                                              naCTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTAACAATAAATTACg                                                                               >  2:364104/2‑251 (MQ=255)
                                                                                                                                                                               aCTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGt                                                                              >  1:631819/1‑251 (MQ=255)
                                                                                                                                                                                   gggTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcat                                                                           >  1:671159/1‑250 (MQ=255)
                                                                                                                                                                                   gggTCGAAGAACGACTTTTTGGCGTTACCCACCGAGACCCCCAGGGCCCCGAAAATCCACCCCACCCAACCAAAAATGGAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcata                                                                          <  2:685187/251‑1 (MQ=255)
                                                                                                                                                                                         gAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGtcatatcat                                                                      >  1:610274/1‑250 (MQ=255)
                                                                                                                                                                                             naCGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGcc                                                                 >  2:363354/2‑251 (MQ=255)
                                                                                                                                                                                                    ttttGGCGTATACCTTCGAGACCCCCGGCCCACCGATAATCGTCACCACCCACCAAAAAATGGAAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCAtt                                                           <  1:42558/250‑1 (MQ=255)
                                                                                                                                                                                                     tttGGCATATACCATCGAGGGCCCAATGGCACCGATAATGGTCCCCAACCTTCCAAAAAGGGAAACCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTc                                                          <  1:502278/250‑1 (MQ=255)
                                                                                                                                                                                                                    ccGAGAGCCCCAGTGCACTAAAAACGGACACAAACCACCCAAAAAGGGCATCCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACGCCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAAGTAAACTGCTGCGAAGGCGAGCTTCGTGACAAGAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAattatt                                           <  1:79260/250‑1 (MQ=255)
                                                                                                                                                                                                                     cGAGAGCCCCAGTGCACCGATAATCGTCACAAACCATCAAAAAAGGACAAACCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAAttatta                                          <  1:297139/251‑1 (MQ=255)
                                                                                                                                                                                                                        gagCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCg                                       >  2:743938/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aTCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGc                    >  1:615425/1‑251 (MQ=255)
                                                                                                                                                                                                                                              cgcccccaaccaTCAAAAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCg                  <  1:634321/246‑1 (MQ=255)
                                                                                                                                                                                                                                              cGTCCCCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGc                    <  1:363354/249‑1 (MQ=255)
                                                                                                                                                                                                                                                  accaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat              >  2:785606/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  accaaccaTCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGATAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGttat              >  1:499861/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             catAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGAGACGGGGATTAACCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTCCGCCATATCATACGCCTGCATTCGATAGTGGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGct    >  2:272470/1‑251 (MQ=255)
                                                                                                                                                                                                                                                             catAAATGGCAATCCCGCCAGTAGAAACCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCACCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGct    >  1:634434/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               tAAATGGCAATCCCGCCAGTAGAGGCCAGGGTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGGGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGctct  >  1:87101/1‑251 (MQ=255)
                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                          
GATGGTTAATACCAATGGATCTTTCAGAATCGGGAAGAAGTAACTTAAATATCCTACGCCAATGACCACCATGGCGATATTGCCGATCCAGCAGGCCAGCCAGTAGAGGACGTTGGTTTGATAACCGAGAAACGGGCCAAAGCAGCGGCGGGCGTAAGCGTAAGAACCACCAGGACTTGGGTCGAGGAACGACATTTTGGCGTATACCATCGAGAGCCCCAGTGCACCGATAATCGTCACCAACCATCCATAAATGGCAATCCCGCCAGTAGAGGCCAGGTTTGCAGGTAACAGAAAAACACCTGACCCCATAATATTCCCCGACACCATCAGGGTGACGGGGATTAAGCCCACTTTGTGAGCATCAGCATCCGAAGACATAATTAAACTCCTGCGAAGGCGAGCTTCGTGACAATAAATTACGTCATATCATACGCCTGCATTCGATAGTTGAATTATTATCGACCGGGTTATCGCTGATGCCGTTATTAATATAGCTCT  >  NC_000913/4336651‑4337151

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