Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,187,619	C→A	T757K (ACG→AAG)	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,187,619	0	C	A	100.0%	250.9 / NA	62	T757K (ACG→AAG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (41/21);  total (41/21)

AACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGC  >  NC_000913/4187471‑4187762                                                                                                                                                     |                                                                                                                                                aaCCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGc                                                                                                                                               >  2:613190/1‑151 (MQ=255) aaCCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGc                                                                                                                                               <  1:80670/151‑1 (MQ=255)      tGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATc                                                                                                                                          >  1:31785/1‑151 (MQ=255)        aCGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCAc                                                                                                                                        >  1:931261/1‑151 (MQ=255)          ggTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACcg                                                                                                                                      >  2:480215/1‑151 (MQ=255)           gTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACcgc                                                                                                                                     >  2:385637/1‑151 (MQ=255)            tCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACcgcg                                                                                                                                    >  2:604924/1‑151 (MQ=255)              aGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGa                                                                                                                                   >  1:256221/1‑150 (MQ=255)              aGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGa                                                                                                                                   >  2:355135/1‑150 (MQ=255)              aGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGa                                                                                                                                   >  1:618161/1‑150 (MQ=255)               ggAAGAGAAGCAGGTTTCCTTCAACAGCATCTACACGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCTTCGAAAAGCCAATCACCGCGAAc                                                                                                                                 <  2:857196/151‑1 (MQ=255)                  agagAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTc                                                                                                                              >  2:540797/1‑151 (MQ=255)                  agagAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTc                                                                                                                              >  1:895418/1‑151 (MQ=255)                    agaAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCg                                                                                                                            >  1:841790/1‑151 (MQ=255)                      aaGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTg                                                                                                                          >  1:462334/1‑151 (MQ=255)                        gCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGa                                                                                                                         >  1:158991/1‑150 (MQ=255)                         cAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGa                                                                                                                         >  2:925368/1‑149 (MQ=255)                         cAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGa                                                                                                                         >  1:885718/1‑149 (MQ=255)                              ttCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTg                                                                                                                  >  1:406658/1‑151 (MQ=255)                                ccTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGaa                                                                                                                <  1:873698/151‑1 (MQ=255)                                  ttCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACg                                                                                                              >  1:721903/1‑151 (MQ=255)                                    cAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGt                                                                                                             >  1:75675/1‑150 (MQ=255)                                             cTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTAc                                                                                                   >  1:912747/1‑151 (MQ=255)                                               aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACtt                                                                                                 >  2:457758/1‑151 (MQ=255)                                               aCATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACt                                                                                                  <  1:816394/150‑1 (MQ=255)                                               aCATGATGGCCGACTCCGGTGCGCGTAGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTAc                                                                                                   <  1:832972/149‑1 (MQ=255)                                                 atgatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTc                                                                                                >  1:488312/1‑150 (MQ=255)                                                  tgatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAt                                                                                              >  1:868141/1‑151 (MQ=255)                                                   gatgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtc                                                                                             <  1:714794/151‑1 (MQ=255)                                                    atgGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtct                                                                                            >  1:431430/1‑151 (MQ=255)                                                      ggCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCAtctc                                                                                           <  1:486957/150‑1 (MQ=255)                                                      ggCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGACGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTcc                                                                                          >  1:456741/1‑151 (MQ=255)                                                          gACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccc                                                                                      >  1:666579/1‑151 (MQ=255)                                                          gACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccc                                                                                      >  1:533593/1‑151 (MQ=255)                                                           aCTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTccaccc                                                                                      >  1:477121/1‑150 (MQ=255)                                                                ggTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTg                                                                                >  2:566040/1‑151 (MQ=255)                                                                    cgcgTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCg                                                                            <  1:868854/151‑1 (MQ=255)                                                                               gCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCGGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCt                                                                   <  2:725230/149‑1 (MQ=255)                                                                                       gATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATAcc                                                         <  1:223329/151‑1 (MQ=255)                                                                                        aTTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCg                                                        <  1:457838/151‑1 (MQ=255)                                                                                         ttCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGc                                                       >  1:65821/1‑151 (MQ=255)                                                                                         ttCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGc                                                       >  1:502521/1‑151 (MQ=255)                                                                                                 cTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGaaa                                                >  2:345176/1‑150 (MQ=255)                                                                                                   tGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACt                                             >  1:678819/1‑151 (MQ=255)                                                                                                     cTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGc                                           >  1:252645/1‑151 (MQ=255)                                                                                                          aTGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGACCTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAAc                                       <  1:916502/150‑1 (MQ=255)                                                                                                            gCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTc                                     <  1:139642/150‑1 (MQ=255)                                                                                                            gCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGCGAAGGTCTGAACGTACTCCCATCCTTCATCTCAAGCCACGGCGCTCGACAAGGTCTGGCGGATCCCGAACTGCAAACTGCGAAc                                       >  2:252469/1‑148 (MQ=255)                                                                                                                ggTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGt                                 <  2:182476/150‑1 (MQ=255)                                                                                                                      aTGGCGAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTg                           >  2:851113/1‑150 (MQ=255)                                                                                                                           gAAGCCGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACtcgt                     <  2:31784/151‑1 (MQ=255)                                                                                                                                cGGATGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCt                <  2:433057/151‑1 (MQ=255)                                                                                                                                    tGGCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGtt            >  2:595072/1‑151 (MQ=255)                                                                                                                                     ggCTCCATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGtt            <  2:273331/150‑1 (MQ=255)                                                                                                                                         ccATCATCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGt       >  2:645495/1‑151 (MQ=255)                                                                                                                                          catcatCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTg      <  2:256185/151‑1 (MQ=255)                                                                                                                                          catcatCGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTg      <  1:110455/151‑1 (MQ=255)                                                                                                                                            tcatcGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGgc    >  1:623851/1‑151 (MQ=255)                                                                                                                                             catcGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGgcg   <  1:194692/151‑1 (MQ=255)                                                                                                                                             catcGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGgc    <  2:938786/150‑1 (MQ=255)                                                                                                                                               tcGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGgcgc  >  1:630689/1‑150 (MQ=255)                                                                                                                                               tcGAAAAGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGgcgc  >  2:893715/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                AACCGTGACGGTCAGGAAGAGAAGCAGGTTTCCTTCAACAGCATCTACATGATGGCCGACTCCGGTGCGCGTGGTTCTGCGGCACAGATTCGTCAGCTTGCTGGTATGCGTGGTCTGATGGCGAAGCCGGATGGCTCCATCATCGAAACGCCAATCACCGCGAACTTCCGTGAAGGTCTGAACGTACTCCAGTACTTCATCTCCACCCACGGTGCTCGTAAAGGTCTGGCGGATACCGCACTGAAAACTGCGAACTCCGGTTACCTGACTCGTCGTCTGGTTGACGTGGCGC  >  NC_000913/4187471‑4187762

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