Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,556,2000GA27.8% 59.5 / 11.8 37S217F (TCC→TTC) rtcARNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base G (22/4);  new base A (0/10);  total (23/14)
Fisher's exact test for biased strand distribution p-value = 3.94e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGA  >  NC_000913/3555951‑3556446
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
aatccccccGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGa                                                                                                                                                                                                                                                        <  2:517922/247‑1 (MQ=255)
                                              aaaGATATCCCCCCCCCCCCGCGTTCCTTGCCGGGTAGCGTTTCCCCTCTTCCACCAACGGGGCCCCACCCCCCTCGGCACTGCCGCGCTTTCCCCCGACAAAAAAAAAGCGTTCGGTAATATTTTCCCTTC‑C‑GACTCAAAGCGAGCCGTTATTCCCCGGCCCCTGGCCGCGCGGCGGGTTATGAATATTCTGTTCATGCAGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                                                                                                                         <  1:256187/250‑1 (MQ=255)
                                               aaaatttcccccccccGCCGCCGTGCTTGCCAGGTGCCTTTTCCCCTCTTCCCCCACCTGTGCCGCACCCCCCTCGCCACTGCCGCGCTTTCCCCCGACGAAAAAAAAGCGTTGGTTGATTTTTTCACTTT‑C‑GCCTTCAAGCGAGACGGTTTTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCAGGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                                                                                                                        <  1:6690/243‑1 (MQ=255)
                                                          cccccccccccttGCTTCCCAGGTGGCTTTCCCCCTCTTCCCCCAACTGTCCCGCACCCCCCTCGCCCCTCCCGCCCTTTCCCCCGCCGACAAAAAAGCTTTCGGTGATTTTTTCCCTTT‑C‑GCCTTCAAGCGAGCCGGTATTCCCCGCCCCCTGGTCGCGCGGCGGGTTATAAATTTTCTTTTCATGCAGGAAAAACCTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                                                                                                             <  1:83194/239‑1 (MQ=255)
                                                          ccACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                                                                                                                             >  2:91678/1‑251 (MQ=255)
                                                                                      tttCCCCTCTTCCCCCAACTGCGCCGCACCCCCCGCGGCACTGACGCGTTTTTCACCGACGACAAAAAAGCTTTCGGTGATATTTCCACTTC‑C‑GACTTCAAGCGAGCCGTTTTTCCCCGGCCCCGGGGCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcac                                                                                                                                                                   <  1:99700/249‑1 (MQ=255)
                                                                                             tctTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                                                                                                                          >  2:535790/1‑251 (MQ=255)
                                                                                              ctTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                                                                                                                         >  1:157365/1‑251 (MQ=255)
                                                                                                  cACCAACTGTGCCGCACCCCCCCCGGCACTGACGCGTTTTCCACCAACGAAAAAAAGGCTTTGGGTGATATTTTCCCTTC‑C‑GATTTCAAGCGAGACGTTATTCCCCGCCCCCTGGTCGCGGGGAAGGTTATAAATATTCTGTTCATGCAGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCt                                                                                                                                                     <  1:535790/251‑1 (MQ=255)
                                                                                                  cACCAACTGGCCCGCACCCCCCTCGGCACTGCCCCTCTTTCCCCCGACGAAAAAAAAGCTTTCGGTAATTTTTTCCCTTT‑C‑GACTTCAGGCGAGACGGTATTCCCCGCCCCCTGGTCGCGCGGCAGGTTATGAATATTCTTTTCATGCTGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCt                                                                                                                                                     <  1:428454/251‑1 (MQ=255)
                                                                                                    ccAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                                                                                                                                   >  1:490726/1‑251 (MQ=255)
                                                                                                        cTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCa                                                                                                                                                >  1:533575/1‑250 (MQ=255)
                                                                                                           tgCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                                                                                                                            >  1:215582/1‑251 (MQ=255)
                                                                                                           ggCCGCACCCCCCCCGGCACTGCCGCGCTTTCCCCCGACGAAAAAAAAGCGTTCGTTGATATTTTCCCTTC‑C‑GCCTACAGGGGAGGCTGTATTCCCCGCCCCCGGGTCGCGCGGCAGTTTATAAATTTTCTGTTCATGC‑GGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                                                                                                                            <  1:265300/250‑1 (MQ=255)
                                                                                                           ggCCGCAACCCCCGCGCCCCTGACGCGTTTTTCACCAACGAAAAAAAAGCGTTCGGTGATATTTTC‑CTTT‑CCGACTTCAAGCGAGACGGTATTCCCCGGCCCGGGGTCGCGCGGCAGGTTATGAATTTTCTTTTCTTGCAGGGAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGtt                                                                                                                                            <  1:513814/250‑1 (MQ=255)
                                                                                                             ccGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGgcgc                                                                                                                                                                                                         >  2:129055/1‑188 (MQ=255)
                                                                                                                                                acgaaAAAAAGGCTTTCGGTGATATTTTCTCTTT‑C‑GACTACAAGCGAGACGGTATTCCCCGGCCCCGGGTCGCGCGGAAGGTTATAAATATTCTGTTCAGGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt                                                                                                       <  1:72901/251‑1 (MQ=255)
                                                                                                                                                  nacAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTg                                                                                                     >  2:255590/2‑251 (MQ=255)
                                                                                                                                                     aaaaaaGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGcc                                                                                                   >  1:553954/1‑250 (MQ=255)
                                                                                                                                                      aaaaaGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                                                                                                                                 >  1:185847/1‑203 (MQ=255)
                                                                                                                                                        aaaGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTCAAACGATGCCACGGGCGAGACTTCCGTTGccacc                                                                                                >  1:527227/1‑250 (MQ=255)
                                                                                                                                                            cGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACAcc                                                                                           >  1:81309/1‑251 (MQ=255)
                                                                                                                                                               tCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGcc                                                                                        >  2:468031/1‑251 (MQ=255)
                                                                                                                                                                gggTGATTTTTTC‑CTTTGC‑GACTACAAGGGAGACGATATTCCCCGGCCCGTGGCCGCGCGGCAGTTTATGAATTTTCTTTTCAGGCAGGAAAAAACCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCt                                                                                       <  1:91678/250‑1 (MQ=255)
                                                                                                                                                                cGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCt                                                                                       >  2:530777/1‑251 (MQ=255)
                                                                                                                                                                     atatTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                                                                                                                         >  2:286488/1‑180 (MQ=255)
                                                                                                                                                                     atatTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGc                                                                                  >  1:70834/1‑251 (MQ=255)
                                                                                                                                                                       atTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGc                                                                                  >  1:235570/1‑249 (MQ=255)
                                                                                                                                                                                           aaGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAACTAATAGAACTTCTCCACGCATCTGCACAATGTTccc                                                                                                                                                          >  1:344092/1‑159 (MQ=255)
                                                                                                                                                                                                agaCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGc                                                         >  2:551913/1‑251 (MQ=255)
                                                                                                                                                                                                     gTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGt                                                    >  2:57883/1‑251 (MQ=255)
                                                                                                                                                                                                                    gggccgccgggtgttttagaaaTATTGTGTCCATGCGGGGAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAAtt                                       <  2:95410/229‑1 (MQ=255)
                                                                                                                                                                                                                     ggTCGCGGGGCAGGTTATAAATTTTTTGTTCATCGAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCt                                    <  1:212701/251‑1 (MQ=255)
                                                                                                                                                                                                                                aGGTTATGAATATTCTGTTCATGCAGGCAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCTGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCACGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACAGCCTCCGGCAAGGTAAAAACCGTGGCGTACCAGCGCGGTTTGCTGATGAAGTCCTATCTTcacc                           >  2:535552/1‑246 (MQ=255)
                                                                                                                                                                                                                                  gttataatttttttGTTCAGGCAGGAAAAACCCCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCcagca                       <  1:468031/238‑1 (MQ=255)
                                                                                                                                                                                                                                                     tGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACATTGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCcggcg    >  1:519239/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        aGGGAAAAACTCCCC‑GACAGGGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACACCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCcggcggc  >  2:114680/1‑249 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       
AATTCCCCCGCGCCCGCCAGCGCCATCGGTAGCACCAGTTGGTCAGCAAGATATTCCCCCACCGCCGCCGTGCTTGCCAGGTAGCGTTTCACCTCTTTCACCAACTGTGCCGCAACCACCTCGGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTT‑C‑GACTTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAACTCCCC‑GCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGTTGCCACCACACCGCCTCCGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGA  >  NC_000913/3555951‑3556446

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