Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,967,132	A→G	H239R (CAC→CGC)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,967,132	0	A	G	100.0%	302.7 / NA	82	H239R (CAC→CGC)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (39/43);  total (39/43)

ATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCG  >  NC_000913/3966987‑3967275                                                                                                                                                  |                                                                                                                                                aTTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTc                                                                                                                                            <  1:736576/151‑1 (MQ=255)  ttGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTc                                                                                                                                            >  2:994916/1‑150 (MQ=255)   tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTc                                                                                                                                            >  2:471290/1‑149 (MQ=255)    gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAgg                                                                                                                                         <  2:202373/151‑1 (MQ=255)    gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAg                                                                                                                                          <  1:730207/150‑1 (MQ=255)     cTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGt                                                                                                                                        >  2:815862/1‑151 (MQ=255)       cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGtt                                                                                                                                       <  1:1103028/150‑1 (MQ=255)       cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGt                                                                                                                                        <  1:1131180/149‑1 (MQ=255)         gagCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCg                                                                                                                                    >  1:914592/1‑151 (MQ=255)            cATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGaa                                                                                                                                 <  2:73733/151‑1 (MQ=255)            cATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGa                                                                                                                                  >  1:497140/1‑150 (MQ=255)               tGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAAt                                                                                                                               <  2:714673/150‑1 (MQ=255)                gCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATg                                                                                                                              <  2:1206240/150‑1 (MQ=255)                 cTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGt                                                                                                                            >  2:987248/1‑151 (MQ=255)                   tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGt                                                                                                                            <  1:436195/149‑1 (MQ=255)                   tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTg                                                                                                                           >  2:1223487/1‑150 (MQ=255)                    aCAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGAt                                                                                                                         <  2:185194/151‑1 (MQ=255)                         caccCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCgag                                                                                                                     <  2:1189167/150‑1 (MQ=255)                           cccGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAg                                                                                                                  <  2:506751/151‑1 (MQ=255)                                aTTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaa                                                                                                             <  2:544460/151‑1 (MQ=255)                                 ttGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCaa                                                                                                             >  2:588827/1‑150 (MQ=255)                                      tgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt                                                                                                       >  1:324021/1‑151 (MQ=255)                                       gCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt                                                                                                       <  1:987802/150‑1 (MQ=255)                                       gCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCt                                                                                                       <  2:396818/150‑1 (MQ=255)                                         tGATGGTTCTGCTGATCGCCGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTg                                                                                                      >  2:1022560/1‑149 (MQ=255)                                          gATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTgg                                                                                                     <  1:987109/149‑1 (MQ=255)                                           aTGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt                                                                                                   >  2:965941/1‑150 (MQ=255)                                            tGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTg                                                                                                  >  2:780360/1‑150 (MQ=255)                                              gTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGc                                                                                               >  2:1136447/1‑151 (MQ=255)                                                 ctgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcac                                                                                             >  1:842099/1‑150 (MQ=255)                                                  tgctgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGcaca                                                                                            >  1:146110/1‑150 (MQ=255)                                                     tgATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                                                                                         >  2:1159459/1‑150 (MQ=255)                                                      gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                        <  1:554756/150‑1 (MQ=255)                                                      gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                        <  2:969585/150‑1 (MQ=255)                                                       aTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                        >  2:433323/1‑149 (MQ=255)                                                         cgacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaaga                                                                                     >  2:43015/1‑150 (MQ=255)                                                          gacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACg                                                                                   >  2:446841/1‑151 (MQ=255)                                                              aaCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTAt                                                                               <  1:1032536/151‑1 (MQ=255)                                                                 gTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCAt                                                                            >  2:756140/1‑151 (MQ=255)                                                                  tCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCAtt                                                                           >  1:621285/1‑151 (MQ=255)                                                                    cGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCt                                                                         <  2:160130/151‑1 (MQ=255)                                                                     ggAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTg                                                                        <  1:471305/151‑1 (MQ=255)                                                                       aagaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCt                                                                      <  1:1005108/151‑1 (MQ=255)                                                                       aagaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCt                                                                      <  2:624943/151‑1 (MQ=255)                                                                        ataattAACCGAGATGTACCGTCTAGTAAAGGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCCTCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCGAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCt                                                                      <  2:874749/145‑1 (MQ=255)                                                                        agaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTc                                                                     >  1:850085/1‑151 (MQ=255)                                                                         gaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCg                                                                    >  2:1202660/1‑151 (MQ=255)                                                                         gaagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCg                                                                    >  1:425778/1‑151 (MQ=255)                                                                            gTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACt                                                                 <  1:1068259/151‑1 (MQ=255)                                                                            gTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACt                                                                 <  2:343366/151‑1 (MQ=255)                                                                             tAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTc                                                                >  1:289661/1‑151 (MQ=255)                                                                              aaCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTcc                                                               <  1:794714/151‑1 (MQ=255)                                                                               aCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTcc                                                               >  1:580919/1‑150 (MQ=255)                                                                                   agaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCAc                                                          <  1:16560/151‑1 (MQ=255)                                                                                   agaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCAc                                                          >  2:44348/1‑151 (MQ=255)                                                                                     aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTc                                                        <  2:868033/151‑1 (MQ=255)                                                                                         aGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCt                                                    >  1:14561/1‑151 (MQ=255)                                                                                                 gTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcgcg                                            >  2:541112/1‑151 (MQ=255)                                                                                                    aaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCtt                                         <  1:257627/151‑1 (MQ=255)                                                                                                     aaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCtt                                         >  1:307893/1‑150 (MQ=255)                                                                                                      aGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGTGCGCGCTTac                                       <  2:574939/151‑1 (MQ=255)                                                                                                           aaGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTacaaca                                   <  1:492962/150‑1 (MQ=255)                                                                                                            aGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACAcc                                 >  1:880393/1‑151 (MQ=255)                                                                                                            aGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACAcc                                 >  2:348407/1‑151 (MQ=255)                                                                                                             gttgttGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCg                                <  1:220458/151‑1 (MQ=255)                                                                                                              ttgttgCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgt                               >  1:742267/1‑151 (MQ=255)                                                                                                                gttgCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgtt                              <  1:591778/150‑1 (MQ=255)                                                                                                                gttgCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgtt                              <  1:1204358/150‑1 (MQ=255)                                                                                                                gttgCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgtt                              <  1:825072/150‑1 (MQ=255)                                                                                                                      tCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCg                        <  2:733698/150‑1 (MQ=255)                                                                                                                      tCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCg                        <  1:124821/150‑1 (MQ=255)                                                                                                                         aCCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGt                    >  1:501157/1‑151 (MQ=255)                                                                                                                          ccTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCg                     <  2:621265/149‑1 (MQ=255)                                                                                                                          ccTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGt                    <  1:409832/150‑1 (MQ=255)                                                                                                                          ccTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTc                   >  1:1109340/1‑151 (MQ=255)                                                                                                                           cTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGGGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCGGGCGCGCGCTTACAGCACCGTTGTTCCGGc                      <  2:119315/147‑1 (MQ=255)                                                                                                                            tttGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTc                   >  1:256070/1‑149 (MQ=255)                                                                                                                              tGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAgg                >  2:91030/1‑150 (MQ=255)                                                                                                                              tGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGt               >  2:452427/1‑151 (MQ=255)                                                                                                                               gACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGt               <  1:541135/150‑1 (MQ=255)                                                                                                                                 cGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTaaa            <  2:108856/151‑1 (MQ=255)                                                                                                                                           tctcGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCg  >  2:234532/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                                ATTGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCG  >  NC_000913/3966987‑3967275

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