Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,186,274	A→T	N309Y (AAC→TAC)	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,274	0	A	T	100.0%	249.2 / NA	65	N309Y (AAC→TAC)	rpoC	RNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (23/42);  total (23/42)

TTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGACTACTCCGGTCGTTCTGTAATCACCGTAG  >  NC_000913/4186126‑4186421                                                                                                                                                     |                                                                                                                                                    tttCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTAc                                                                                                                                                   >  2:1010530/1‑151 (MQ=255)   tCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACgg                                                                                                                                                 >  1:673726/1‑151 (MQ=255)    cgcgACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACgg                                                                                                                                                 <  2:542236/150‑1 (MQ=255)        aCTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGgtcgtc                                                                                                                                            <  1:994519/151‑1 (MQ=255)        aCTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGgtcgtc                                                                                                                                            <  2:533682/151‑1 (MQ=255)           tCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGtcgtcg                                                                                                                                           <  1:589551/149‑1 (MQ=255)           tCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGtcgtcg                                                                                                                                           <  1:448906/149‑1 (MQ=255)              gACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTc                                                                                                                                      <  1:613760/151‑1 (MQ=255)               aCCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCg                                                                                                                                     >  2:1033094/1‑151 (MQ=255)                ccTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGt                                                                                                                                    <  1:588097/151‑1 (MQ=255)                ccTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGt                                                                                                                                    <  1:480822/151‑1 (MQ=255)                   gAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGc                                                                                                                                  <  2:1055298/150‑1 (MQ=255)                   gAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCg                                                                                                                                 >  2:536486/1‑151 (MQ=255)                   gAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCg                                                                                                                                 >  1:167863/1‑151 (MQ=255)                    aaCGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCg                                                                                                                                 >  1:745353/1‑150 (MQ=255)                            gTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATc                                                                                                                              <  2:796570/145‑1 (MQ=255)                              aTCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTc                                                                                                                      <  1:185738/151‑1 (MQ=255)                                cgtcgCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCt                                                                                                                     <  2:790768/150‑1 (MQ=255)                                cgtcgCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCt                                                                                                                     >  2:55403/1‑150 (MQ=255)                                    gcgTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACaa                                                                                                                <  2:264719/151‑1 (MQ=255)                                      gTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAg                                                                                                               >  2:1050855/1‑150 (MQ=255)                                        cATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGt                                                                                                            <  1:946061/151‑1 (MQ=255)                                         aTTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTc                                                                                                           <  1:467745/151‑1 (MQ=255)                                            aaCCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCct                                                                                                         <  1:352112/150‑1 (MQ=255)                                             aCCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCctct                                                                                                       <  1:911309/151‑1 (MQ=255)                                                gTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGa                                                                                                     >  2:57553/1‑150 (MQ=255)                                                  aacaacCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGaaa                                                                                                   <  2:281045/150‑1 (MQ=255)                                                  aacaacCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAAt                                                                                                  >  2:536751/1‑151 (MQ=255)                                                     aacCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCTGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACACCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGaaa                                                                                                   <  1:940666/147‑1 (MQ=255)                                                       ccGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTg                                                                                             <  1:680860/151‑1 (MQ=255)                                                       ccGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTg                                                                                             <  1:133749/151‑1 (MQ=255)                                                             gAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCg                                                                                         <  1:10107/149‑1 (MQ=255)                                                             gAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGAc                                                                                       <  1:312261/151‑1 (MQ=255)                                                                  gTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGAt                                                                                  >  1:937408/1‑151 (MQ=255)                                                                   tCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGAt                                                                                  <  1:982917/150‑1 (MQ=255)                                                                      gctgGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCaa                                                                               <  2:86166/150‑1 (MQ=255)                                                                       ctgGATTTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATc                                                                                 >  1:529480/1‑147 (MQ=255)                                                                           aTCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTa                                                                          >  2:498957/1‑150 (MQ=255)                                                                               ggCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTaaa                                                                        <  2:457884/148‑1 (MQ=255)                                                                                gCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAgg                                                                    <  1:201831/151‑1 (MQ=255)                                                                                 cTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAggg                                                                   >  2:1056500/1‑151 (MQ=255)                                                                                   gcgcCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAggg                                                                   <  1:57533/149‑1 (MQ=255)                                                                                      ccGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGtt                                                              <  2:674487/151‑1 (MQ=255)                                                                                       cGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGttt                                                             <  2:301365/151‑1 (MQ=255)                                                                                       cGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGttt                                                             >  2:876902/1‑151 (MQ=255)                                                                                          aCATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCg                                                          >  2:828768/1‑151 (MQ=255)                                                                                            atcatcGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTc                                                        <  1:379905/151‑1 (MQ=255)                                                                                               atcGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTc                                                        >  1:12809/1‑148 (MQ=255)                                                                                                    acgtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCt                                                 <  1:536717/150‑1 (MQ=255)                                                                                                    acgtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCt                                                 <  2:708786/150‑1 (MQ=255)                                                                                                    acgtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAAAAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCt                                                 <  1:622071/150‑1 (MQ=255)                                                                                                      gtaACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCt                                              <  1:55077/151‑1 (MQ=255)                                                                                                        aaCGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCg                                            >  1:145733/1‑151 (MQ=255)                                                                                                         aCGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCgg                                           <  2:843466/151‑1 (MQ=255)                                                                                                         aCGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCgg                                           >  1:521405/1‑151 (MQ=255)                                                                                                            aaaaaCGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTaa                                        <  1:917825/151‑1 (MQ=255)                                                                                                              aaaCGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGc                                      >  2:1071542/1‑151 (MQ=255)                                                                                                                aCGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCgt                                    <  1:1043564/151‑1 (MQ=255)                                                                                                                 cGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCgtg                                   <  2:632091/151‑1 (MQ=255)                                                                                                                    aTGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGtt                                 <  2:589777/150‑1 (MQ=255)                                                                                                                    aTGCTGCAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTg                                >  1:447233/1‑151 (MQ=255)                                                                                                                      gctgcAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGac                              <  1:426085/151‑1 (MQ=255)                                                                                                                        tgcAGGAAGCGGTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGact                             >  2:1005063/1‑150 (MQ=255)                                                                                                                                  ggTAGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGACTACTCCGGTCGt                  <  1:570475/151‑1 (MQ=255)                                                                                                                                     aGACGCCCTGCTGGATTACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGACTACTCCGGTCGTTCt               >  2:332773/1‑151 (MQ=255)                                                                                                                                                   aaaaCGGTCGTCGCGGTCGTGCGATAACCGGTGCTAACAAGCGTCCTCTGAAGTCTTTGGCCGACATGATCAAAGGAAACAAGGGTCGTTTCCGTAAGAACCTGCGCGGTAAGCGTTGTGACTAGTCCGGTCGTTCTGTAAACACCGTAg  >  1:991489/3‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TTTCGCGACTTCTGACCTGAACGATCTGTATCGTCGCGTCATTAACCGTAACAACCGTCTGAAACGTCTGCTGGATCTGGCTGCGCCGGACATCATCGTACGTAACGAAAAACGTATGCTGCAGGAAGCGGTAGACGCCCTGCTGGATAACGGTCGTCGCGGTCGTGCGATCACCGGTTCTAACAAGCGTCCTCTGAAATCTTTGGCCGACATGATCAAAGGTAAACAGGGTCGTTTCCGTCAGAACCTGCTCGGTAAGCGTGTTGACTACTCCGGTCGTTCTGTAATCACCGTAG  >  NC_000913/4186126‑4186421

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