Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,718	0	G	C	43.3%	87.6 / 32.1	69	P292R (CCG→CGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base G (12/26);  new base C (29/0);  total (43/26)
Fisher's exact test for biased strand distribution p-value = 1.03e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCT  >  NC_000913/2984568‑2984867                                                                                                                                                       |                                                                                                                                                       cccccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCg                                                                                                                                                        <  2:1228974/148‑1 (MQ=255) ccaccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCg                                                                                                                                                        <  2:463357/151‑1 (MQ=255)    ccaTTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGgcg                                                                                                                                                     <  1:1331300/151‑1 (MQ=255)         gacattgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAAtt                                                                                                                                               <  1:515618/151‑1 (MQ=255)            attgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACAACAGCG‑CCAATccc                                                                                                                                            >  2:258521/1‑151 (MQ=255)             ttgacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGCTCGCCTACAACGGCG‑CCATTccc                                                                                                                                            >  2:1047247/1‑150 (MQ=255)               gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCG‑CCATTccc                                                                                                                                            >  1:1029466/1‑148 (MQ=255)               gacCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCAt                                                                                                                                          <  1:1039494/150‑1 (MQ=255)                      cGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCAGACACCAGCG‑CAAATCCCATCAATgac                                                                                                                                   >  2:414183/1‑148 (MQ=255)                       gCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCCTCacttcc                                                                                                                                   >  2:1342561/1‑144 (MQ=255)                          cATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGa                                                                                                                                <  2:893439/149‑1 (MQ=255)                           aTCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAATTCCCATCAAAGCCGGGTCt                                                                                                                             >  1:125031/1‑151 (MQ=255)                              ccACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTAACCCAACTGCCAGCCTACACGCTCCAGGCAACGCCGGGTCGCATACACCGGCG‑CAATTTCCAACAATACCCGAACTa                                                                                                                            >  1:792226/1‑149 (MQ=255)                                aCTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACAcc                                                                                                                        <  1:587205/151‑1 (MQ=255)                                     aGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCCTCCAAGCCGGATCTACACCaa                                                                                                                      >  1:1190187/1‑148 (MQ=255)                                        aTCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACACCCGCG‑CCAATTCCACCCATGCCGGATCTACACCAACGCTgg                                                                                                                >  2:782000/1‑151 (MQ=255)                                          cTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAATTTCCCTCAATGCCGGGTCTACACCCACCCTgg                                                                                                                >  2:601999/1‑149 (MQ=255)                                           ttGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CCATTCCCATCAATAGCGGATCTACAACAAACCTgg                                                                                                                >  1:1275269/1‑148 (MQ=255)                                            tGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGc                                                                                                               <  2:927256/148‑1 (MQ=255)                                            tGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCa                                                                                                              <  1:44371/149‑1 (MQ=255)                                               cAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAAATCCCATCAATGCCGGGTCTACACCCACGCTGGCAAATg                                                                                                          >  1:471303/1‑150 (MQ=255)                                               cAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCACCCACCCGCG‑GAAATCCCATCAATGGCCGGTCAACACCAACCCTGGCaaa                                                                                                            >  1:346867/1‑148 (MQ=255)                                                  ccGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATgcgc                                                                                                       <  2:333608/150‑1 (MQ=255)                                                   cGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCGCG‑CCAATCCCCTCAATGCCCGATCTACACCAACGCTGGGAAAAGcgcg                                                                                                      >  1:723231/1‑150 (MQ=255)                                                       aGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGc                                                                                                 <  2:777074/151‑1 (MQ=255)                                                               tgCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCCGCG‑CAAATCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGcc                                                                                            >  1:1119517/1‑148 (MQ=255)                                                                     aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGt                                                                                   <  1:396868/151‑1 (MQ=255)                                                                     aaaCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCG‑CGATTTCCATCAAATCCGGGTCTTCACCAAACCTGGCAAAAGCGCGGAAGCGGGCCACCACcgg                                                                                    >  1:374428/1‑150 (MQ=255)                                                                     aaaCGCTTCATTAGCCCCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATTCCCCGGCG‑CAATTCCCCTCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGt                                                                                   <  1:1198377/151‑1 (MQ=255)                                                                       aCGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACAACCGCG‑CAATTCCCCTCAATACCGGATCTACACCCACCCTGGGCAAAGCGCGGAATCGGGGCAGCACggg                                                                                    >  1:283067/1‑148 (MQ=255)                                                                         gCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCCCATAAACCCGCG‑CCATTCCCCTCAATGCCGGGTCTACACCAACCCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGaaa                                                                                >  1:7603/1‑150 (MQ=255)                                                                             cATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCCTACAACCGCG‑CAATTCCCATCAATGCCGGATCTCCACCCACACTGGCAAAAGCGCGAAAGCGGGCCCACACgggggaat                                                                               >  2:810376/1‑142 (MQ=255)                                                                              aTTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAAc                                                                          <  2:283062/151‑1 (MQ=255)                                                                              aTTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTTCCCTCAAAGCCCGATCTCCACCAACCCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAATTCaa                                                                           >  2:1191972/1‑150 (MQ=255)                                                                                tAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACAACCGCG‑CCATTCCCATCAAAGCCCGGTCTACACCAACGCCGGCAAAAGCGCGAATGCGGGGCAGCACGGGGAAAATCAAAGc                                                                        >  2:1351368/1‑151 (MQ=255)                                                                                tAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTAAACCCGCG‑CAATTCCCACCAATGCCGGAACTTCACCAAAACCGGCGAAATCGCGAATGCGGGCCCACAAGGGGAAATTTAAAGc                                                                        >  1:1098985/1‑151 (MQ=255)                                                                                 aGCCTCGATAAGATCGACTTCAGCCGACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CACTTCCCATCAATGCCCGATCTACACCAACGCCGGGAAAAGCGCGAATGCGGGGCAGCACGGGGAAATTCAACGcc                                                                       >  2:17432/1‑150 (MQ=255)                                                                                 aGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTCCACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAAc                                                                          <  2:1194670/148‑1 (MQ=255)                                                                                 aGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACCCCCGCG‑CAATTCCCATCAATGCCGGATCTACACCAACACTGGCAAATGCGCGAATGCGGGCCAGCACGGGGCAATTTaa                                                                           >  1:1082056/1‑147 (MQ=255)                                                                                 aGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCG‑CAAATCCCATCAATTCCGGATCTACACCAAAGCTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAATTCaa                                                                           >  2:1166881/1‑147 (MQ=255)                                                                                     tCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CAATTCCCATCAATGCCCGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCCGCACGGGGAAATTCAACGCTCGTg                                                                   >  1:535841/1‑151 (MQ=255)                                                                                     tCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACAACCCCG‑CAATTCCCACCCATCCCCGATCTACACCCACCCCGGCCAAAGGCCGAATGCGGGCCCGCCCGGGTAAATTTaa                                                                           >  2:936372/1‑143 (MQ=255)                                                                                           aGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACAACCGCG‑CAAATCCCCCCAATACCCGATCTTCCCCAACGCCGGCAAAAGCGCGCAAGCGGGCCAGCACCGGGAAAATTAAAACCCGTGCttt                                                               >  2:983543/1‑149 (MQ=255)                                                                                           aGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGTGCATACACCGGCG‑CAAATCCCCTCAATGCCGGAACTACACCAACGCCGGCAAATACGCGGAAGCGGGCCAGCCCCGGGAAATTAAAAGCTCCTGCttt                                                               >  2:429135/1‑149 (MQ=255)                                                                                            gATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCCTACACCCCCGCCAATTCCC‑TCAAAGCCCGGGCTACACCAACCCCGGCCAAAGCGCGAATGCCGGGCAGCACCGGGAAATTCAAcaccccggcttt                                                               >  1:681069/1‑137 (MQ=255)                                                                                                    tCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtc                                                     <  1:721098/150‑1 (MQ=255)                                                                                                     cAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtca                                                    <  1:1142723/150‑1 (MQ=255)                                                                                                      aGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGGCGCATACACCGGCG‑CAATTTCCATCAATGCCGGATCTACACCAACGCCGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAATTCAACGCTCGTGCTTTTGCTTCGCtc                                                     >  1:683688/1‑148 (MQ=255)                                                                                                        ccATCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCtcatca                                                 <  2:736292/150‑1 (MQ=255)                                                                                                           aCTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACAACCGCG‑CAAATCCCCTCAATACCGGATCTACACCAACACTGGCAAAAGCGCGAATGCGGGCCAGCACGGGGAAATTTAAACCTCGTGCTTTGGCTTCGCtca                                                    >  2:638892/1‑144 (MQ=255)                                                                                                            cTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATAAACCCGCG‑CCATTCCCATCCAAGCCGGAACTACACCAACACCGGCCAATGAGCGGAAGGGGGGCAGCACCGGGAAATTCAACGCTCGTGCTTTTGCTTCGCccatcatcact                                            >  1:742670/1‑149 (MQ=255)                                                                                                                cAGCCTACACGCTCCAGGCAACGGCGGGTCGCAAACACCCGCGCCAATACCCTCCAAT‑CCGGATCTACACCCAACCCGGCCAAAGCGCGAAAGCGGGGCAGCACCGGGAAATTCAAAGCTCCTGCTTTTGCTTCGCCCATCATCCTTAcc                                         >  2:1260278/1‑150 (MQ=255)                                                                                                                   ccTACACGCTCCAGGCAACGGCGGGTCGCATACTCCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAg                                     <  2:1323662/151‑1 (MQ=255)                                                                                                                   ccTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAg                                     <  2:854503/151‑1 (MQ=255)                                                                                                                      acacGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCCGCACGGGGAAATTCAACGCCCCTGCTTTTGCTTCGCTCATCATCATTACCGCAGcc                                   >  1:1350759/1‑149 (MQ=255)                                                                                                                      acacGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCt                                   <  1:1132488/150‑1 (MQ=255)                                                                                                                      acacGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCt                                   <  2:451157/150‑1 (MQ=255)                                                                                                                       cacGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTgc                                 <  2:1129993/151‑1 (MQ=255)                                                                                                                          gCTCCAGGCAACGGCGGGTCGCATACAACGGCG‑CAAATCCCATCAATGCCGGATCTACACCAACGCTGGCAAAAGCGCGAAAGCGGGCCAGCACGGGGAAAATCAACGCCCCGGCTTTTGCTTCGCTCCTCATCATTACCGCAGCTGCGcc                              >  2:892212/1‑151 (MQ=255)                                                                                                                           cTCCAGGCAACGGCGGGTCGCATACACCCGCG‑CCATTCCCCTCAATGCCGGAGCTACACCCACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAATTCAACGCTCGTGCTTTTGCTTCGCTCATCATCATTTCCCCAGCTGCGcc                              >  2:541962/1‑150 (MQ=255)                                                                                                                           cTCCAGGCAACGGCGGGGCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGGAAATTCAACGCTCGTGCTTTTGCTTCGCTCATCATCATTACCGCAGCTGCGcc                              >  1:480002/1‑150 (MQ=255)                                                                                                                                    aCGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTTGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTAt                     <  1:697088/150‑1 (MQ=255)                                                                                                                                      ggcggGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACGCCAACGCTGGCAAATGCGCGAATGCGGGCCAACACGGGTAAATTCAACGCTCGTGCTTTCGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATgg                   >  2:491652/1‑150 (MQ=255)                                                                                                                                      ggcggGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGga                  <  1:10363/151‑1 (MQ=255)                                                                                                                                       gcggGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgat                 >  2:1279244/1‑151 (MQ=255)                                                                                                                                          ggTCGCCTACACGGGCG‑CGATTCCCATCAAAGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGgatgat              <  2:1008853/151‑1 (MQ=255)                                                                                                                                             cGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGc            >  2:774780/1‑150 (MQ=255)                                                                                                                                              gCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAt          >  1:1336513/1‑151 (MQ=255)                                                                                                                                                aTACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCAtt         >  1:549898/1‑150 (MQ=255)                                                                                                                                                      cGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt  <  1:1351398/151‑1 (MQ=255)                                                                                                                                                      cGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCt  <  1:542647/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                       CCACCATTGACATTGACCCGACGCTCATCCCACTCAAGCATCTTGCCAACCGAAAGCGCCTGTGCAGCAAACGCTTCATTAGCCTCGATAAGATCGACTTCAGCCAACTGCCAGCCTACACGCTCCAGGCAACGGCGGGTCGCATACACCGGCG‑CAATTCCCATCAATGCCGGATCTACACCAACGCTGGCAAATGCGCGAATGCGGGCCAGCACGGGTAAATTCAACGCTCGTGCTTTGGCTTCGCTCATCATCATTACCGCAGCTGCGCCATCGTTTATGGATGATGCATTACCCGCT  >  NC_000913/2984568‑2984867

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