Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,556,207	0	T	C	34.5%	76.8 / 25.6	55	S215G (AGT→GGT) ‡	rtcA	RNA 3'‑terminal phosphate cyclase
Reads supporting (aligned to +/- strand):  ref base T (35/1);  new base C (0/19);  total (35/20)
Fisher's exact test for biased strand distribution p-value = 7.13e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556073‑3556345                                                                                                                                         |                                                                                                                                           ggCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCg                                                                                                                                       >  1:186442/1‑140 (MQ=255)   cACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGcc                                                                                                                                     >  1:909970/1‑140 (MQ=255)        aCGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAgtgt                                                                                                                                >  1:1561942/1‑140 (MQ=255)           cgcTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGc                                                                                                                             >  2:738153/1‑140 (MQ=255)             cTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGa                                                                                                                           >  2:44323/1‑140 (MQ=255)                  ccccGACAAAAAAAAAGCGTTGGGTGATTTTTCCACTTCCGAC‑TTCAAGCGAGACGGTATTCCCCGCCCCCGGGTCGCGCGGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCa                                                                                                                      <  2:978042/138‑1 (MQ=255)                                 ggCGTTCGGTGATATTTTCACTTTCGCC‑TTCACGCGAGAGGGTGTCACCGGCCCCTTGGGCGCGGGACAGTTAATAAATTTTCTGTTCAGGGAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATgg                                                                                                       <  2:1018029/139‑1 (MQ=255)                                     ttCGGTGATATTTTCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGcgcg                                                                                                   >  1:30983/1‑140 (MQ=255)                                      tgggTGATTTTTTCACTTTCGAC‑TTAAAGCGAGACGTTATTCCCCGCCCCCTGGCCGCGCGGGAGGTTATGAATATTCTGTTCATGCAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCgg                                                                                                  <  2:1336835/138‑1 (MQ=255)                                        ggTGATTTTTTCTTTTTCGAC‑TTCAAGCGAGCCGGTATTCCCGGCCCCCGGGTCGCGGGGCAGGTTATAAATTTTCTTTTCATGAAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                <  2:763982/140‑1 (MQ=255)                                        ggTAATTTTTTCCTTTTCGAC‑TTAAAGCGAGACGTTATTCCCCGCCCCCTGGCCGCGCGGCAGTTTATAATTTTTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGca                                                                                                <  2:870372/140‑1 (MQ=255)                                            attttttCACTTTCGAC‑TTCAAGGGAGACGGTATTCCCCGCCCCCTGGCCGCGCGGCGGGTTATGAATATTCTGTTCATG‑AGGGAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:127324/137‑1 (MQ=255)                                             gttttcCACTTTCGAT‑TTAAAGCGAGACGTTATTCCCCGCCCCCTGGTCGCGGGGCAGGTTATAAATTTTCTGTCCATGGAGGGAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:894303/138‑1 (MQ=255)                                              tttttCCCTTTTGGCTTTCAGGCGAGACGGTTTTCCCGGGCCCCGGGCCGCGCGGCAGTTTAAAAATTTTGTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCa                                                                                           <  2:260999/139‑1 (MQ=255)                                                tttccttttCGCC‑TTAAAGCGAGAGGTTATTCCCCGGCCCCTGGCCGGGCGGAAGTTTAGAAATATTCTGTCCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCt                                                                                        <  2:1459557/134‑1 (MQ=255)                                                  tCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTaa                                                                                      >  1:1434320/1‑140 (MQ=255)                                                       ttCGAC‑TTAAAGGGAGCCGTTATTCCCCGCCCCCGGGCCGCGCGGCGGGTAATAAATTTTTTTTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                 <  2:1388106/140‑1 (MQ=255)                                                       ttCAAC‑TTCAAGCAGGAGGTTTTTCCCCGGCCCTTGGTGGGGGGGAAGTTTATAATTTTCCTGTTCAGGAAGAAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGaa                                                                                 <  2:405180/140‑1 (MQ=255)                                                        tCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAc                                                                                >  1:1015511/1‑140 (MQ=255)                                                         cGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAAcc                                                                               >  1:894303/1‑139 (MQ=255)                                                                aaaGGGAGAGGGTATTCCCCGGCCCCTGGTCGGGGGGCGGTTTATAAATTTTTTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                         <  2:1623275/139‑1 (MQ=255)                                                                aaaGCGAGACGTTTTTCCCCGGCCCCGGGCCGCGGGGCGGTTTATAAATTTTTTTTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCa                                                                         <  2:1015511/139‑1 (MQ=255)                                                                     gggccGGTATCCCCCGCCCCCTGGCCGCGCGGCAGTTTATGAATATTCTTTTCAGGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  2:612566/136‑1 (MQ=255)                                                                     gagaCGGTATCCCCCGGCCCCTGGCCGCGGGGCAGGTTATAAATATTCTGTTCATGAAGGAAAAAC‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  2:1561942/140‑1 (MQ=255)                                                                     gagaCGGTATCCCCCGCCCCCTGGCCGCGCGGCAGTTTATAAATATTCTGTTCA‑GGAGGGAAAAACCCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCAt                                                                    <  2:909970/140‑1 (MQ=255)                                                                        aCGGTATTCCCCGGCCCCTGGTCGCGCGGCAGTTTATGAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                 <  1:738153/140‑1 (MQ=255)                                                                        aCGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTg                                                                 >  2:56865/1‑140 (MQ=255)                                                                            tatcccccGCCCCCTGGCCGCGCGGCAGTTTATAAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGcaca                                                             <  2:30983/135‑1 (MQ=255)                                                                                 cccGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGtt                                                        >  1:284142/1‑140 (MQ=255)                                                                                  ccGGCCCGGGGTCGCGCGGCAGTTTAAGAATATTCTTTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTc                                                       <  2:1552479/140‑1 (MQ=255)                                                                                     gCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcg                                                    >  2:782534/1‑140 (MQ=255)                                                                                      ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:998546/1‑140 (MQ=255)                                                                                      ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  1:1255633/1‑140 (MQ=255)                                                                                      ccccTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCcgc                                                   >  2:240670/1‑140 (MQ=255)                                                                                         cTGGTCGCGGGGCAGTTTATAAATTTTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCt                                                <  2:5324/140‑1 (MQ=255)                                                                                           ggTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCg                                              >  2:286123/1‑140 (MQ=255)                                                                                             ccgggcgGCAGTTTATGAATTTTCTGTTCATGCAGGAAAAAA‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGcc                                            <  1:768801/136‑1 (MQ=255)                                                                                                cgcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAg                                         >  2:1043213/1‑140 (MQ=255)                                                                                                 gcgGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGt                                        >  2:1562655/1‑140 (MQ=255)                                                                                                     cAGGTTATGAATATTCTGTTAATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCa                                    <  2:284142/140‑1 (MQ=255)                                                                                                      ggtttATGAATATTCTGTTCATGAAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCTGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   <  2:1534670/137‑1 (MQ=255)                                                                                                      aaaagaTGAATATTCTGTTCAGGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCaa                                   <  2:877333/135‑1 (MQ=255)                                                                                                       gtttATGAATATTTTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAg                                  <  2:870158/138‑1 (MQ=255)                                                                                                          tataAATATTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGgtg                               <  1:190484/140‑1 (MQ=255)                                                                                                               atatTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTaaa                          >  1:1274715/1‑140 (MQ=255)                                                                                                                tatTCTGTTCATGCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAAc                         <  1:101395/140‑1 (MQ=255)                                                                                                                  ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCTGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACgc                       >  1:1534670/1‑139 (MQ=255)                                                                                                                  ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  2:1265638/1‑140 (MQ=255)                                                                                                                  ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  1:1430860/1‑140 (MQ=255)                                                                                                                  ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  1:1057715/1‑140 (MQ=255)                                                                                                                  ttCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGa                       >  2:805047/1‑140 (MQ=255)                                                                                                                   tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  2:1265586/1‑140 (MQ=255)                                                                                                                   tCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGAt                      >  2:377536/1‑140 (MQ=255)                                                                                                                     tGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGc                    >  2:1174659/1‑140 (MQ=255)                                                                                                                      gTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGCTgtc                   >  1:877333/1‑138 (MQ=255)                                                                                                                       ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  1:1240471/1‑140 (MQ=255)                                                                                                                       ttCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCa                  >  1:1405496/1‑140 (MQ=255)                                                                                                                            gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             >  1:939595/1‑140 (MQ=255)                                                                                                                            gCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             >  2:293874/1‑140 (MQ=255)                                                                                                                            gCAGGAAAAAC‑CCCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGc             <  2:1075950/140‑1 (MQ=255)                                                                                                                                 gAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc        >  1:1574519/1‑140 (MQ=255)                                                                                                                                 gAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCACCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGAc        >  1:758664/1‑140 (MQ=255)                                                                                                                                        cTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGt  >  2:1236792/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                           GGCACTGACGCGCTTTTCACCGACGACAAAAAAGCGTTCGGTGATATTTTCACTTTCGAC‑TTCAAGCGAGACGGTATTACCCGGCCCCTGGTCGCGCGGCAGGTTATGAATATTCTGTTCATGCAGGGAAAAA‑CTCCCCGCCAGTGTAGCGATTTCACGCTCAGCAACATGGCGCGGCACACCAGCTAATAGAACTTCTCCACGCATCTGCACAATGTTCCCGCGCTCGCCAAGTTGCAAGGTGTTAAACGATGCCACGGGCGAGACTTCCGT  >  NC_000913/3556073‑3556345

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