Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	4,235,408	Δ1 bp	intergenic (‑2175/‑499)	lysC ← / → yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	4,235,407	0	A	.	98.1%	224.9 / ‑0.6	55	intergenic (‑2174/‑500)	lysC/yjbE	lysine‑sensitive aspartokinase 3/extracellular polysaccharide production threonine‑rich protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (23/29);  total (24/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGA  >  NC_000913_3_bme_pgi/4235259‑4235548                                                                                                                                                     |                                                                                                                                              tCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                               <  1:654203/149‑1 (MQ=255) gcaagggcTTCGATCAATGGGGCGTGGAACTGGGTAAGGAAGTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGCGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGa                                                                                                                                               <  1:3477/142‑1 (MQ=255)  cAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGat                                                                                                                                              <  2:409266/149‑2 (MQ=255)    aTGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGATTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGatca                                                                                                                                            <  1:370457/149‑4 (MQ=255)       ccTTGGATCAAGGGGGCGGGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcg                                                                                                                                        <  2:307070/149‑1 (MQ=255)         ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                      >  2:178141/1‑149 (MQ=255)         ttCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCAtcgtc                                                                                                                                      >  1:713688/1‑149 (MQ=255)                 agggggCGGGGAACTGGTTAAAGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCCACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTa                                                                                                                              <  1:555685/147‑1 (MQ=255)                 agggggCGGGGAACTGGGGAAGGAACTGGCGACGGAACTTTTGCGGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTGGCCACGCTTCCACGCAAGGGTTTGTCGCGCATCTGTATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTa                                                                                                                              <  1:555460/147‑1 (MQ=255)                    gggCGTGGAACTGGGTAAGGGACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGGATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGc                                                                                                                           <  1:512498/149‑1 (MQ=255)                    gggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGc                                                                                                                           <  1:637284/149‑1 (MQ=255)                    gcgCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGc                                                                                                                           <  1:531433/147‑1 (MQ=255)                     ggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCGTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGcc                                                                                                                          <  1:703536/149‑1 (MQ=255)                     ggCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGGTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGcc                                                                                                                          <  1:700528/149‑1 (MQ=255)                      tcggggAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                         <  2:691039/145‑1 (MQ=255)                      gCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                         >  1:526643/1‑149 (MQ=255)                      gCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGGATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCg                                                                                                                         <  1:68270/149‑1 (MQ=255)                            aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGCAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGCATCAg                                                                                                                   >  1:105033/1‑149 (MQ=255)                            aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAg                                                                                                                   <  2:575235/149‑1 (MQ=255)                            aaCTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTACATACGCGAAGGAAAGCTTG‑ATCATCGTCGCTATGTAGGCCGGATacg                                                                                                                   >  2:407317/1‑147 (MQ=255)                               tGGGTAAGGAACTGGCGAAGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCg                                                                                                                <  1:278794/149‑1 (MQ=255)                                  gTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGGGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGC‑TGGATCATCGTCGATATGTAGGCCGGATAAGGCGTTc                                                                                                             >  2:94096/1‑149 (MQ=255)                                    aaGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAc                                                                                                           >  2:15251/1‑149 (MQ=255)                                      ggAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTGGGCCGGATAAGGCGTTCAcgc                                                                                                         >  1:63993/1‑149 (MQ=255)                                      ggAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgc                                                                                                         >  2:327076/1‑149 (MQ=255)                                       gAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                        <  2:456306/149‑1 (MQ=255)                                       gAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCCATATGTAGGCCGGATAAGGCGTTCAcgcc                                                                                                        >  2:272439/1‑149 (MQ=255)                                               gACGGAACTTTTGCCGATTGTATCGGGCAAGGAAGGGGCGAGTGGTCGAGACGCTTCCACGCAAGGGCTTGTCGCGCCTCTGCATGCGCGAAGGAAGGCTAG‑ATCAGCGTCGATATGTAGGACCGAAAAAGCCTTCACCCCGCCTTCgg                                                                                                >  2:194951/1‑149 (MQ=255)                                                aCGGAAGTTTTGCCGGTTGGATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                               <  2:364708/149‑1 (MQ=255)                                                aCGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                               >  1:328771/1‑149 (MQ=255)                                                  ggAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGATTACACGCAAGGGCTTGTCGCGCATCTCCATGCGCGCAGGAAAGATTG‑ATCCTCGTCGATATGTAAGACGGCTCAGGCGTTCGAGCAGCATCCGGccc                                                                                             >  1:333419/1‑147 (MQ=255)                                                   gAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAc                                                                                            >  1:302451/1‑149 (MQ=255)                                                       ttttGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGAt                                                                                        <  2:713688/149‑1 (MQ=255)                                                          tGCCGGTTGTACCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGcc                                                                                     <  1:702588/149‑1 (MQ=255)                                                             cGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTCGTCGCGCTTCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGa                                                                                  <  1:653911/149‑1 (MQ=255)                                                              ggTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGAt                                                                                 <  1:178141/149‑1 (MQ=255)                                                                   tttCGGGCAAGGAAGGCGCGAGTGGGCGGGACGCTTCCACGCAAGGGTTTGTCGCGCATCTGCAGGCGCGAAGGAAAGCTTG‑ATCATTGTCGATATGTGGGCCGGATAAGGCGTTCACGCCGCTTCCGGTAACCGATGCCTGATGCGAc                                                                            <  2:89707/147‑1 (MQ=255)                                                                     tCGGGCAGGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGAcgc                                                                          <  1:340577/149‑1 (MQ=255)                                                                           aaggaaggGGCGGGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTcgc                                                                    <  2:631449/149‑1 (MQ=255)                                                                               aaggGGCGAGTGGGCGCGACGGTTCCACGCAAGGGCTGGTCGCGTATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCt                                                                <  2:96431/149‑1 (MQ=37)                                                                                aggGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCtt                                                               >  1:167248/1‑149 (MQ=255)                                                                                            tCGCGACGCTCCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                   <  1:379111/149‑1 (MQ=255)                                                                                            gaccgACGCTTCCACGCCAGGGTTGTTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ACCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACa                                                   <  1:381015/146‑1 (MQ=25)                                                                                                    cTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATg                                           <  2:167248/149‑1 (MQ=255)                                                                                                      tCCACGCACGGGCTGGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCCTCGTCGATATGTAGGCCGGATAAGGCGTCCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGcc                                         <  2:368544/149‑1 (MQ=25)                                                                                                         aCGCAAGGGCTTGTCGCGCATCTGCATACCTGCAGCAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGat                                      >  1:726386/1‑149 (MQ=25)                                                                                                           gCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGatct                                    >  2:317406/1‑147 (MQ=255)                                                                                                            cAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGACTACAGGTCGATGCCGATATg                                   >  1:488045/1‑149 (MQ=255)                                                                                                                   ttGTCGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                            >  1:76766/1‑149 (MQ=255)                                                                                                                   ttGTAGCGCATCTGCATGCGCGAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCg                            >  1:212412/1‑149 (MQ=255)                                                                                                                          gcATTTTCATGCGCGAAGCAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCgg                     <  1:87878/149‑1 (MQ=255)                                                                                                                              cagCATGCCCGAAGGAAAGCTTGACT‑ATCGTCGATATGTAGGCCGGATAAGGCGTTAACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAAt                 <  2:237949/147‑1 (MQ=25)                                                                                                                                   tGCGTGAAGGCAAG‑TTTCATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATa            <  2:62159/149‑1 (MQ=25)                                                                                                                                      gggAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt         <  2:726386/147‑1 (MQ=255)                                                                                                                                      gcgAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACAt         >  1:658917/1‑149 (MQ=255)                                                                                                                                        gAAGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCCTCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGACTACAGGTCGATGCCGATATGTACTTCGTATTCGGAAATTAATACATAg       >  2:77648/1‑149 (MQ=37)                                                                                                                                          aGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCTTTCCCGCCGCATCCGGTAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCa     <  1:536046/149‑1 (MQ=255)                                                                                                                                          aGGAAAGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCa     >  2:438436/1‑149 (MQ=255)                                                                                                                                             aaaGCTTG‑ATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGTTGGCGATATGTACATCGTATTCGGCAATTAATACATAGCACga  >  1:136199/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                              TCAATGCCTTCGATCAATGGGGCGTGGAACTGGGTAAGGAACTGGCGACGGAACTTTTGCCGGTTGTATCGGGCAAGGAAGGGGCGAGTGGTCGCGACGCTTCCACGCAAGGGCTTGTCGCGCATCTGCATGCGCGAAGGAAAGCTTGAATCATCGTCGATATGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAACCGATGCCTGATGCGACGCGGTCGCGTCTTATCAGGCCTACAGGTCGATGCCGATATGTACATCGTATTCGGCAATTAATACATAGCACGA  >  NC_000913_3_bme_pgi/4235259‑4235548

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