Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0129713,561,6960AG28.2% 118.3 / 36.9 71G21G (GGT→GGC) ECD_RS25800type I toxin‑antitoxin system toxin Ldr family protein
Reads supporting (aligned to +/- strand):  ref base A (35/16);  new base G (8/12);  total (43/28)
Fisher's exact test for biased strand distribution p-value = 3.34e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGTCACGTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATTAAGCCCCCGCAACATTGCGCCTCACC‑GGAATCCTTCCGGCTTGACGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGTCTC  >  NC_012971/3561401‑3561986
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                   
ccAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGTCACGTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCa                                                                                                                                                                                                                                                                                                  >  1:191894/1‑299 (MQ=255)
                             aCATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGTCACGTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCTCTTGCCCCTATTTCCTTCTTCTTACATAACCGGCCGCCGCTTTCCACCCCTGCTGCCTGGTTTAGGCTCAAATTGCCATCACGGCAGCCTGACGCACGCATGGAACGTTAATTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGAAATGAACGGAGCCGCTaaa                                                                                                                                                                                                                                                                             >  2:460526/1‑291 (MQ=14)
                             aCATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGTCACGTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCa                                                                                                                                                                                                                                                                     >  1:476978/1‑299 (MQ=255)
                                                ccTTTCCGGCTTGCCGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGTCACGTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCAAGAAGGCCATGCCCAGTTcc                                                                                                                                                                                                                                                  >  2:409835/1‑298 (MQ=255)
                                                                   aTTCTACGACTTTTTTTCTGTTCCTGAATCCGTCACTTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCTCTTGCAACTGTTGCATTTTTCTTACTGAAGAGGCCGTGGAGGTTTACCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAgc                                                                                                                                                                                                                                  <  1:24932/296‑1 (MQ=11)
                                                                      cTACGACTGTTTTTCTGTTCCGGAATCCGTCACGTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCACGGCACATTACTTCCGCTTGTTCAGACAGTTCACGATCATACTGGCAAGAATACCAGCAATGAAAGGAGCCGCTAAATCATGCAAGAAGGCCATGCACAGTTCTGCGAGCGTCATATAgccgc                                                                                                                                                                                                                               >  1:479305/1‑296 (MQ=255)
                                                                                                      cGTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGNCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCtt                                                                                                                                                                                            >  2:550406/1‑299 (MQ=39)
                                                                                                        tCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAgc                                                                                                                                                                                                                                  >  1:4976/1‑259 (MQ=37)
                                                                                                        tCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTATTCAGCCAGTTCACGATCATACTGGCAAGAATACCAACAATAAACGAAGCCGCTAACTCATGCCGGACGGCAATGCCCAGTTCTGCGAGCGTCATATAGCCGAATGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTa                                                                                                                                                                                         >  1:675951/1‑300 (MQ=11)
                                                                                                             aCTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGtt                                                                                                                                                                                    >  1:523175/1‑300 (MQ=255)
                                                                                                              ctctGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGttt                                                                                                                                                                                   >  1:606561/1‑300 (MQ=255)
                                                                                                                  gtgCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCCGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCAGCTAAATCATGCAAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGAt                                                                                                                                                                               >  1:101681/1‑300 (MQ=255)
                                                                                                                         gcATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAACGTctcctagccacatgta                                                                                                                                                                                                                          >  2:708642/1‑235 (MQ=255)
                                                                                                                              ccGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGNAGCCTGACGCCCGCATGGCACGTTACTTCCGCNTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAgc                                                                                                                                                                                                                                  <  2:570141/237‑1 (MQ=37)
                                                                                                                              ccGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAgc                                                                                                                                                                                                                                  >  1:570141/1‑237 (MQ=37)
                                                                                                                                   gAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAg                                                                                                                                                                                                                                   <  1:294873/231‑1 (MQ=37)
                                                                                                                                   gAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAg                                                                                                                                                                                                                                   <  1:294901/231‑1 (MQ=37)
                                                                                                                                   gAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAg                                                                                                                                                                                                                                   >  2:294873/1‑231 (MQ=37)
                                                                                                                                   gAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAg                                                                                                                                                                                                                                   >  2:294901/1‑231 (MQ=37)
                                                                                                                                     atatTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCAAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCgt                                                                                                                                                              >  1:732944/1‑298 (MQ=255)
                                                                                                                                          atCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGtgttg                                                                                                                                                       <  2:476978/300‑1 (MQ=255)
                                                                                                                                                         ttGCATTTTTCTTACGAAAAATGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGTGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGGTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATAcc                                                                                                                                          <  1:239565/298‑1 (MQ=255)
                                                                                                                                                            cgcagtaagaaCGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGCTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAgg                                                                                                                                                  <  1:705718/277‑1 (MQ=11)
                                                                                                                                                            cgcagtaagaaCGAAAAAGGCCGCAGAGTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGACACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGCTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAgg                                                                                                                                                  >  2:705718/11‑287 (MQ=11)
                                                                                                                                                             aTTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGNCGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCaa                                                                                                                                     <  1:165905/299‑1 (MQ=255)
                                                                                                                                                              tttttCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACNCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACg                                                                                                                                   <  1:550406/300‑1 (MQ=255)
                                                                                                                                                              tttttCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACg                                                                                                                                   <  2:523175/300‑1 (MQ=255)
                                                                                                                                                              tttttCCTACGAAAAAGGCCGCAGAGTTTCACCCCTTCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCACGGCACGTTACCTCCGCTTGTTCAGCTAGTTCGCGATCATACTGTTAAGAATGCCAGCAGTGGCCGGAGCCGCTAGATCATGCCATAATGCTAGGTCTAGGTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACg                                                                                                                                   <  2:373246/300‑1 (MQ=14)
                                                                                                                                                                     tgtggagaagatacaggagTTTCCCCCCTGCGGCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTc                                                                                                                                                     <  2:688456/259‑1 (MQ=255)
                                                                                                                                                                                       gTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGCTCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCAAGCaa                                                                                                          >  1:69403/1‑300 (MQ=255)
                                                                                                                                                                                       gTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAGATTATGCACGCCTACAATGCCCCGTTCTTCGAGTGTGGTGTGGGTGCATGTATTGTAATGACAACGTTTAGTGGATATTCTTGTATTGTTTGGATGCAAGAATGGCCCCAGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCaa                                                                                                          >  1:283438/1‑300 (MQ=14)
                                                                                                                                                                                             ccccTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTTAGCCATTTCCGGATCATCCTGGCCGGCATACCTGCAATGACCGGAGCCGCTCAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGTGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAAt                                                                                                     >  1:27092/1‑299 (MQ=17)
                                                                                                                                                                                             ccccTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAAt                                                                                                     >  1:27188/1‑299 (MQ=255)
                                                                                                                                                                                             ccccTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTCGTCCAGCCAGTTCACGATCCTAATGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGTAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAAt                                                                                                     >  1:27039/1‑299 (MQ=255)
                                                                                                                                                                                                    ggCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCC‑CATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGccacca                                                                                              >  1:350890/1‑298 (MQ=255)
                                                                                                                                                                                                    ggCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCNGCAATGACCGGAGCCGCTAAATCATGCCAGAANGCCATGCCCAGTTCTGCGAACGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATAGCCCCCGTGGTGTTGTCAGGTGCATACCTGCaa                                                                                                                                     >  2:558363/1‑260 (MQ=255)
                                                                                                                                                                                                    ggCCTGGTCCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCaa                                                                                                                                     <  1:558363/260‑1 (MQ=255)
                                                                                                                                                                                                           tccGGGCGCAAATTGCCGTTACGGCAGCCTGATACCCGCATGACACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCAGACTGGCAAGAAGACCAGCAAGGACCGGAGCCGCTAAACCGTGCCA‑AAGGCCATGCCCAGTTCTGCGAGCGTCATATAgc                                                                                                                                                                                                                                  <  2:285214/157‑1 (MQ=11)
                                                                                                                                                                                                                   cAAATTGCCATTACGGCAGCCTGACGCCTGCCTGGCATGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATATTGGCAAGCATACCAGCCATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTACAGCAACCAATGCCACCAAGGATTAAGCCCCCGc                                                                              >  2:209746/1‑300 (MQ=17)
                                                                                                                                                                                                                   cAAATTGCCATTACGGCAGCATGACGCCCGCATGGCACGTTACTTACGATTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCACCGACCGGCGCCGCTAAATTATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGACTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTACAGCAACCAATGCCACCAGGGATTAAGCACCCGc                                                                              >  2:209980/1‑300 (MQ=17)
                                                                                                                                                                                                                          ccATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATTAAGCCCCCGCAACAtt                                                                        >  1:12928/1‑299 (MQ=255)
                                                                                                                                                                                                                          ccATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCAACAGGGATTAAGCCCCCGCAACAtt                                                                        >  1:251460/1‑299 (MQ=255)
                                                                                                                                                                                                                           cATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGa                                                                                                                                                                                                                                              <  1:70281/132‑1 (MQ=18)
                                                                                                                                                                                                                           cATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGa                                                                                                                                                                                                                                              <  1:66770/132‑1 (MQ=18)
                                                                                                                                                                                                                           cATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGa                                                                                                                                                                                                                                              >  2:70281/1‑132 (MQ=18)
                                                                                                                                                                                                                           cATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGa                                                                                                                                                                                                                                              >  2:66770/1‑132 (MQ=40)
                                                                                                                                                                                                                                      ccTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAgc                                                                                                                                                                                                                                  <  2:397368/133‑1 (MQ=18)
                                                                                                                                                                                                                                      ccTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAgc                                                                                                                                                                                                                                  <  2:397417/133‑1 (MQ=18)
                                                                                                                                                                                                                                      ccTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAgc                                                                                                                                                                                                                                  >  1:397368/1‑133 (MQ=40)
                                                                                                                                                                                                                                      ccTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGACCGCCTCGCCATGCCGACGCTTCACAGCTATTACTGTATTGCTTGGATGCAAGAATGGCCAACGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTATATACAGCAACCAATGCAAACAGGGATTAAGCCCACGCAACATTGCGCCTCACC‑gaa                                                          >  1:397417/1‑298 (MQ=14)
                                                                                                                                                                                                                                          aCGCTCGCACGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGc                                                                                                                                               <  1:717521/212‑1 (MQ=255)
                                                                                                                                                                                                                                                 cATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCCGTTCCGCGAGCGTCATATAGCAGAATGGATTGTAATGAAAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTATCCAGCAACCAATGCCAACAGGGATTAAGCCCCCGCAACATTGCGCCTCAAC‑GGAATCCTTCCGGc                                               >  1:252827/1‑300 (MQ=17)
                                                                                                                                                                                                                                                     acACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGCTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGCATTGTCTGGATTCAAGACTAGCCCCCGTGGTGTTGTCAGGTGc                                                                                                                                               >  2:717521/2‑201 (MQ=11)
                                                                                                                                                                                                                                                        cGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGCTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACg                                                                                                                                   >  1:33381/1‑210 (MQ=255)
                                                                                                                                                                                                                                                        cGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTATCTCCAGCAACCAATGCAACCAGGGATTAAGCCCCCGCAACATTGCGCCTCACC‑GGAATCCTTCCGGCTTGaagc                                        >  1:584640/1‑297 (MQ=255)
                                                                                                                                                                                                                                                        cGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGCTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACg                                                                                                                                   <  2:33381/210‑1 (MQ=255)
                                                                                                                                                                                                                                                                 cGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGCCAACGTTTCGCGGCTATTCTTGCATTGTTTGGATTCAAGACTAGCCCCCGTGATGTTGTCa                                                                                                                                                    >  2:302658/1‑184 (MQ=11)
                                                                                                                                                                                                                                                                  gCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATTAAGCCCCCGCAACATTGCGCCTCACC‑GGAATCCTTCCGGCTTGACGCTGATTCTa                                <  2:606561/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tGTTCAGCCAGTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGCTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGc                                                                                                                                               <  1:105714/185‑1 (MQ=255)
                                                                                                                                                                                                                                                                     tGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGACTGGCCCCCGTGGTGTTGTCAGGTGc                                                                                                                                               >  2:105631/1‑185 (MQ=255)
                                                                                                                                                                                                                                                                               gTTCACGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGTTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTTTAATGACAACGTTTCGCTGCTATTCTTGTATTGTTTGGATACAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATAc                                                                                                                                           <  2:380179/179‑1 (MQ=17)
                                                                                                                                                                                                                                                                               gTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATAGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGGGCGAGGGTACCCCCTCCACAGCAACCAATACAACCAAGAAGTAAGCACACGCACCATTGCGCC‑CACCCAGAAGCCTTCCGGATTGTCACCGATTCTAAGACGGttttt                     >  2:27308/1‑296 (MQ=14)
                                                                                                                                                                                                                                                                               gTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGGTTGGATGCAAGAATGGCCCCCgtggtg                                                                                                                                                          <  2:363499/164‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    cGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTg                                                                                                                                                >  2:48820/1‑169 (MQ=255)
                                                                                                                                                                                                                                                                                    cGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTg                                                                                                                                                <  1:48820/169‑1 (MQ=255)
                                                                                                                                                                                                                                                                                    cGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGCTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTg                                                                                                                                                >  2:460866/1‑169 (MQ=255)
                                                                                                                                                                                                                                                                                    cGATCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGCTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTg                                                                                                                                                <  1:460866/169‑1 (MQ=255)
                                                                                                                                                                                                                                                                                      aTCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCCAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTAt                                                                                                                                                                                        >  1:613143/1‑127 (MQ=255)
                                                                                                                                                                                                                                                                                      aTCATACTGGCAAGAATGCCAGCAATGACCGGAGCCGCCAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTAt                                                                                                                                                                                        <  2:613143/127‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              ggCAAGAATGCCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATGGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATgg                                                                                                                                                                     <  2:645274/138‑1 (MQ=255)
                                                                                                                                                                                                                                                                                              ggCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATTAAGCCCCCGCAACATTGCGCCTCACC‑GGAATCCTTCCGGCTTGACGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGtctc  >  1:515986/1‑300 (MQ=255)
                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                   
