Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,396,7720TC76.9% 66.3 / 14.7 40F236F (TTT→TTC) amiBN‑acetylmuramoyl‑L‑alanine amidase B
Reads supporting (aligned to +/- strand):  ref base T (4/5);  new base C (16/14);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGA  >  NC_000913/4396528‑4397016
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
aaCCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGc                                                                                                                                                                                                                                                <  1:104191/251‑1 (MQ=255)
 aCCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCg                                                                                                                                                                                                                                               <  1:241975/251‑1 (MQ=255)
      cgcgcAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCCTCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGttttt                                                                                                                                                                                                                                          >  2:87239/1‑250 (MQ=255)
       gcgcAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTaa                                                                                                                                                                                                                                         <  1:41779/251‑1 (MQ=255)
          cAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAAccc                                                                                                                                                                                                                                      >  2:59725/1‑251 (MQ=255)
            aTCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGTTCAATGACGAACCGATGTTCAAAAGCGTTTTAACCCgg                                                                                                                                                                                                                                    >  2:244754/1‑250 (MQ=255)
             tacGTTTAAAACGGAAAGTAACAGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTCCCGGGGGGAAAAAAGACACCATCGCCATCGCGCGTAAATTGCCTACCTTCCTCAATGACCATCCGATGTTCAAACGCGTTTTCACCCGTg                                                                                                                                                                                                                                   >  2:161106/3‑251 (MQ=255)
               cGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGATAAATGTCACCATCGCCATCGCGCGTAAATTGCGTATTTTGCTCAATGTGGATCCGATGTTCAAAGGCGTTTTAACCCGTGAc                                                                                                                                                                                                                                 <  1:246689/251‑1 (MQ=255)
               agTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTCACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGCCCCTTGCGCTATTGGCCCCGGTAGGACGCGGGGAAAAAATGTCCCCATAGCCATCGCGCGGAAATTTCGTACTTTTATCACTGATGATCCGATGTTCAAAGGCGTGTTAAACCGTgg                                                                                                                                                                                                                                  >  1:9962/2‑249 (MQ=255)
                    aaaaCGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACgggg                                                                                                                                                                                                                             <  1:232641/250‑1 (MQ=255)
                               tAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCCTGACGGGGATTACTTTAtct                                                                                                                                                                                                                 >  2:118181/1‑251 (MQ=255)
                                     cACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGa                                                                                                                                                                                                           <  1:87239/251‑1 (MQ=255)
                                         aCGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCAGCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATg                                                                                                                                                                                                         <  1:217807/249‑1 (MQ=255)
                                           gggTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCCATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTTTCTCGGTGATGGGgc                                                                                                                                                                                                     >  2:216464/1‑251 (MQ=255)
                                                 tATCAGAAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGGGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGCTACGGGGATTACTTTATCTCGGTGATGGGGCGCCGTg                                                                                                                                                                                               >  2:179015/1‑251 (MQ=255)
                                                               cGGTTACGCGTCTGACTGCGCGAGCGACGGCTGACACTGGCGATAAAATTATCATCGCTATCGATGCCGGACACGGCGGTCA‑TGACCCTGGCGCTATCGCCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACGTTGCTCAATGACTATCCGATGTTTAAAGGCGTTTGAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTaat                                                                                                                                                                                 <  2:68016/251‑2 (MQ=255)
                                                                        gTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGcc                                                                                                                                                                        <  1:14661/251‑1 (MQ=255)
                                                                        gTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGAGCTATCGGCCCCGGCGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAAACCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGcc                                                                                                                                                                        >  1:179879/1‑251 (MQ=255)
                                                                                   cgcgcgACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCgt                                                                                                                                                             <  2:223252/251‑1 (MQ=255)
                                                                                                            aaaTTATCATCGCTATTGATGCCGGACACGGCGGTCAAGAACACCGG‑GCCATCGGCCCCGGTGGTACGCGGGGGAAAACTGACACCATCGCCATCGCGCGTAAATTGCATATTTTGCTCAATGACGAAACGACGGTGAAAAGCGGTTTAACCCGTGACGGGGATTAATTTATCTCGGGGATGGGGCGCCGTGCTGTGGCACATCAGCAAACCGCCAATTTCCTCGTGTCGATTCAAGATGAAGCCCCACCg                                                                                                                                    >  1:267445/1‑251 (MQ=255)
                                                                                                             aaTTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGCTGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCg                                                                                                                                    >  1:158469/1‑250 (MQ=255)
                                                                                                                   catcGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATAGGCCCCGGTGGGACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTACTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGGGATGGGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGc                                                                                                                              >  1:177719/1‑250 (MQ=255)
                                                                                                                          aTTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCCTCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACCCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCCCAGTGCGAc                                                                                                                      >  2:151860/1‑251 (MQ=255)
                                                                                                                                gCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCAGCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGGGATGTGGCACTTAAGCAAAGCTCCAATTTCCTCGTGTCGATGCACGGTGATGCCTCAGCGAACCGCAGTGCGACTGgcg                                                                                                                 <  2:263802/250‑1 (MQ=255)
                                                                                                                                           cggTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATggg                                                                                                     <  1:197029/251‑1 (MQ=255)
                                                                                                                                            ggTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGt                                                                                                    <  1:42238/251‑1 (MQ=255)
                                                                                                                                                         ggCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAAcc                                                                                         <  1:183268/250‑1 (MQ=255)
                                                                                                                                                                  ggCCACGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCac                                                                               >  2:62520/1‑250 (MQ=255)
                                                                                                                                                                             tACGCGGGAGAAAAATGTCACCCTCGCAATCGCGCGTAAATTGCGTAATTTGCTCAATGACCCTCCGATGTTCCAAGCCGTTTTAACCCGTGACGGGGATTACTTTCTCTCGGAGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCTATTTCCTCGTGTCGATTAACGCTGATGCCGCACCGCACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTACCAGTCGCGCCAACAGTGAAATgg                                                                    >  2:249089/1‑251 (MQ=255)
                                                                                                                                                                                      gaAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAgctggc                                                           <  1:216464/251‑1 (MQ=255)
                                                                                                                                                                                           aTGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGag                                                      <  1:118181/251‑1 (MQ=255)
                                                                                                                                                                                             gTCACCATCGCCATCGCGCGTAAATTGCGTACGTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGagc                                                     <  2:234923/250‑1 (MQ=255)
                                                                                                                                                                                                aCCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGagcagc                                                  >  2:3356/1‑250 (MQ=255)
                                                                                                                                                                                                aCCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCAGATGTTTACAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCCACAGTGAAATGGCCAGCTGGCTGGagcagc                                                  >  1:187196/1‑250 (MQ=255)
                                                                                                                                                                                                        catcgcGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAAc                                         <  1:45480/251‑1 (MQ=255)
                                                                                                                                                                                                          tcgcGCGTAAATTGCGTCCTTTGCTCAATGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAAc                                         >  1:220091/1‑249 (MQ=255)
                                                                                                                                                                                                                                    aaTGACGATCCGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCgggggg             >  2:216459/1‑248 (MQ=255)
                                                                                                                                                                                                                                             ccGATGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGAGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCCCGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACCGGGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGGGGGGCGGGGGAtgagcaggc    >  2:258133/1‑244 (MQ=255)
                                                                                                                                                                                                                                             ccGATGTGCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGTCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGAATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAAGGGCCAGCTGGCTGGAGCAGCACGAGAAACAGGCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGc    <  1:249089/251‑1 (MQ=255)
                                                                                                                                                                                                                                                aTGTTCAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGGGATGTGCTGGCGa  >  1:52640/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
AACCGGCGCGCAATCCGTTTAAAACGGAAAGTAACCGCACTACGGGTGTTATCAGCAGTAATACGGTAACGCGTCCGGCAGCGCGCGCGACGGCTAACACTGGCGATAAAATTATCATCGCTATTGATGCCGGACACGGCGGTCA‑GGACCCTGGCGCTATCGGCCCCGGTGGTACGCGGGAGAAAAATGTCACCATCGCCATCGCGCGTAAATTGCGTACTTTGCTCAATGACGATCCGATGTTTAAAGGCGTTTTAACCCGTGACGGGGATTACTTTATCTCGGTGATGGGGCGCAGTGATGTGGCACGTAAGCAAAACGCCAATTTCCTCGTGTCGATTCACGCTGATGCCGCACCGAACCGCAGTGCGACTGGCGCTTCCGTATGGGTGCTCTCTAACCGTCGCGCCAACAGTGAAATGGCCAGCTGGCTGGAGCAGCACGAGAAACAGTCGGAGCTGCTGGGTGGGGCGGGTGATGTGCTGGCGA  >  NC_000913/4396528‑4397016

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: ‑