Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,186,152	T→G	L268R (CTG→CGG)	rpoC →	RNA polymerase, beta prime subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,186,152	0	T	G	100.0%	235.9 / NA	61	L268R (CTG→CGG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (33/28);  total (33/28)

AGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATC  >  NC_000913/4186005‑4186297                                                                                                                                                    |                                                                                                                                                  aGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTa                                                                                                                                                <  1:562541/151‑1 (MQ=255)   cGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTAtc                                                                                                                                              <  1:284089/151‑1 (MQ=255)      aTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTAtcgtc                                                                                                                                           >  2:404414/1‑151 (MQ=255)      aTCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTAtcgtc                                                                                                                                           <  2:811789/151‑1 (MQ=255)         aaaCTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTcgcg                                                                                                                                        >  1:695459/1‑151 (MQ=255)          aaCTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGt                                                                                                                                       >  1:473041/1‑151 (MQ=255)           aCTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTc                                                                                                                                      >  1:680835/1‑151 (MQ=255)              gctgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCAtt                                                                                                                                   >  1:1082084/1‑151 (MQ=255)               ctgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCCtt                                                                                                                                   >  1:60243/1‑150 (MQ=255)                tgGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTa                                                                                                                                  >  2:1010631/1‑150 (MQ=255)                     gCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTa                                                                                                                            <  2:674694/151‑1 (MQ=255)                       gttcgttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTaac                                                                                                                          >  1:330169/1‑151 (MQ=255)                        ttcgttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACCGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTaac                                                                                                                          >  1:1037114/1‑150 (MQ=255)                          cgttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTaacaa                                                                                                                        <  2:395413/150‑1 (MQ=255)                           gttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAAcc                                                                                                                      <  1:518369/151‑1 (MQ=255)                            ttcAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCg                                                                                                                     >  1:17996/1‑151 (MQ=255)                               aGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCt                                                                                                                  >  1:642437/1‑151 (MQ=255)                                 tCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTg                                                                                                                 >  1:5313/1‑150 (MQ=255)                                 tCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTg                                                                                                                 >  2:631544/1‑150 (MQ=255)                                  cTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGa                                                                                                                >  1:737921/1‑150 (MQ=255)                                   tGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGaa                                                                                                               >  2:173612/1‑150 (MQ=255)                                    ggTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGaaa                                                                                                              <  2:751267/150‑1 (MQ=255)                                    ggTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGaa                                                                                                               <  1:48955/149‑1 (MQ=255)                                       aaCAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTc                                                                                                          >  2:371438/1‑151 (MQ=255)                                         cAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTACGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctg                                                                                                        >  1:209032/1‑151 (MQ=255)                                          aaaCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctgc                                                                                                       >  2:396839/1‑151 (MQ=255)                                          aaaCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctgc                                                                                                       >  2:639137/1‑151 (MQ=255)                                           aaCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTctgct                                                                                                      <  1:573561/151‑1 (MQ=255)                                                    ggATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTgg                                                                                              <  2:935775/150‑1 (MQ=255)                                                       tgatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTg                                                                                           >  2:698348/1‑150 (MQ=255)                                                        gatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTgc                                                                                          <  2:338401/150‑1 (MQ=255)                                                              gACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGa                                                                                    <  2:414624/150‑1 (MQ=255)                                                                ccGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGAcat                                                                                 <  1:403896/151‑1 (MQ=255)                                                                 cGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGAcatc                                                                                <  1:410180/151‑1 (MQ=255)                                                                  gTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGAcatca                                                                               <  1:338547/151‑1 (MQ=255)                                                                    tCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGAcatcat                                                                              <  2:197957/150‑1 (MQ=255)                                                                       gCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcg                                                                            <  2:98507/149‑1 (MQ=255)                                                                        ccGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcgtac                                                                         >  2:601202/1‑151 (MQ=255)                                                                                gccgccAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACg                                                                   <  1:375713/149‑1 (MQ=255)                                                                                  cgccAGGTCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGaaa                                                                >  2:35204/1‑150 (MQ=255)                                                                                     cAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACg                                                            >  1:912715/1‑151 (MQ=255)                                                                                     cAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACg                                                            >  1:997419/1‑151 (MQ=255)                                                                                      aGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGt                                                           <  2:1115017/151‑1 (MQ=255)                                                                                            gCGTCCGCTGGTTCCGCTGGAGGGTTGCCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCCGCGCCGGACATCATCGTCCGTAACGAAAAACGAAtgct                                                      >  1:946657/1‑150 (MQ=255)                                                                                                ccGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAg                                                  <  2:501951/150‑1 (MQ=255)                                                                                                  gCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGa                                                >  1:1134839/1‑150 (MQ=255)                                                                                                   cTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGa                                                >  1:264131/1‑149 (MQ=255)                                                                                                   cTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGa                                                >  2:494899/1‑149 (MQ=255)                                                                                                     ggTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCg                                            <  2:597295/151‑1 (MQ=255)                                                                                                      gTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCgg                                           <  2:1035651/151‑1 (MQ=255)                                                                                                                       tCGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGAt                          >  2:824469/1‑151 (MQ=255)                                                                                                                        cGTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGAt                          <  2:626905/150‑1 (MQ=255)                                                                                                                         gTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAAc                                                             >  1:123411/1‑114 (MQ=255)                                                                                                                         gTTTCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAAc                                                             <  2:123432/114‑1 (MQ=255)                                                                                                                          tttCGCGACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAAc                       >  2:1061521/1‑151 (MQ=255)                                                                                                                             cgcgACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACgg                     <  2:954996/150‑1 (MQ=255)                                                                                                                             cgcgACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgt                    <  2:1098774/151‑1 (MQ=255)                                                                                                                               cgACTTCTGACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtcg                  >  2:549641/1‑151 (MQ=255)                                                                                                                                       gACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCgg            <  2:399654/149‑1 (MQ=255)                                                                                                                                       gACCTGAACGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTc          >  2:913594/1‑151 (MQ=255)                                                                                                                                               cGATCGGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATc  <  1:511249/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  AGCGTATCAAACTGCTGGAAGCGTTCGTTCAGTCTGGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATC  >  NC_000913/4186005‑4186297

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