Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,186,186	G→C	L279L (CTG→CTC)	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,186	0	G	C	100.0%	397.5 / NA	100	L279L (CTG→CTC)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (52/48);  total (52/48)

GGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTT  >  NC_000913/4186040‑4186329                                                                                                                                                   |                                                                                                                                                ggTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCaa                                                                                                                                               <  1:1560876/149‑1 (MQ=255)    aaCAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTc                                                                                                                                          >  1:913005/1‑151 (MQ=255)    aaCAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTc                                                                                                                                          >  2:1527111/1‑151 (MQ=255)      cAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTctg                                                                                                                                        >  2:1530140/1‑151 (MQ=255)      cAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTctg                                                                                                                                        <  1:690695/151‑1 (MQ=255)       aaaCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTctgc                                                                                                                                       >  2:426848/1‑151 (MQ=255)         aCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTctgctg                                                                                                                                     <  1:1321380/151‑1 (MQ=255)         aCCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTctgctg                                                                                                                                     <  2:882727/151‑1 (MQ=255)          ccAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTgg                                                                                                                                    <  2:273238/151‑1 (MQ=255)             gagTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATc                                                                                                                                 >  1:314605/1‑151 (MQ=255)                 ggATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGc                                                                                                                             <  2:811770/151‑1 (MQ=255)                 ggATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGc                                                                                                                             <  2:660348/151‑1 (MQ=255)                 ggATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGc                                                                                                                             >  1:37098/1‑151 (MQ=255)                   atgatCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTg                                                                                                                           <  1:1327151/151‑1 (MQ=255)                        ccTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCg                                                                                                                      >  2:959346/1‑151 (MQ=255)                           gACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGa                                                                                                                    <  2:199988/150‑1 (MQ=255)                           gACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGAc                                                                                                                   <  2:1578936/151‑1 (MQ=255)                           gACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGAc                                                                                                                   >  1:271432/1‑151 (MQ=255)                            aCCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGAca                                                                                                                  <  1:1166299/151‑1 (MQ=255)                              cGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGAcatc                                                                                                                >  2:985861/1‑151 (MQ=255)                               gTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGAcatca                                                                                                               <  2:1441749/151‑1 (MQ=255)                               gTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGAcatc                                                                                                                >  2:12579/1‑150 (MQ=255)                               gTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGCTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGCACGATCTGTATCGACGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGAcctc                                                                                                                >  1:1186132/1‑150 (MQ=255)                                  cTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcg                                                                                                            >  2:625520/1‑151 (MQ=255)                                   tGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcgt                                                                                                           >  2:76039/1‑151 (MQ=255)                                      cGGTACTGCCGCCAGATCTGCGTCCGCTGGTTACGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcgtacg                                                                                                        >  1:982974/1‑151 (MQ=255)                                       ggTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcgtacgt                                                                                                       <  2:1487737/151‑1 (MQ=255)                                        gttacgccgccAAATCTTCGTCCGCCGGTTCCGCTGGATGGTGGTCGTTTCGCCACTTCTGACCTTAACGATCTGTATTGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATcgtacgta                                                                                                      <  1:1486043/146‑1 (MQ=255)                                          aCTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTaa                                                                                                     <  1:183475/150‑1 (MQ=255)                                          aCTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTaa                                                                                                     <  2:871644/150‑1 (MQ=255)                                          aCTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAAc                                                                                                    >  2:1317477/1‑151 (MQ=255)                                           cTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAAc                                                                                                    <  2:207573/150‑1 (MQ=255)                                            tGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACg                                                                                                   >  1:523248/1‑150 (MQ=255)                                               cgccAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGaaa                                                                                                >  1:1486782/1‑150 (MQ=255)                                                gccAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGaaaa                                                                                               >  2:1452209/1‑150 (MQ=255)                                                  cAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACg                                                                                            >  1:986060/1‑151 (MQ=255)                                                  cAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACg                                                                                            <  2:299004/151‑1 (MQ=255)                                                    gATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGt                                                                                           >  1:1357507/1‑150 (MQ=255)                                                      tCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtg                                                                                        >  1:952200/1‑151 (MQ=255)                                                        tGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtgct                                                                                      >  2:984518/1‑151 (MQ=255)                                                         gCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTa                                                                                          <  2:1120339/146‑1 (MQ=255)                                                         gCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtgctg                                                                                     >  2:1434798/1‑151 (MQ=255)                                                         gCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtgc                                                                                       <  1:506019/149‑1 (MQ=255)                                                           gTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCa                                                                                   <  1:906740/151‑1 (MQ=255)                                                            tCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAtgctgc                                                                                    >  1:847048/1‑149 (MQ=255)                                                             ccGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAgg                                                                                 <  1:1311884/151‑1 (MQ=255)                                                              cGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAgg                                                                                 >  1:1106633/1‑150 (MQ=255)                                                              cGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGa                                                                                <  2:398754/151‑1 (MQ=255)                                                                 tGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGc                                                                             >  2:1516444/1‑151 (MQ=255)                                                                  ggTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGGCATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCg                                                                            <  1:779377/151‑1 (MQ=255)                                                                  ggTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCg                                                                            >  1:552269/1‑151 (MQ=255)                                                                  ggTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCg                                                                            <  1:314833/151‑1 (MQ=255)                                                                    ttCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCgg                                                                           >  1:1530674/1‑150 (MQ=255)                                                                      ccGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGt                                                                          >  2:485577/1‑149 (MQ=255)                                                                       cGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACATCCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGa                                                                       <  1:1052538/151‑1 (MQ=255)                                                                        gCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGAc                                                                      <  2:215293/151‑1 (MQ=255)                                                                          tGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGc                                                                    >  2:1405627/1‑151 (MQ=255)                                                                           ggATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGc                                                                    <  1:985549/150‑1 (MQ=255)                                                                             aTGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCct                                                                 <  2:482090/151‑1 (MQ=255)                                                                               ggtggtCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCctg                                                                <  2:1060821/150‑1 (MQ=255)                                                                                gtggtCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCctgct                                                              >  2:1325950/1‑151 (MQ=255)                                                                                  ggtCGTTTCGCGACTTCTGGCCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCctgctg                                                             <  1:439853/150‑1 (MQ=255)                                                                                     cGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGAt                                                          >  1:795009/1‑150 (MQ=255)                                                                                      gTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGAt                                                          >  2:1265793/1‑149 (MQ=255)                                                                                      gTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATa                                                         >  2:908058/1‑150 (MQ=255)                                                                                        ttCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAAc                                                       <  2:498520/150‑1 (MQ=255)                                                                                          cgcgACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgt                                                    <  2:1315362/151‑1 (MQ=255)                                                                                           gcgACTTCTGACCTGAACGATCTGTATCGTAGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtc                                                   >  1:763500/1‑151 (MQ=255)                                                                                            cgACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtc                                                   <  1:1500471/150‑1 (MQ=255)                                                                                             gACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGGCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtcgt                                                 <  1:1429182/151‑1 (MQ=255)                                                                                              aCTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtcgtc                                                <  1:1080462/151‑1 (MQ=255)                                                                                               cTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGgtcgtc                                                <  2:1004602/150‑1 (MQ=255)                                                                                                 tCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTcgc                                              <  2:566170/150‑1 (MQ=255)                                                                                                   tGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGt                                           >  2:837665/1‑151 (MQ=255)                                                                                                    gACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTcgcg                                             <  1:381679/148‑1 (MQ=255)                                                                                                    gACCTGAACGAGCTGTATCGTCGCGTCATTAACCGTAACACCCGTCTCAAACGTCGGCTGGGTCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTAGGCTGCAGGAAGGGGTAGACGCCCGGCTGGATAACGGTCGTCGCGGTc                                          <  1:981700/151‑1 (MQ=255)                                                                                                     aCCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCg                                         >  1:1030318/1‑151 (MQ=255)                                                                                                     aCCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCg                                         <  2:122562/151‑1 (MQ=255)                                                                                                     aCCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCg                                         >  2:1160909/1‑151 (MQ=255)                                                                                                     aCCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCg                                         >  1:1166360/1‑151 (MQ=255)                                                                                                         gAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCg                                     >  2:850440/1‑151 (MQ=255)                                                                                                             gATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATc                                  >  2:1121514/1‑150 (MQ=255)                                                                                                               tCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCAcc                               >  1:1549802/1‑151 (MQ=255)                                                                                                                cTGTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAACCGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCAcc                               >  1:938101/1‑150 (MQ=255)                                                                                                                  gTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATc                                  <  1:1122097/145‑1 (MQ=255)                                                                                                                  gTATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGt                            >  2:1083842/1‑151 (MQ=255)                                                                                                                   tATCGTCGCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGtt                           >  1:1136660/1‑151 (MQ=255)                                                                                                                      cgtcgCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCt                         <  2:344507/150‑1 (MQ=255)                                                                                                                      cgtcgCGTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTa                        <  2:268563/151‑1 (MQ=255)                                                                                                                         cgcgTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAAc                      >  2:629021/1‑150 (MQ=255)                                                                                                                         cgcgTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACa                     <  2:593328/151‑1 (MQ=255)                                                                                                                          gcgTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACa                     >  1:1021083/1‑150 (MQ=255)                                                                                                                           cgTCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAg                   <  1:49964/151‑1 (MQ=255)                                                                                                                             tCATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGc                  >  1:300641/1‑150 (MQ=255)                                                                                                                              cATTAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGt                >  2:278822/1‑151 (MQ=255)                                                                                                                                ttAACCGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTcc              >  1:755919/1‑151 (MQ=255)                                                                                                                                    ccGTAACAACCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTg          >  2:1156916/1‑151 (MQ=255)                                                                                                                                        aacaacCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAAt      <  2:613708/151‑1 (MQ=255)                                                                                                                                           aacCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCtt   >  2:847756/1‑151 (MQ=255)                                                                                                                                            acCGTCTCAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCttt  <  1:656274/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                GGTAACAAACCAGAGTGGATGATCCTGACCGTTCTGCCGGTACTGCCGCCAGATCTGCGTCCGCTGGTTCCGCTGGATGGTGGTCGTTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTT  >  NC_000913/4186040‑4186329

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