Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,185,708	T→C	L120P (CTG→CCG)	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,185,708	0	T	C	100.0%	279.8 / NA	77	L120P (CTG→CCG)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (36/41);  total (36/41)

CTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACG  >  NC_000913/4185560‑4185851                                                                                                                                                     |                                                                                                                                                cTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATcgccgc                                                                                                                                               >  1:998962/1‑151 (MQ=255)           aGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTAt                                                                                                                                     <  2:466194/151‑1 (MQ=255)               cAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGt                                                                                                                                 >  2:542786/1‑151 (MQ=255)               cAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGt                                                                                                                                 >  2:156309/1‑151 (MQ=255)                   cgcCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTctgc                                                                                                                             >  1:648766/1‑151 (MQ=255)                    gcCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTctgct                                                                                                                            <  2:1074385/151‑1 (MQ=255)                    gcCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTctgct                                                                                                                            >  1:509251/1‑151 (MQ=255)                    gcCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTctgct                                                                                                                            >  2:771275/1‑151 (MQ=255)                     ccTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGGCTAAAGTACGCCATGAGCGTATGGGCCACATCGAACTGGCTTGCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGTCGTCCCGTATCGGTctgctg                                                                                                                           <  1:662561/151‑1 (MQ=255)                      cTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCtgctgc                                                                                                                          >  2:638884/1‑151 (MQ=255)                        gAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTgctgct                                                                                                                         <  1:496907/150‑1 (MQ=255)                         aaaCACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTc                                                                                                                        <  1:587158/150‑1 (MQ=255)                          aaCACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCg                                                                                                                       >  1:808179/1‑150 (MQ=255)                          aaCACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGa                                                                                                                      <  2:113713/151‑1 (MQ=255)                              ccGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCGCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATg                                                                                                                  <  1:165231/151‑1 (MQ=255)                                gTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGcc                                                                                                                >  2:355235/1‑151 (MQ=255)                                gTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGcc                                                                                                                <  1:661511/151‑1 (MQ=255)                                  ggCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCg                                                                                                               <  1:380156/150‑1 (MQ=255)                                     gTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTgc                                                                                                           >  2:1114331/1‑151 (MQ=255)                                     gTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTgc                                                                                                           >  2:99807/1‑151 (MQ=255)                                       cATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTgcgc                                                                                                         >  2:981097/1‑151 (MQ=255)                                       cATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTgcgc                                                                                                         <  2:497558/151‑1 (MQ=255)                                        aTCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGcgcg                                                                                                        >  1:430549/1‑151 (MQ=255)                                            gtgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGcgcg                                                                                                        <  2:460187/147‑1 (MQ=255)                                            gtgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGata                                                                                                     <  1:1101556/150‑1 (MQ=255)                                             tgAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGCCTGCTGCTCGATATGCCGCTGCGCGATATc                                                                                                   >  2:660362/1‑151 (MQ=255)                                               agaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCg                                                                                                  >  2:1044305/1‑150 (MQ=255)                                                gaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGaa                                                                                                <  2:823953/151‑1 (MQ=255)                                                 aaGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGaa                                                                                                <  2:1011373/150‑1 (MQ=255)                                                 aaGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAAc                                                                                               >  2:143651/1‑151 (MQ=255)                                                  aGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAAcg                                                                                              <  1:223658/151‑1 (MQ=255)                                                  aGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAAcg                                                                                              <  2:1051799/151‑1 (MQ=255)                                                   gTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAAcgc                                                                                             <  2:1217026/151‑1 (MQ=255)                                                        gcgTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgtact                                                                                        >  2:648344/1‑151 (MQ=255)                                                             gAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCgtactgtac                                                                                    <  1:252280/150‑1 (MQ=255)                                                              aaGTGACCCAGACTAAAGTACGCAGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACtt                                                                                  >  1:92802/1‑151 (MQ=255)                                                               aGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACttt                                                                                 <  1:1158747/151‑1 (MQ=255)                                                                 tGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTg                                                                                >  2:544643/1‑150 (MQ=255)                                                                             aaGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGt                                                                   <  1:487735/151‑1 (MQ=255)                                                                                tACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTAt                                                                >  1:64847/1‑151 (MQ=255)                                                                                tACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTAt                                                                >  1:998231/1‑151 (MQ=255)                                                                                  cGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCg                                                              <  1:337360/151‑1 (MQ=255)                                                                                   gCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCg                                                              >  1:553420/1‑150 (MQ=255)                                                                                   gCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCg                                                              >  2:1034863/1‑150 (MQ=255)                                                                                   gCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGa                                                             <  1:390853/151‑1 (MQ=255)                                                                                      gTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAgg                                                          <  2:922278/151‑1 (MQ=255)                                                                                      gTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAgg                                                          <  1:1041645/151‑1 (MQ=255)                                                                                      gTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAg                                                           <  1:1143246/150‑1 (MQ=255)                                                                                       tGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGc                                                         >  2:645598/1‑151 (MQ=255)                                                                                           cGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGt                                                      >  1:258166/1‑150 (MQ=255)                                                                                                gggCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGATGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGAcc                                                >  1:852283/1‑151 (MQ=255)                                                                                                    cacaTCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAAcc                                            >  2:103742/1‑151 (MQ=255)                                                                                                     acaTCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCt                                           >  2:435648/1‑151 (MQ=255)                                                                                                      caTCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTg                                          >  2:722686/1‑151 (MQ=255)                                                                                                      caTCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTg                                          <  2:70220/151‑1 (MQ=255)                                                                                                         cGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGaa                                       <  2:900841/151‑1 (MQ=255)                                                                                                               ggCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTcag                                 <  2:500730/151‑1 (MQ=255)                                                                                                                gCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTcagc                                <  2:345773/151‑1 (MQ=255)                                                                                                                     cccGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATc                           <  2:263078/151‑1 (MQ=255)                                                                                                                       cGACTGCGCACATCTGGTTCCTGAAATCGCCGCCGTCCTGTATAGGTCTCCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGAt                            >  2:738488/1‑148 (MQ=255)                                                                                                                            gcgcACATCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCctgactg                    <  2:808212/151‑1 (MQ=255)                                                                                                                               cacaTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCctgactga                   >  1:691995/1‑149 (MQ=255)                                                                                                                                acaTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCctgactga                   >  2:951768/1‑148 (MQ=255)                                                                                                                                 caTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAagag               <  2:116385/151‑1 (MQ=255)                                                                                                                                 caTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAaga                <  2:92483/150‑1 (MQ=255)                                                                                                                                  aTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGc              >  2:1141314/1‑151 (MQ=255)                                                                                                                                  aTCTGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGc              >  1:1201221/1‑151 (MQ=255)                                                                                                                                     tGGTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGt           >  1:29131/1‑151 (MQ=255)                                                                                                                                      ggTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCa             <  1:355224/148‑1 (MQ=255)                                                                                                                                       gTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTAt         <  2:648784/151‑1 (MQ=255)                                                                                                                                       gTTCCTGAAATCGCCGCCGTCCCGTATCGGTCTGCTACTCGATGTGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTAt         <  1:651603/151‑1 (MQ=255)                                                                                                                                          ccTGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCt       <  2:817142/150‑1 (MQ=255)                                                                                                                                            tGAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTgg     >  2:690817/1‑150 (MQ=255)                                                                                                                                             gAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTgg     <  2:50144/149‑1 (MQ=255)                                                                                                                                             gAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTgg     <  1:143645/149‑1 (MQ=255)                                                                                                                                             gAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGAc   <  1:1130247/151‑1 (MQ=255)                                                                                                                                             gAAATCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGAc   <  2:1173167/151‑1 (MQ=255)                                                                                                                                              aaaTCGCCGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGAcg  >  2:903072/1‑151 (MQ=255)                                                                                                                                                     |                                                                                                                                                CTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCGAACGCGTACTGTACTTTGAATCCTATGTGGTTATCGAAGGCGGTATGACCAACCTGGAACGTCAGCAGATCCTGACTGAAGAGCAGTATCTGGACG  >  NC_000913/4185560‑4185851

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