Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,926,366	0	A	C	20.7%	123.8 / 10.6	58	G30G (GGT→GGG)	mioC	flavoprotein MioC
Reads supporting (aligned to +/- strand):  ref base A (35/11);  new base C (0/12);  total (35/23)
Fisher's exact test for biased strand distribution p-value = 1.52e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

TGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACG  >  NC_000913/3926229‑3926500                                                                                                                                           |                                                                                                                                       tGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAAc                                                                                                                                       >  1:972458/1‑139 (MQ=255) tGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAAc                                                                                                                                       <  2:972458/139‑1 (MQ=255)   ttCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCg                                                                                                                                    >  1:777593/1‑140 (MQ=255)     ccTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGcc                                                                                                                                  >  2:540024/1‑140 (MQ=255)            tGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAACCCCGCCTCTTCCa                                                                                                                           <  2:636758/140‑1 (MQ=255)             gCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAg                                                                                                                          >  2:610096/1‑140 (MQ=255)             gCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAg                                                                                                                          >  1:440751/1‑140 (MQ=255)             gCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAg                                                                                                                          >  1:283537/1‑140 (MQ=255)                  aTAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCtttt                                                                                                                     <  1:649099/140‑1 (MQ=255)                   tAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTc                                                                                                                    >  1:729350/1‑140 (MQ=255)                       aaggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGCCCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGcc                                                                                                                <  1:209931/140‑1 (MQ=255)                        aggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCa                                                                                                               >  2:259607/1‑140 (MQ=255)                         ggagaaaggTTGTCCGGAATATCTCCGGCACCGTGGGGGGAGGTTATAACCAGCCAGTTCCCT‑GAGGCAGGTAATTCTTCTAACAGCGGACCGTGCAGCGTTTGGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAg                                                                                                              <  2:692957/140‑1 (MQ=255)                           agaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAggg                                                                                                            >  2:694165/1‑139 (MQ=255)                            gaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgtg                                                                                                           <  1:957780/140‑1 (MQ=255)                            gaaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGgtg                                                                                                           <  1:578850/140‑1 (MQ=255)                             aaaggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGccctgtct                                                                                                          >  1:928982/1‑134 (MQ=255)                              aaggTTGCCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGtt                                                                                                         <  1:96099/140‑1 (MQ=255)                               aggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCttt                                                                                                                      <  1:574616/126‑1 (MQ=255)                               aggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCttt                                                                                                                      >  2:574616/1‑126 (MQ=255)                                ggTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTcc                                                                                                       >  1:984972/1‑139 (MQ=255)                                    gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGcc                                                                                                                <  1:479570/127‑1 (MQ=255)                                    gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGcc                                                                                                                >  2:479570/1‑127 (MQ=255)                                    gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGcc                                                                                                                >  2:224536/1‑127 (MQ=255)                                    gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGcc                                                                                                                <  1:224536/127‑1 (MQ=255)                                    gTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTa                                                                                                   >  1:86400/1‑140 (MQ=255)                                     tCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTAc                                                                                                  <  1:100878/140‑1 (MQ=255)                                     tCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTAc                                                                                                  >  2:259869/1‑140 (MQ=255)                                       cGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACat                                                                                                >  1:150618/1‑140 (MQ=255)                                               ctcCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCa                                                                                        >  1:592835/1‑140 (MQ=255)                                               ctcCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCa                                                                                        >  1:498102/1‑140 (MQ=255)                                                 ccGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCAcc                                                                                      >  1:539593/1‑140 (MQ=255)                                                    gccccGTGGGGGGAGCTGAAAACCAGCCAGACCCCG‑GAGGCAGGAAAATCTTCAAACAGCGGCCCGTGCAGCTTTTGGGGGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCAccgcc                                                                                   <  2:592835/137‑1 (MQ=255)                                                              tggAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGGAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGggctg                                                                         >  1:334902/1‑140 (MQ=255)                                                              gggAGCTAATACCCAGCCAGATCCCT‑GAGGCAGGAAAATTTTCTAACAGCGGACCGTGCAGTTTTTCGGGGGAAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGtgctg                                                                         <  2:86400/139‑1 (MQ=255)                                                              gggAGCGAATACCCACCCAGATCCCG‑GAGGCAGGAAAATCTTCTAACAGCGCCCTGTGCAGCGTTTGGGTGAAAAACCCCGCCTCTTCCAGCTTTCCAGCCAGGTGTTCTGCTACATATTCGGCGCCGCCGAGGGtgctg                                                                         <  2:777593/139‑1 (MQ=255)                                                                gAGCTAATAACCA‑CCAAATCCCTGGAGGCGGGAAAATTTTCTAACAGCGGCCGGTGCGGCTTTGCGGTGAAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGcc                                                                       <  2:283537/140‑1 (MQ=255)                                                                 ggCTGATACCCAGCCAGCTCCCT‑GGGGCAGGAAAATTTTCTAACAGCGGACCGTGCGGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCg                                                                      <  2:440751/139‑1 (MQ=255)                                                                   cTGATACCCAGCCAGACCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGCCCGTGCGGCGTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCt                                                                    <  2:498102/140‑1 (MQ=255)                                                                      aTAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGAt                                                                 >  2:615433/1‑140 (MQ=255)                                                                             gccagATCCCG‑GAGGAAGGAAATTTTTCTAACAGCGGCCCGTCCAGCGTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTg                                                          <  1:118216/140‑1 (MQ=255)                                                                                         gAGGAAGGAAAATTTTCAAACAGCGGCCCGTGCAGCTTTTGGGTGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCAt                                               <  1:259869/140‑1 (MQ=255)                                                                                             cAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAAc                                           >  2:375066/1‑140 (MQ=255)                                                                                               ggTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCg                                         >  1:220581/1‑140 (MQ=255)                                                                                               ggTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCg                                         >  2:553698/1‑140 (MQ=255)                                                                                               ggTAAATCTTCTAACAGCGCAACGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGGGATATCTGCCATAAAccc                                         >  1:196023/1‑139 (MQ=255)                                                                                                ttAAATCTTCAAACAGCGGCCCGTGCGGGTTTTGGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                        <  2:340015/139‑1 (MQ=255)                                                                                                gTAAATCTTCTAACAGCGGCCGGTGCAGCGTTTCGGGGGAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                        <  1:828453/140‑1 (MQ=255)                                                                                                atAACTCTTCTACCGGCGCCCCGGCCGGCTTTTGGGGGATAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGc                                        <  2:196023/139‑1 (MQ=255)                                                                                                 tAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGCGATATCTGCCATAAACCGcc                                       >  1:263204/1‑140 (MQ=255)                                                                                                  aatttttcaAACGGGGGACGGGGCAGTTTTTGGGTGAAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCa                                      <  2:539593/135‑1 (MQ=255)                                                                                                         ctAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttat                               >  2:558315/1‑140 (MQ=255)                                                                                                          tAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTttatt                              >  1:757862/1‑140 (MQ=255)                                                                                                               gCGGCCCGTGCAGCTTTTGGGGGAAAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAagag                         <  2:757862/140‑1 (MQ=255)                                                                                                                cGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGt                        >  1:436431/1‑140 (MQ=255)                                                                                                                cGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGt                        >  2:329296/1‑140 (MQ=255)                                                                                                                     cGTGCGGCGTTTGGGTGGTAAACCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGt                   <  1:259607/140‑1 (MQ=255)                                                                                                                                 cGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCtgt       >  1:966562/1‑140 (MQ=255)                                                                                                                                 cGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCtgt       >  1:296859/1‑140 (MQ=255)                                                                                                                                   ggggaaaacCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGa     <  2:729350/135‑1 (MQ=255)                                                                                                                                    ggaaaaacCCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGaa    <  2:220581/136‑1 (MQ=255)                                                                                                                                      gtAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAAcg  >  2:188093/1‑140 (MQ=255)                                                                                                                                      gtAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAAcg  >  2:522005/1‑140 (MQ=255)                                                                                                                                           |                                                                                                                                       TGTTCCTGCAATGCTTCATAGAAAGGAGAAAGGTTGTCCGGAATATCTCCGGCACCGTGGGTGGAGCTGATAACCAGCCAGATCCCT‑GAGGCAGGTAAATCTTCTAACAGCGGACCGTGCAGCGTTTCGGTGGTAAAACCCGCCTCTTCCAGCTTTTCAGCCAGGTGTTCTGCTACATATTCGGCACCGCCGAGGGTGCTGCCGCTGATAAGAGTGATATCTGCCATAAACCGCCACCTTTATTAAGAGTGGCGTATTGTACGCTGTGAACG  >  NC_000913/3926229‑3926500

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