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%	289.8 / NA	76	G36D (GGC→GAC)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/39);  total (37/39)

GGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGG  >  NC_000913/3440777‑3441062                                                                                                                                                   |                                                                                                                                            ggTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCaa                                                                                                                                          >  2:990274/1‑150 (MQ=255)  gTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCaa                                                                                                                                          >  1:699937/1‑149 (MQ=255)     gAGCAGGATTTCCAGGATCTCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGcc                                                                                                                                      >  2:154621/1‑150 (MQ=255)     gAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGcc                                                                                                                                      >  1:464095/1‑150 (MQ=255)      aGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCAc                                                                                                                                    >  1:504866/1‑151 (MQ=255)        cAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGc                                                                                                                                  >  1:908473/1‑151 (MQ=255)        cAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAAACTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGCGCGTTACCCAGAGTATGGTCAAAGCCAc                                                                                                                                    >  2:782899/1‑149 (MQ=255)           gATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGctct                                                                                                                               >  2:1016211/1‑151 (MQ=255)               tCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTaaa                                                                                                                            <  2:326241/150‑1 (MQ=255)                  aGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCt                                                                                                                        >  2:307902/1‑151 (MQ=255)                  aGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCt                                                                                                                        >  2:783180/1‑151 (MQ=255)                      tatCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCa                                                                                                                      >  1:425642/1‑149 (MQ=255)                       atCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAg                                                                                                                    <  1:245580/150‑1 (MQ=255)                          ttCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACAACATCAATCTCAACCTCGGTCACCGCGCAACCCAGCATCGATGAGAGCAGAATAAGGCGCAGTGCGTTACCCAGAGCATGGTCAAAGACACGCTCTAAAGGCTCAAGGGt                                                                                                                 >  1:328381/1‑150 (MQ=255)                            ccTGAACGCCGTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCa                                                                                                               <  1:908982/150‑1 (MQ=255)                                  cGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgc                                                                                                        <  2:568092/151‑1 (MQ=255)                                   gCCTTCTTTGTGGCTGGACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGGGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcg                                                                                                       >  1:408281/1‑151 (MQ=255)                                   gCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcg                                                                                                       >  1:628914/1‑151 (MQ=255)                                     cttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtg                                                                                                     >  2:909287/1‑151 (MQ=255)                                      ttcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtgc                                                                                                    >  2:1009669/1‑151 (MQ=255)                                       tcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtg                                                                                                     <  2:777180/149‑1 (MQ=255)                                       tcttTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGgcgtg                                                                                                     <  1:107685/149‑1 (MQ=255)                                          ttGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTc                                                                                                 <  1:341231/150‑1 (MQ=255)                                          ttGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTc                                                                                                 <  2:706146/150‑1 (MQ=255)                                           tGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCg                                                                                                <  2:1009980/150‑1 (MQ=255)                                           tGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCg                                                                                                <  1:611516/150‑1 (MQ=255)                                            ggTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGa                                                                                               <  1:547505/150‑1 (MQ=255)                                             gTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAAc                                                                                             <  2:949695/151‑1 (MQ=255)                                               gCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTc                                                                                           >  2:293688/1‑151 (MQ=255)                                                 tGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCAc                                                                                         >  1:457029/1‑151 (MQ=255)                                                  gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGCTGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCAc                                                                                         >  1:852759/1‑150 (MQ=255)                                                  gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAGGGGTCACCTTGGCGTGCGTCGAACTCACt                                                                                        <  2:906032/151‑1 (MQ=255)                                                  gTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACt                                                                                        <  2:463898/151‑1 (MQ=255)                                                            agtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGatat                                                                              >  2:71348/1‑151 (MQ=255)                                                            agtaCACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGata                                                                               <  2:504873/150‑1 (MQ=255)                                                               acacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCa                                                                            >  2:881633/1‑150 (MQ=255)                                                                cacCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAc                                                                          >  1:245157/1‑151 (MQ=255)                                                                  ccATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAAcc                                                                         >  2:67655/1‑150 (MQ=255)                                                                   cATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAg                                                                       <  1:250396/151‑1 (MQ=255)                                                                      cAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgc                                                                     <  1:588663/150‑1 (MQ=255)                                                                       aaTCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgc                                                                   >  2:512274/1‑151 (MQ=255)                                                                         tctcAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCgg                                                                 >  1:122539/1‑151 (MQ=255)                                                                            cAACCTCGGTCTCCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGgcgca                                                                  <  1:220741/147‑2 (MQ=255)                                                                            cAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGcgcg                                                                  <  1:102823/147‑1 (MQ=255)                                                                            cAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGttt                                                              <  1:121583/151‑1 (MQ=255)                                                                                  cGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGa                                                          <  2:163729/149‑1 (MQ=255)                                                                                   ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAc                                                       >  2:125544/1‑151 (MQ=255)                                                                                   ggTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAc                                                       <  1:775813/151‑1 (MQ=255)                                                                                    gTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAAct                                                      <  1:301614/151‑1 (MQ=255)                                                                                      cACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAActc                                                     <  2:130737/150‑1 (MQ=255)                                                                                         cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTc                                                 <  1:253250/151‑1 (MQ=255)                                                                                         cgcgCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTc                                                 <  1:74533/151‑1 (MQ=255)                                                                                          gcgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTca                                                <  1:239020/151‑1 (MQ=255)                                                                                           cgcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcac                                               <  1:906848/151‑1 (MQ=255)                                                                                            gcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcaca                                              <  2:563897/151‑1 (MQ=255)                                                                                            gcAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTcac                                               >  2:900965/1‑150 (MQ=255)                                                                                                cccGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAAc                                          <  2:156357/151‑1 (MQ=255)                                                                                                   ggCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCt                                       <  2:130301/151‑1 (MQ=255)                                                                                                   ggCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCt                                       <  2:8907/151‑1 (MQ=255)                                                                                                    gCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTg                                      <  2:48817/151‑1 (MQ=255)                                                                                                      aTCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGc                                     >  2:365480/1‑150 (MQ=255)                                                                                                                gCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctc                          >  2:748778/1‑151 (MQ=255)                                                                                                                 cAGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctct                         <  2:1015530/151‑1 (MQ=255)                                                                                                                  aGAATACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCctctc                        <  2:212745/151‑1 (MQ=255)                                                                                                                     aTACGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCttt                     >  2:604215/1‑151 (MQ=255)                                                                                                                       aCGGCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTgg                   >  2:893/1‑151 (MQ=255)                                                                                                                         ggCGCAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTa                 <  1:531614/151‑1 (MQ=255)                                                                                                                             cAGTGCGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTaa             <  2:709503/151‑1 (MQ=255)                                                                                                                                  cGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCttt        >  2:573399/1‑151 (MQ=255)                                                                                                                                  cGTTACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCttt        >  2:144736/1‑151 (MQ=255)                                                                                                                                    ttACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTAc      >  2:741968/1‑151 (MQ=255)                                                                                                                                     tACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt     >  2:636494/1‑151 (MQ=255)                                                                                                                                     tACCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACt     >  2:548804/1‑151 (MQ=255)                                                                                                                                      aCCCAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTAGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACtt    <  2:408286/151‑1 (MQ=255)                                                                                                                                       cccAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTg   <  1:1016196/151‑1 (MQ=255)                                                                                                                                        ccAGAGTATGGTCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTgg  >  2:17043/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                            GGTTGAGCAGGATTTCCAGGATATCTTCCTGAACGCCTTCTTTGGTGCTGTACTCATGTAGTACACCATCAATCTCAACCTCGGTCACCGCGCAACCCGGCATCGATGAGAGCAGAATACGGCGCAGTGCGTTACCCAGAGTATGGCCAAAGCCACGCTCTAAAGGCTCAAGGGTCACCTTGGCGTGCGTCGAACTCACTTGCTCGATATCAACCAGGCGCGGTTTTAGAAACTCTGTCACAGAACCCTGCATTGTGTCCTCTCTTTGGTACTAAGCTTTACTTGG  >  NC_000913/3440777‑3441062

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