CCAATGCCACCAGGGATAAAGCCCCCGCAACATTGCGCCTCACCGGAACCTTTCCGGCTTGCCGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGTCACGTCTTCACTCTGTGCTGCGCATTCCGCAGAATATTATCCGCTTGCAACCGTTGCATTTTTCTTACGAAAAAGGCCGCAGAGTTTCACCCCTGCGGCCTGGTTCGGGCGCAAATTGCCATTACGGCAGCCTGACGCCCGCATGGCACGTTACTTCCGCTTGTTCAGCCAGTTCACGATCATACTGGCAAGAATACCAGCAATGACCGGAGCCGCTAAATCATGCCAGAAGGCCATGCCCAGTTCTGCGAGCGTCATATAGCCGCCTGTATTGTAATGACAACGTTTCGCGGCTATTCTTGTATTGTTTGGATGCAAGAATGGCCCCCGTGGTGTTGTCAGGTGCATACCTGCAACGTGCGGGGGTTTTCTCTCTCCAGCAACCAATGCCACCAGGGATTAAGCCCCCGCAACATTGCGCCTCACC‑GGAATCCTTCCGGCTTGACGCTGATTCTACGACTGTTTTTCTGTTCCGGAATCCGTCTC  >  NC_012971/3561401‑3561986

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