Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,270,605	G→A	10.3%	intergenic (+60/‑121)	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,605	0	G	A	10.3%	192.1 / 10.2	68	intergenic (+60/‑121)	mepS/pdeN	peptidoglycan DD‑endopeptidase/peptidoglycan LD‑carboxypeptidase/putative c‑di‑GMP phosphodiesterase PdeN
Reads supporting (aligned to +/- strand):  ref base G (33/28);  new base A (5/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.68e-01

CCAGCAGTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGAC  >  NC_000913/2270457‑2270753                                                                                                                                                     |                                                                                                                                                     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:929254/1‑150 (MQ=255)              tattatTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTg                                                                                                                                        >  2:1021247/1‑150 (MQ=255)              tattatTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTg                                                                                                                                        <  2:314568/150‑1 (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                                                                                                                                                    >  2:1080519/1‑122 (MQ=255)                              aatgaaCCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGc                                                                                                                                                    <  1:1080519/122‑1 (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                                                                                                           >  1:744428/1‑150 (MQ=255)                                           ggAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTc                                                                                                           >  2:1309540/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                                <  1:1089676/132‑1 (MQ=255)                                                                                                                                     aaaaaaGCACTGCCTAAGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGaa                                >  2:1089676/4‑135 (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:247241/1‑150 (MQ=255)                                                                                                                                           gCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCa           >  2:1251051/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)                                                                                                                                                    aGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGAc  <  1:230497/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCAGCAGTGGTGTTATTATTTCCAGCATGAATGAACCGTACTGGAAGAAGCGTTACAACGAAGCACGCCGGGTTCTCAGCCGCAGCTAATAAACCGTTTGGATGCAATCCCTTGGCTATCCTGACGAGTTAACTGAAAGCACTGCTTAGGCAGTGCTTTTTTGTTTTCATTCATCAGAGAAAATGATGTTTCCGCGTCTTGATCCAGGCTATAGTCCGGTCATTGTTATCTTTTAAATGTTGTCGTAATTTCAGGAAATTAACGGAATCATGTTCATACGCGCTCCCAATTTTGGAC  >  NC_000913/2270457‑2270753

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: ‑