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%	304.1 / NA	81	H239R (CAC→CGC)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (36/45);  total (36/45)

TGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGT  >  NC_000913/3966989‑3967277                                                                                                                                                |                                                                                                                                                  tGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCa                                                                                                                                             <  1:84826/150‑1 (MQ=255)  gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAg                                                                                                                                            <  1:421598/150‑1 (MQ=255)  gCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAg                                                                                                                                            <  1:99028/150‑1 (MQ=255)     cAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGtt                                                                                                                                         <  2:296090/150‑1 (MQ=255)        agCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCgg                                                                                                                                     >  2:718487/1‑151 (MQ=255)              gCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATgg                                                                                                                               <  1:526641/151‑1 (MQ=255)                 tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTCCCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGt                                                                                                                              <  1:426293/149‑1 (MQ=255)                 tACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGt                                                                                                                              <  2:443533/149‑1 (MQ=255)                  aCAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGa                                                                                                                            <  1:906476/150‑1 (MQ=255)                  aCAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGa                                                                                                                            <  1:455991/150‑1 (MQ=255)                    aaCCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCg                                                                                                                         >  2:520075/1‑151 (MQ=255)                     accaccCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCga                                                                                                                        <  1:653665/151‑1 (MQ=255)                         cccGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCgaga                                                                                                                      >  2:61728/1‑149 (MQ=255)                                 gtgtgCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAAc                                                                                                             <  2:129199/150‑1 (MQ=255)                                        gATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTgg                                                                                                       <  2:142079/149‑1 (MQ=255)                                        gATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGtt                                                                                                     <  2:1067240/151‑1 (MQ=255)                                           ggTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAg                                                                                                  <  2:494571/151‑1 (MQ=255)                                           ggTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAg                                                                                                  <  2:85972/151‑1 (MQ=255)                                            gTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAg                                                                                                  >  2:634159/1‑150 (MQ=255)                                            gTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGc                                                                                                 >  1:605542/1‑151 (MQ=255)                                             ttCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGc                                                                                                 >  1:1023047/1‑150 (MQ=255)                                                    gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                          >  2:830986/1‑150 (MQ=255)                                                    gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                                                                                           <  1:670780/149‑1 (MQ=255)                                                    gATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaaga                                                                                           <  1:856289/149‑1 (MQ=255)                                                     aTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaag                                                                                        <  2:792216/151‑1 (MQ=255)                                                     aTCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaa                                                                                          >  2:265721/1‑149 (MQ=255)                                                       cgacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACaagaaaga                                                                                       <  1:191444/150‑1 (MQ=255)                                                        gacgaACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACg                                                                                     >  2:1028064/1‑151 (MQ=255)                                                            aaCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTAt                                                                                 <  2:31605/151‑1 (MQ=255)                                                             aCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATc                                                                                <  2:894153/151‑1 (MQ=255)                                                             aCGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATc                                                                                <  1:68214/151‑1 (MQ=255)                                                              cGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCa                                                                               <  1:139916/151‑1 (MQ=255)                                                               gTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGCAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCAt                                                                              >  2:770914/1‑151 (MQ=255)                                                                 ccGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTc                                                                            >  1:1194399/1‑151 (MQ=255)                                                                  cGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCt                                                                           <  1:1034843/151‑1 (MQ=255)                                                                  cGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCt                                                                           <  2:59432/151‑1 (MQ=255)                                                                   ggAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTg                                                                          <  2:288391/151‑1 (MQ=255)                                                                        aagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCg                                                                      >  1:704116/1‑150 (MQ=255)                                                                        aagTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCg                                                                      >  2:348798/1‑150 (MQ=255)                                                                             aCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCa                                                                <  1:601437/151‑1 (MQ=255)                                                                             aCCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCa                                                                <  1:551148/151‑1 (MQ=255)                                                                                gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTc                                                                       >  2:335152/1‑141 (MQ=255)                                                                                gagaTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTc                                                                       <  1:335073/141‑1 (MQ=255)                                                                                   aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTc                                                          <  2:457280/151‑1 (MQ=255)                                                                                   aTGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTc                                                          <  2:640027/151‑1 (MQ=255)                                                                                    tGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCg                                                         >  1:726210/1‑151 (MQ=255)                                                                                      cAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTc                                                       >  1:1035094/1‑151 (MQ=255)                                                                                        gCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTCATCATTCTGCTCGACTCCATCACTCGTCTg                                                     >  1:98834/1‑151 (MQ=255)                                                                                         cGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTgg                                                    >  1:509185/1‑151 (MQ=255)                                                                                         cGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTgg                                                    >  2:939887/1‑151 (MQ=255)                                                                                           tCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcg                                                  >  1:303500/1‑151 (MQ=255)                                                                                             tGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcg                                                >  2:421941/1‑151 (MQ=255)                                                                                              ggTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcgc                                               >  2:393495/1‑151 (MQ=255)                                                                                               gTAAAAGGTGAAGTTGTTGCTTCTGCCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcg                                                <  1:177802/149‑1 (MQ=255)                                                                                               gTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcgcg                                              >  2:956986/1‑151 (MQ=255)                                                                                                tAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGgcgcgcgc                                             >  2:1187720/1‑151 (MQ=255)                                                                                                 aaaaGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGAACGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCt                                            >  1:257234/1‑151 (MQ=255)                                                                                                    aGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTac                                         >  1:481726/1‑151 (MQ=255)                                                                                                       agAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTacaac                                      >  1:1152723/2‑151 (MQ=255)                                                                                                        gAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTacaaca                                     <  1:295216/151‑1 (MQ=255)                                                                                                         aaGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAacac                                    >  1:136635/1‑151 (MQ=255)                                                                                                         aaGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAacac                                    <  2:387205/151‑1 (MQ=255)                                                                                                           gttgttGCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCg                                  <  1:197102/151‑1 (MQ=255)                                                                                                              gttgCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgtt                                <  1:830796/150‑1 (MQ=255)                                                                                                                 gCTTCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCgttgtt                             <  1:1029284/150‑1 (MQ=255)                                                                                                                   ttCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCg                          >  1:1063698/1‑151 (MQ=255)                                                                                                                    tCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCgg                         >  2:1026191/1‑151 (MQ=255)                                                                                                                    tCTACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCg                          <  2:1058511/150‑1 (MQ=255)                                                                                                                      tACCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCg                       >  1:588809/1‑151 (MQ=255)                                                                                                                       aCCTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGAAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGCCTGGCGCGCGCTTACAACACCGTTGTTCCg                          <  1:233766/147‑1 (MQ=255)                                                                                                                        ccTTTGACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTc                     >  1:470238/1‑151 (MQ=255)                                                                                                                             gACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCATAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGt                 >  1:317474/1‑150 (MQ=255)                                                                                                                             gACGAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGt                 <  1:718325/150‑1 (MQ=255)                                                                                                                                gAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTaa               >  1:565291/1‑149 (MQ=255)                                                                                                                                gAACCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAg             <  1:416749/151‑1 (MQ=255)                                                                                                                                 aaCCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAgt            <  2:438206/151‑1 (MQ=255)                                                                                                                                 aaCCCGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAgt            <  2:502072/151‑1 (MQ=255)                                                                                                                                   cccGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAgtgt          >  2:12848/1‑151 (MQ=255)                                                                                                                                    ccGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGtt         <  2:1062898/151‑1 (MQ=255)                                                                                                                                     cGCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTg        >  1:497331/1‑151 (MQ=255)                                                                                                                                      gCATCTCGCCGCGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGa       <  1:232691/151‑1 (MQ=255)                                                                                                                                            cgccgcGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCggt  <  1:356646/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  TGCTCAGAGCATTGCTTACAACCACCCGGATTGTGTGCTGATGGTTCTGCTGATCGACGAACGTCCGGAAGAAGTAACCGAGATGCAGCGTCTGGTAAAAGGTGAAGTTGTTGCTTCTACCTTTGACGAACCCGCATCTCGCCACGTTCAGGTTGCGGAAATGGTGATCGAGAAGGCCAAACGCCTGGTTGAGCACAAGAAAGACGTTATCATTCTGCTCGACTCCATCACTCGTCTGGCGCGCGCTTACAACACCGTTGTTCCGGCGTCAGGTAAAGTGTTGACCGGT  >  NC_000913/3966989‑3967277

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