Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2070TC23.7% 145.4 / 10.6 60S215G (AGT→GGT) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (42/3);  new base C (0/14);  total (43/17)
Fisher's exact test for biased strand distribution p-value = 5.11e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGC  >  NC_000913/3555958‑3556456
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
ccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTc                                                                                                                                                                                                                                                           >  1:552597/1‑251 (MQ=255)
    gcccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCGCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCTGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑g                                                                                                                                                                                                                                                      <  1:200085/251‑1 (MQ=255)
    ccccgccAGCGCCATCGGTATCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCCTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCCCTTTTCTCCTCCTTCCTTTTTTCCTTCTCTCATATTTTCACTTTCGACTTCAAGCGATACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑g                                                                                                                                                                                                                                                      <  1:449625/250‑1 (MQ=255)
      ccgccAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑Gcc                                                                                                                                                                                                                                                    <  1:469756/251‑1 (MQ=255)
       cgccAGCGCCATCGGTAGCACCCGTTGGTCAGCACGATATTCCCCCACCGCCGCCGCGCTTGCCAGGCAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTAGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGCTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑Gcc                                                                                                                                                                                                                                                    >  2:50337/1‑250 (MQ=255)
                     gTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTc                                                                                                                                                                                                                                      >  1:513656/1‑250 (MQ=255)
                             aGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGc                                                                                                                                                                                                                             >  1:320903/1‑251 (MQ=255)
                               ttGGTCAGCAAGTTATCCCCCCCCCGCCCCCGTGCTTGCCAGGTAGCGTTCCCCCTCTTTCACCAACTGGCCCCCACCCCCCTCGGCACTGCCGCGTTTTCCACCGACAAAAAAAAAGCGTTGGGTGATTTTTTCTCTTCCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGG‑TCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCa                                                                                                                                                                                                                            <  1:176055/250‑1 (MQ=255)
                                           atatTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGc                                                                                                                                                                                                                >  1:477538/1‑250 (MQ=255)
                                                             cgTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACtt                                                                                                                                                                                             >  2:122823/1‑251 (MQ=255)
                                                              gTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTtc                                                                                                                                                                                            >  1:317012/1‑251 (MQ=255)
                                                                       aGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAACAACTCCCC‑GCCAGTGTAGCGATTTCATGCTCAGCAACATGGCGCGGCTCACCAGCTAATAGAACTGCTCCACGCAt                                                                                                                                                                                   >  2:147776/1‑251 (MQ=255)
                                                                        ggTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCGCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACCCCAGCTAATAGAACTACTCCACGCATc                                                                                                                                                                                  >  2:8804/1‑251 (MQ=255)
                                                                                ttCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAAt                                                                                                                                                                          >  2:255564/1‑251 (MQ=255)
                                                                                   aCCTCTTTCACCAACTGTGCCGCAACAACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                                                                                                                                       >  1:379079/1‑251 (MQ=255)
                                                                                                   gtgCCGCACCCACCTCGGCACTGACGCGCTTTTCACAGACGACAAAAAAGCGTTCGGGGAGATTTTCACTTTCGACTTCAAGCGAGACGACATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATAAATATTCTGTACATGCAGGGAAAAACACCCC‑GCCAGTGTAGCGATTTAAAGCACAGCAACATGGCGCGGCACACCAGCTAATAGCACTTCACCACACCACTGCACTATATTCCAGCGCTCGCCAAgg                                                                                                                                                       >  2:515690/1‑250 (MQ=255)
                                                                                                               ccTCGGCACTGCCGCTCTTTCCACCGACGAAAAAAAAGCGTTCGTTGATTTTTTCACTTTCGACTTAAAGCGAGACGGTATTCCCCGGCCCCGGG‑TCGCGCGGCAGGTTATAAATTTTCTGTTCATGCAGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTa                                                                                                                                           <  1:550843/251‑1 (MQ=255)
                                                                                                                     cACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                                                                                                                                     >  1:376754/1‑251 (MQ=255)
                                                                                                                            gcgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGCATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACgg                                                                                                                              >  1:477555/1‑251 (MQ=255)
                                                                                                                              gcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc                                                                                                                            >  2:91267/1‑251 (MQ=255)
                                                                                                                                     aCCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGACGCCACGGGCGAGACtt                                                                                                                     >  2:621281/1‑251 (MQ=255)
                                                                                                                                      ccGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCAAGGGCGAGACTTc                                                                                                                    >  2:475570/1‑251 (MQ=255)
                                                                                                                                       cgacgaCAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc                                                                                                                   >  1:460185/1‑251 (MQ=255)
                                                                                                                                       cgacgaAAAAAATGCGTGTGTTGATATTTTCCCTTCCGCCTAAAAGCGAGAGGTTATTCCCCGCCCCCTGG‑CCGCGCGCCAGGTTATGAATATTCTGTCCATGCAGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTcc                                                                                                                   <  1:122823/251‑1 (MQ=255)
                                                                                                                                          cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTCCCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTTTTCATGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGtt                                                                                                                <  1:209451/251‑1 (MQ=255)
                                                                                                                                          cgacAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGCTTTCACGCTCAGCAACATGGCGCGGCACACCAGCTACTAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGCGACTTCCGtt                                                                                                                >  2:210663/1‑251 (MQ=255)
                                                                                                                                             aaaaaaaGCTTTGGGTGATTTTTTCCCTTTCGACTTCAAGCGAGACGTTATT‑CCCCGCCCCTGGGTCGCGCGGCAGGTTATAAATTTTCTTTTCATGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcc                                                                                                             <  1:475570/250‑1 (MQ=255)
                                                                                                                                                    gCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCacac                                                                                                      >  1:185652/1‑251 (MQ=255)
                                                                                                                                                      gTTCGGGGATATTTTCCCTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCg                                                                                                    >  2:507073/1‑251 (MQ=255)
                                                                                                                                                       ttCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGc                                                                                                   >  1:136531/1‑251 (MQ=255)
                                                                                                                                                            tGATATTTTCACTTTCGACTTCAAGCGCGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGATATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCACGGGGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCg                                                                                              >  1:434898/1‑251 (MQ=255)
                                                                                                                                                               tatTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGc                                                                                            >  1:579361/1‑250 (MQ=255)
                                                                                                                                                                     cACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTa                                                                                      >  2:557124/1‑250 (MQ=255)
                                                                                                                                                                         ttCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAAc                                                                                 >  2:41700/1‑251 (MQ=255)
                                                                                                                                                                           caccgtCAAGCGAGGCGTTCTTCCCCGGCCCCTGG‑TCGCGCGGGGGGTTAGAAATTTGCGGTCCAGGC‑GGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCg                                                                               <  2:524100/246‑1 (MQ=255)
                                                                                                                                                                            gACTTCCAGCGAGACGGAATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCACGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGGGTTAAACGATGCCACGGGCGAGACTTCCGTGGCCACCACACCGCCTCCGGCACGGTAAAAACCGt                                                                              >  2:599295/1‑251 (MQ=255)
                                                                                                                                                                              cTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTgg                                                                            >  1:211827/1‑251 (MQ=255)
                                                                                                                                                                               ttCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGACTATTCAGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCCAGTTGCAAGGAGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCCGGGTAAAAACCGTggg                                                                           >  2:132908/1‑250 (MQ=255)
                                                                                                                                                                                     ggagggggTATTCCCCGGCCCTGGG‑GCGCGCGGCAGGTTAGGAATTTTCTTTTCATGCAGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACa                                                                     <  1:91267/245‑1 (MQ=255)
                                                                                                                                                                                      gagaCTGTATTCCCCGCCCCCGGG‑CCGCGCGGATGTTTATAAATTTTTTTTTAATGCAGGGAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAg                                                                    <  1:86925/251‑1 (MQ=255)
                                                                                                                                                                                        gaCGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTACCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCg                                                                  >  1:274280/1‑251 (MQ=255)
                                                                                                                                                                                          cGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTg                                                                >  1:201105/1‑251 (MQ=255)
                                                                                                                                                                                             tATTCCCCGCCCCCTGG‑TCGCGCGGCAGGTTATAAATATTCTGTTCATGC‑GGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGt                                                              <  1:209978/250‑1 (MQ=255)
                                                                                                                                                                                                 ccccGGCCCCTGG‑CCGCGCGGCAGGTTATGAATATTCTTTGCAGGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGc                                                          <  1:503825/249‑1 (MQ=255)
                                                                                                                                                                                                 aCCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCt                                                         >  1:35810/1‑251 (MQ=255)
                                                                                                                                                                                                 aCCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCt                                                         >  1:413864/1‑251 (MQ=255)
                                                                                                                                                                                                 aCCCGGCCCCTGG‑TCGCGCGGCAGCTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCACGTTGCAAGGTGTTAAAAGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCCGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCt                                                         >  2:324188/1‑251 (MQ=255)
                                                                                                                                                                                                  cccGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCtg                                                        >  1:93053/1‑251 (MQ=255)
                                                                                                                                                                                                    cGGCCCCTGT‑TCGCGGGGCAGTTTATAAATATTCTTTACAGGAAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCtgat                                                      <  1:16641/251‑1 (MQ=255)
                                                                                                                                                                                                     ggCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGATTGCtgatg                                                     >  2:96871/1‑251 (MQ=255)
                                                                                                                                                                                                                         gtatgtcatattcgtgatgaTGCAGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCGGCAACATGGCGCGGCACACCAGCTAATAGAATTTCTCCACGCACCTGCACAAGGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCCCGGTCGAGACTTCCGTTGCCACCACCCCGCCTCCGGCAGGGTAAAAACCGTGGTGGAACAGCGTGGGTTGCTGATGAATTCCTATTTTCGCcagc                                  <  2:212831/232‑1 (MQ=255)
                                                                                                                                                                                                                            tATGAATATTCTGTTCATGAAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCcagcag                                <  1:217286/250‑1 (MQ=255)
                                                                                                                                                                                                                             ataaatgttcttgcaacggggaaaaaaaCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGATCTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCg                              <  2:274280/228‑1 (MQ=255)
                                                                                                                                                                                                                                        gTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCAcc                   >  2:601828/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCcgg                >  1:617903/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGAGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCcgc                >  2:91866/1‑250 (MQ=255)
                                                                                                                                                                                                                                               cAGGGAAACACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGGGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTACTATTTTCGCCAGCAGCGTCTCCAGAACCcggcgg             >  2:589863/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   gAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATaa         >  1:603671/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   gAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTACAGCACCCGGCGGATaa         >  1:210640/1‑250 (MQ=255)
                                                                                                                                                                                                                                                   aaaaccatCCCG‑GCCAGTGTTGCGATTTCACGCTCAGCCACATGTCGCGGCACGCCAGCTAATAGAACTGCTCCACGCCGTTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGTGCGAGACTTCCTTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAAATCCTAGTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGAGaaa        <  2:47699/244‑1 (MQ=255)
                                                                                                                                                                                                                                                         cTCCCC‑GCCAGTGTAGCGATTTCACGCTCCGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGc  >  1:133279/1‑251 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                          
CCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGACTTCAAGCGAGACGGTATTACCCGGCCCCTGG‑TCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGC  >  NC_000913/3555958‑3556456

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