Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,440,923	C→T	G36D (GGC→GAC)	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,923	0	C	T	100.0%	278.2 / NA	72	G36D (GGC→GAC)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (30/42);  total (30/42)

TGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTA  >  NC_000913/3440780‑3441066                                                                                                                                                |                                                                                                                                                tGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGcc                                                                                                                                          >  2:738723/1‑151 (MQ=255)   aGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCAc                                                                                                                                        <  2:76474/151‑1 (MQ=255)   aGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCAc                                                                                                                                        <  2:122586/151‑1 (MQ=255)    gCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACg                                                                                                                                       <  2:588497/151‑1 (MQ=255)    gCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACg                                                                                                                                       <  2:172521/151‑1 (MQ=255)     cAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGc                                                                                                                                      <  2:362398/151‑1 (MQ=255)         aTTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGctct                                                                                                                                   >  2:490153/1‑150 (MQ=255)             ccAGGATATCTTCCTTAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAg                                                                                                                               <  2:48370/150‑1 (MQ=255)             ccAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAgg                                                                                                                              <  2:696615/151‑1 (MQ=255)             ccAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAgg                                                                                                                              >  2:32264/1‑151 (MQ=255)               aGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCt                                                                                                                            <  2:32403/151‑1 (MQ=255)                 gATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTc                                                                                                                           >  1:466694/1‑150 (MQ=255)                    atCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAgg                                                                                                                       >  2:23865/1‑151 (MQ=255)                    atCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAgg                                                                                                                       >  1:589757/1‑151 (MQ=255)                        tCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCTCGCTCTAAAGGCTCAAGGGt                                                                                                                     <  1:674017/149‑1 (MQ=255)                        tCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGt                                                                                                                     <  2:950806/149‑1 (MQ=255)                         ccTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCa                                                                                                                   <  1:14148/150‑1 (MQ=255)                          cTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCAcc                                                                                                                 >  2:735417/1‑151 (MQ=255)                           tGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCt                                                                                                                >  2:691630/1‑151 (MQ=255)                               cGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                            <  2:505089/151‑1 (MQ=255)                                gCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcg                                                                                                           >  1:190077/1‑151 (MQ=255)                                 ccTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgt                                                                                                          <  2:792072/151‑1 (MQ=255)                                    tcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtgc                                                                                                        <  1:679660/150‑1 (MQ=255)                                    tcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtgc                                                                                                        <  1:20705/150‑1 (MQ=255)                                        tGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCg                                                                                                    >  1:695300/1‑150 (MQ=255)                                        tGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCg                                                                                                    >  2:532456/1‑150 (MQ=255)                                          gTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAAc                                                                                                 <  2:359068/151‑1 (MQ=255)                                               gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCa                                                                                              <  2:296146/149‑1 (MQ=255)                                               gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCAc                                                                                             <  1:730466/150‑1 (MQ=255)                                               gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACt                                                                                            >  2:750045/1‑151 (MQ=255)                                               gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACt                                                                                            <  2:763037/151‑1 (MQ=255)                                               gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACt                                                                                            >  2:552120/1‑151 (MQ=255)                                                    cATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCCAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGc                                                                                         <  1:118496/149‑1 (MQ=255)                                                       gtagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                    <  1:392897/151‑1 (MQ=255)                                                       gtagtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGat                                                                                    <  2:690338/151‑1 (MQ=255)                                                            acacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCaa                                                                               <  2:571139/151‑1 (MQ=255)                                                                 aTCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAgg                                                                          >  2:784607/1‑151 (MQ=255)                                                                    aaTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                       >  2:209659/1‑151 (MQ=255)                                                                    aaTCTCAACCTCGGTCACCGCGCAACCCAGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGATTATGGTCAAGGCCCCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCATCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                       <  1:982952/151‑1 (MQ=255)                                                                        tcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtt                                                                   >  1:89274/1‑151 (MQ=255)                                                                         cAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                      <  2:756784/147‑1 (MQ=255)                                                                         cAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGtta                                                                  <  1:198487/151‑2 (MQ=255)                                                                             cTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAg                                                               <  1:440485/150‑1 (MQ=255)                                                                                ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAc                                                           >  2:648229/1‑151 (MQ=255)                                                                                   cACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                         <  2:442994/150‑1 (MQ=255)                                                                                      cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTc                                                     >  2:780386/1‑151 (MQ=255)                                                                                       gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTca                                                    <  1:749756/151‑1 (MQ=255)                                                                                       gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTca                                                    <  2:193572/151‑1 (MQ=255)                                                                                       gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTca                                                    <  1:738671/151‑1 (MQ=255)                                                                                        cgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcac                                                   >  2:22225/1‑151 (MQ=255)                                                                                              ccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAcc                                             <  2:239967/151‑1 (MQ=255)                                                                                               cGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAccc                                            <  2:236705/151‑1 (MQ=255)                                                                                                   aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGc                                         >  1:299132/1‑150 (MQ=255)                                                                                                   aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGc                                         >  2:3150/1‑150 (MQ=255)                                                                                                     cGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCAt                                       <  1:864945/150‑1 (MQ=255)                                                                                                        tGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCGCTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATtgtg                                   >  2:726652/1‑151 (MQ=255)                                                                                                           gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTcc                                <  1:64944/151‑1 (MQ=255)                                                                                                           gagCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTcc                                >  2:123206/1‑151 (MQ=255)                                                                                                              cAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctct                             >  1:208160/1‑151 (MQ=255)                                                                                                               aGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctc                            >  2:244698/1‑151 (MQ=255)                                                                                                                gAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctct                           <  2:863942/151‑1 (MQ=255)                                                                                                                 aaTACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCtt                          <  2:24625/151‑1 (MQ=255)                                                                                                                            gTGCGTTACCCAGAGTATGGTCAAAGCCCCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAg                <  1:234294/150‑1 (MQ=255)                                                                                                                               cGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCttt            <  1:609607/151‑1 (MQ=255)                                                                                                                               cGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCttt            <  1:405157/151‑1 (MQ=255)                                                                                                                                gTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCCCCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTa           <  1:310723/151‑1 (MQ=255)                                                                                                                                  tACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt         >  2:501873/1‑151 (MQ=255)                                                                                                                                   aCCCAGAGTATGGTCAAAGCCGCGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt         <  1:733722/150‑1 (MQ=255)                                                                                                                                   aCCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACtt        >  2:791296/1‑151 (MQ=255)                                                                                                                                      cAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg      >  1:728428/1‑150 (MQ=255)                                                                                                                                        gagTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGt   >  1:669862/1‑151 (MQ=255)                                                                                                                                          gTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTa  >  1:275822/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                TGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGGAGTA  >  NC_000913/3440780‑3441066

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