Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,440,950	G→A	T27I (ACC→ATC)	rpoA ←	RNA polymerase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,440,950	0	G	A	100.0%	312.0 / NA	96	T27I (ACC→ATC)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (47/49);  total (47/49)

GCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGT  >  NC_000913/3440823‑3441083                                                                                                                                |                                                                                                                                      gCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcg                                                                                                                            <  2:208738/139‑1 (MQ=255)     tACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgtgcg                                                                                                                        >  1:158168/1‑139 (MQ=255)          aggtagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGTCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGaa                                                                                                                   >  2:775625/3‑139 (MQ=255)               gtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              <  1:775625/139‑1 (MQ=255)               gtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              >  2:1228440/1‑139 (MQ=255)               gtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              <  2:158168/139‑1 (MQ=255)                  cacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTg                                                                                                           >  1:548533/1‑139 (MQ=255)                  cacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTg                                                                                                           >  1:513216/1‑139 (MQ=255)                    ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgt                                                                                                                           >  2:1066567/1‑121 (MQ=255)                    ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcgt                                                                                                                           <  1:1066567/121‑1 (MQ=255)                    ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  1:111690/139‑1 (MQ=255)                    ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  1:358916/139‑1 (MQ=255)                     cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcg                                                                                                                            <  2:106231/119‑1 (MQ=255)                     cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcg                                                                                                                            <  2:337791/119‑1 (MQ=255)                     cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcg                                                                                                                            >  1:337791/1‑119 (MQ=255)                     cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGgcg                                                                                                                            >  1:106231/1‑119 (MQ=255)                     cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTc                                                                                                        >  2:1067128/1‑139 (MQ=255)                       tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              >  2:794552/1‑131 (MQ=255)                       tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCAc                                                                                                              <  1:794552/131‑1 (MQ=255)                       tCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGa                                                                                                      >  2:685794/1‑139 (MQ=255)                        cAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                                     >  2:1273203/1‑139 (MQ=255)                         aaTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                                                                                                    <  1:1058894/139‑1 (MQ=255)                          aTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                                                   >  1:1158862/1‑139 (MQ=255)                          aTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                                                   >  1:357350/1‑139 (MQ=255)                               aaCCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc                                                                                              >  1:130885/1‑139 (MQ=255)                                 ccTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  1:1067227/1‑126 (MQ=255)                                 ccTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  2:1067227/126‑1 (MQ=255)                                 ccTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAg                                                                                            >  1:22688/1‑139 (MQ=255)                                 ccTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAg                                                                                            >  1:609219/1‑139 (MQ=255)                                 ccTCGGTCACCGCGAAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCCAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAg                                                                                            >  1:1091453/1‑139 (MQ=255)                                  cTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAgg                                                                                           >  2:1245328/1‑139 (MQ=255)                                    cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         <  2:609219/139‑1 (MQ=255)                                    cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         <  2:1091453/139‑1 (MQ=255)                                     ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                                        >  2:410552/1‑139 (MQ=255)                                     ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                                        >  1:704473/1‑139 (MQ=255)                                      gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  2:691089/121‑1 (MQ=255)                                      gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  1:691089/1‑121 (MQ=255)                                        cACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGt                                                                                     <  2:130885/139‑1 (MQ=255)                                                     ggCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTg                                                                        <  2:513216/139‑1 (MQ=255)                                                          cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACGTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  2:843467/139‑1 (MQ=255)                                                          cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  1:479191/139‑1 (MQ=255)                                                          cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  2:22688/139‑1 (MQ=255)                                                          cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCCCGGTTTTAGAAACTCTGTcaca                                                                   <  1:685794/139‑1 (MQ=255)                                                              gagaGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAc                                                               <  1:237100/139‑1 (MQ=255)                                                              gagaGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAc                                                               <  1:1273203/139‑1 (MQ=255)                                                                 agCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCt                                                            >  1:957049/1‑139 (MQ=255)                                                                          cGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTgtgt                                                   <  1:1228440/139‑1 (MQ=255)                                                                             cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCCAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                >  2:875631/1‑139 (MQ=255)                                                                             cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                                        <  2:47155/99‑1 (MQ=255)                                                                             cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                                        >  1:47155/1‑99 (MQ=255)                                                                             cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                <  2:957049/139‑1 (MQ=255)                                                                             cgcAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCct                                                <  1:1245328/139‑1 (MQ=255)                                                                                aGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctc                                             >  1:540474/1‑139 (MQ=255)                                                                                  tGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCtt                                           <  1:1067128/139‑1 (MQ=255)                                                                                  tGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCtt                                           <  2:1158862/139‑1 (MQ=255)                                                                                       tACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTa                                      <  2:704473/139‑1 (MQ=255)                                                                                              agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       <  2:701729/83‑1 (MQ=255)                                                                                              agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                                       >  1:701729/1‑83 (MQ=255)                                                                                              agTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCt                               <  1:410552/139‑1 (MQ=255)                                                                                                tATGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCttt                             >  2:809954/1‑139 (MQ=255)                                                                                                 aTGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCt                               <  1:937690/136‑1 (MQ=255)                                                                                                 aTGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCt                               >  2:937690/1‑136 (MQ=255)                                                                                                  tGGCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGATCCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTAc                           >  2:826125/1‑139 (MQ=255)                                                                                                   ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   >  1:606796/1‑98 (MQ=255)                                                                                                   ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                                                   <  2:606796/98‑1 (MQ=255)                                                                                                   ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTACTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  <  1:1086284/115‑1 (MQ=255)                                                                                                   ggCCAAAGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTACTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  >  2:1086284/1‑115 (MQ=255)                                                                                                        aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtttt                                                                                  <  1:1196752/78‑1 (MQ=255)                                                                                                        aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtttt                                                                                  >  2:1196752/1‑78 (MQ=255)                                                                                                        aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                     >  2:16822/1‑139 (MQ=255)                                                                                                        aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                     >  1:729226/1‑139 (MQ=255)                                                                                                        aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                     >  2:51391/1‑139 (MQ=255)                                                                                                        aaGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAg                     >  1:1065954/1‑139 (MQ=255)                                                                                                         aGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                                          >  1:10080/1‑85 (MQ=255)                                                                                                         aGCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                                          <  2:10080/85‑1 (MQ=255)                                                                                                          gCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTa                   <  2:1128607/139‑1 (MQ=255)                                                                                                          gCCACGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTa                   <  1:826125/139‑1 (MQ=255)                                                                                                              cGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         <  2:403334/65‑1 (MQ=255)                                                                                                              cGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         >  1:403334/1‑65 (MQ=255)                                                                                                              cGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  <  2:637704/104‑1 (MQ=255)                                                                                                              cGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  >  1:637704/1‑104 (MQ=255)                                                                                                              cGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGc               >  1:1284877/1‑139 (MQ=255)                                                                                                              cGCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGc               >  2:237978/1‑139 (MQ=255)                                                                                                               gCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACATAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTGGGAGTAAAGCt              <  1:875631/139‑1 (MQ=255)                                                                                                               gCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCt              <  1:809954/139‑1 (MQ=255)                                                                                                               gCTCTAAAGGCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCCCTCTTTGGTACTAAGATTTACCTGGAGTAAAGCt              <  1:166453/139‑1 (MQ=255)                                                                                                                       ggCTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATc      <  2:1065954/139‑1 (MQ=255)                                                                                                                         cTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  <  2:495946/93‑1 (MQ=255)                                                                                                                         cTCAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTc                                                  >  1:495946/1‑93 (MQ=255)                                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         <  1:705736/36‑1 (MQ=255)                                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCt                                                                                                         >  2:705736/1‑36 (MQ=255)                                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAccc                                                             <  2:612132/80‑1 (MQ=255)                                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAccc                                                             >  1:612132/1‑80 (MQ=255)                                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt  <  2:234246/139‑1 (MQ=255)                                                                                                                           cAAGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGgt  >  2:160875/1‑139 (MQ=255)                                                                                                                             aGGATCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcg                                                                                         <  2:837311/50‑1 (MQ=255)                                                                                                                                |                                                                                                                                      GCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTAAAGCTCGACGATCAGGT  >  NC_000913/3440823‑3441083

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