Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,270,602	T→C	10.6%	intergenic (+57/‑124)	mepS → / → pdeN	peptidoglycan DD‑endopeptidase/peptidoglycan LD‑carboxypeptidase/putative c‑di‑GMP phosphodiesterase PdeN

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,270,602	0	T	C	10.6%	179.8 / 11.0	66	intergenic (+57/‑124)	mepS/pdeN	peptidoglycan DD‑endopeptidase/peptidoglycan LD‑carboxypeptidase/putative c‑di‑GMP phosphodiesterase PdeN
Reads supporting (aligned to +/- strand):  ref base T (33/26);  new base C (5/2);  total (38/28)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CCAGCAGTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTT  >  NC_000913/2270457‑2270748                                                                                                                                                  |                                                                                                                                                   ccaGCAGTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAgg                                                                                                                                                >  1:56980/1‑150 (MQ=255)  cagcagTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGc                                                                                                                                               <  1:816280/150‑1 (MQ=255)      agTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTg                                                                                                                                           <  1:1221045/150‑1 (MQ=255)            gtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCttttt                                                                                                                                     <  2:807986/150‑1 (MQ=255)            gtTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGAAGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttca                                                                                                                                     <  2:560832/150‑3 (MQ=255)             ttattaTTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCCTAAGCAGTGCtttcag                                                                                                                                    >  1:1429222/1‑147 (MQ=255)              tattatTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTg                                                                                                                                   >  2:1021247/1‑150 (MQ=255)              tattatTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTg                                                                                                                                   <  2:314568/150‑1 (MQ=255)              tattatTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTg                                                                                                                                   >  2:929254/1‑150 (MQ=255)              tattatTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTg                                                                                                                                   <  1:1021247/150‑1 (MQ=255)                 tattTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGttt                                                                                                                                <  1:296814/150‑1 (MQ=255)                   tttCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTttc                                                                                                                              <  1:567699/150‑1 (MQ=255)                       cAGCATGAATGAACCGTACTGGAAGAAGCGCTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTttcattc                                                                                                                          >  2:420589/1‑150 (MQ=255)                         gCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTtcattcat                                                                                                                        >  2:1008118/1‑150 (MQ=255)                            tgaatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCag                                                                                                                     >  1:937590/1‑150 (MQ=255)                             gaatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCaga                                                                                                                    <  2:145843/150‑1 (MQ=255)                             gaatgaaACGTACTGGAAGAAGCGTTACAAAGAAGCACGCCGGGTTCTCAGCCGCATCTAATAAAACGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCaga                                                                                                                    >  2:677777/1‑150 (MQ=255)                              aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGc                                                                                                                                               <  1:1080519/122‑1 (MQ=255)                              aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGc                                                                                                                                               >  2:1080519/1‑122 (MQ=255)                              aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCagag                                                                                                                   >  2:1386117/1‑150 (MQ=255)                                 gaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGaaa                                                                                                                >  1:268308/1‑150 (MQ=255)                                   aCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAat                                                                                                              <  2:937590/150‑1 (MQ=255)                                      gTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatgat                                                                                                           <  1:86388/150‑1 (MQ=255)                                      gTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAatgat                                                                                                           <  1:91507/150‑1 (MQ=255)                                         cTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGtt                                                                                                        >  2:347269/1‑150 (MQ=255)                                           ggAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTc                                                                                                      >  2:1309540/1‑150 (MQ=255)                                           ggAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTc                                                                                                      >  1:744428/1‑150 (MQ=255)                                              agaagCGCTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCcgc                                                                                                   <  1:420589/150‑1 (MQ=255)                                                 agCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTc                                                                                                <  2:205368/150‑1 (MQ=255)                                                      tACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGAt                                                                                           >  2:451518/1‑150 (MQ=255)                                                        cAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATcc                                                                                         <  1:805833/150‑1 (MQ=255)                                                          aCGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAg                                                                                       >  1:261806/1‑150 (MQ=255)                                                            gAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGAGTCTTGATACAGGc                                                                                     >  1:781529/1‑150 (MQ=255)                                                               gCACGCCGGGTGCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCtat                                                                                  <  2:168310/150‑1 (MQ=255)                                                                           ctcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCAtt                                                                      <  2:413354/150‑1 (MQ=255)                                                                            tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGTGTCTTGATCCAGGCTATAGTCCGGTCATTg                                                                     >  1:451215/1‑150 (MQ=255)                                                                            tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTg                                                                     <  2:56980/150‑1 (MQ=255)                                                                            tcAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTg                                                                     >  2:1378927/1‑150 (MQ=255)                                                                                      cTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTa                                                           >  2:319554/1‑150 (MQ=255)                                                                                       taataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTaa                                                          >  1:544405/1‑150 (MQ=255)                                                                                        aataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTaaa                                                         >  2:404381/1‑150 (MQ=255)                                                                                         ataaACCGTTTGGATGCAATCCCTTGGCTATCCTGACGACTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTAATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAt                                                        >  1:1230040/1‑150 (MQ=255)                                                                                            aaCCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAAtgtt                                                     >  1:697249/1‑150 (MQ=255)                                                                                               cGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTc                                                  <  2:235755/150‑1 (MQ=255)                                                                                                gTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCg                                                 <  1:628458/150‑1 (MQ=255)                                                                                                  ttGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTa                                               <  1:677777/150‑1 (MQ=255)                                                                                                       tGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTc                                          >  2:1223745/1‑150 (MQ=255)                                                                                                       tGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTc                                          <  2:650905/150‑1 (MQ=255)                                                                                                           aTCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGa                                      >  1:1019240/1‑150 (MQ=255)                                                                                                              ccTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAAt                                   >  2:539031/1‑150 (MQ=255)                                                                                                               ctaggCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAAtt                                  <  1:404028/147‑1 (MQ=255)                                                                                                                         cTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCa                        <  1:1008118/150‑1 (MQ=255)                                                                                                                         cTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCa                        >  2:164760/1‑150 (MQ=255)                                                                                                                              gAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTc                   <  1:929254/150‑1 (MQ=255)                                                                                                                                  aacaaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATAc               >  1:1136958/7‑150 (MQ=255)                                                                                                                                   acaaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATAcg              >  2:525123/6‑150 (MQ=255)                                                                                                                                   aaCTGAAAGCGCTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCAAAcg              >  1:300494/1‑150 (MQ=255)                                                                                                                                   aaCTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATAcg              >  2:230497/1‑150 (MQ=255)                                                                                                                                     aaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGaa                           >  2:1089676/4‑135 (MQ=255)                                                                                                                                     aaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGaa                           <  1:1089676/132‑1 (MQ=255)                                                                                                                                     aaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATAcgcg            >  2:1350275/4‑150 (MQ=255)                                                                                                                                       gAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCt          >  2:1222428/1‑150 (MQ=255)                                                                                                                                         aaGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTcc        >  2:183568/1‑150 (MQ=255)                                                                                                                                           gCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCa      >  2:1251051/1‑150 (MQ=255)                                                                                                                                           gCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCa      >  2:247241/1‑150 (MQ=255)                                                                                                                                               tGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAAttt  <  2:162893/150‑1 (MQ=255)                                                                                                                                                  ttAGGCAGTGCTTTTTTGTCTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAAt                          <  1:155971/123‑1 (MQ=255)                                                                                                                                                  ttAGGCAGTGCTTTTTTGTATTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAAt                          >  2:155971/1‑123 (MQ=255)                                                                                                                                                  |                                                                                                                                                   CCAGCAGTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTT  >  NC_000913/2270457‑2270748

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