Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	3,440,086	C→A	G315V (GGC→GTC)	rpoA ←	RNA polymerase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	3,440,086	0	C	A	100.0%	255.1 / NA	72	G315V (GGC→GTC)	rpoA	RNA polymerase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (40/32);  total (40/32)

TGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGC  >  NC_000913_3_bme_pgi/3439953‑3440231                                                                                                                                      |                                                                                                                                                  tGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACg                                                                                                                                   <  1:509754/150‑1 (MQ=255)  ggCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACgg                                                                                                                                  >  2:169432/1‑150 (MQ=255)    ctgctgTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTTGCAAGTTTTCCAGGCGCATGACCAGTGACAGTCCACGGGa                                                                                                                                <  2:382510/150‑1 (MQ=255)     tgctgTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAAAGTGTGATTAGGCCAGCCGATAGCGAGGCTTTCCGGGGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGaa                                                                                                                               >  1:598647/1‑150 (MQ=255)        tgTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGAc                                                                                                                                           >  1:234485/1‑135 (MQ=255)        tgTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGAc                                                                                                                                           <  2:234485/135‑1 (MQ=255)           tGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAACCCTTAGCCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGGCAGTTTTCCACGTGCATGACCAGAGACAGTCCACGGGAAGCCAGc                                                                                                                         <  2:562856/150‑1 (MQ=255)           tGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGc                                                                                                                         <  2:107132/150‑1 (MQ=255)                   aGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCttt                                                                                                                 >  1:617445/1‑150 (MQ=255)                   aGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCttt                                                                                                                 >  2:208039/1‑150 (MQ=255)                    gTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTa                                                                                                                >  1:528170/1‑150 (MQ=255)                      tGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAt                                                                                                              <  1:208039/150‑1 (MQ=255)                                tctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAaga                                                                                                    >  2:140835/1‑150 (MQ=255)                                tctTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAaga                                                                                                    >  4:376275/1‑150 (MQ=255)                                   tACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCag                                                                                                                                               >  2:217893/1‑104 (MQ=255)                                   tACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCag                                                                                                                                               <  1:217893/104‑1 (MQ=255)                                      gATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGCGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtttt                                                                                              >  3:361356/1‑150 (MQ=255)                                       aTGGCGCATGACCTTATCCTTCTCAGTAAAGCCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAttttt                                                                                             <  2:598647/150‑1 (MQ=255)                                        tGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtttttt                                                                                            >  2:272097/1‑150 (MQ=255)                                         ggCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTa                                                                                           <  1:169432/150‑1 (MQ=255)                                          gcgcATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                                          <  2:796322/134‑1 (MQ=255)                                          gcgcATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCa                                                                                                          >  1:796322/1‑134 (MQ=255)                                           cgcATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTAcc                                                                                         >  1:4852/1‑150 (MQ=255)                                                gACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGt                                                                                    >  1:727615/1‑150 (MQ=255)                                                  ccTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTa                                                                                  <  3:376275/150‑1 (MQ=255)                                                  ccTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTa                                                                                  <  4:494723/150‑1 (MQ=255)                                                    ttATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAgg                                                                                <  1:32036/150‑1 (MQ=255)                                                      aTCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCg                                                                              >  2:758403/1‑150 (MQ=255)                                                          ttCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtc                                                                                                             <  3:461709/115‑1 (MQ=255)                                                          ttCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtc                                                                                                             >  4:461709/1‑115 (MQ=255)                                                          ttCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGtttt                                                                          >  2:777982/1‑150 (MQ=255)                                                           tctcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTaa                                                                                                               <  2:298618/112‑1 (MQ=255)                                                           tctcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCGTTaa                                                                                                               >  1:298618/1‑112 (MQ=255)                                                            ctcAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTaa                                                                        <  2:528170/150‑1 (MQ=255)                                                                gTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGCAGCCAGCACGTCTTTAATCTCAGTAAGAGAtttttt                                                                                            >  2:252307/1‑126 (MQ=255)                                                                gTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAtttttt                                                                                            <  1:252307/126‑1 (MQ=255)                                                                gTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAg                                                                    >  2:97095/1‑150 (MQ=255)                                                                 tAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCa                                                                                                                           >  2:196309/1‑94 (MQ=255)                                                                 tAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCa                                                                                                                           <  1:196309/94‑1 (MQ=255)                                                                 tAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAttttt                                                                                             >  1:463669/1‑124 (MQ=255)                                                                 tAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGAttttt                                                                                             <  2:463669/124‑1 (MQ=255)                                                                      ccTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAAc                                                              <  1:140835/150‑1 (MQ=255)                                                                            ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAg                                                                                                                                         <  3:159041/69‑1 (MQ=255)                                                                            ccTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAg                                                                                                                                         >  4:159041/1‑69 (MQ=255)                                                                                  aTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGttt                                                                           <  1:355904/125‑1 (MQ=255)                                                                                  aTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGttt                                                                           >  2:355904/1‑125 (MQ=255)                                                                                  aTCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTg                                                  >  1:609708/1‑150 (MQ=255)                                                                                    ccGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTa                                                >  3:57108/1‑150 (MQ=255)                                                                                     cGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTAc                                               <  1:97095/150‑1 (MQ=255)                                                                                     cGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTAc                                               <  3:209474/150‑1 (MQ=255)                                                                                       gTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCa                                             >  4:167195/1‑150 (MQ=255)                                                                                              gTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCAcc                                      >  1:297110/1‑150 (MQ=255)                                                                                                  gCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGata                                  >  2:852086/1‑150 (MQ=255)                                                                                                            ccGGTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAg                        >  2:167457/1‑150 (MQ=255)                                                                                                               gTGGCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCATCCTCGGTACGCTGTACCAGATCTCCGATATAGTGGATAGCtt                     <  4:361356/150‑1 (MQ=255)                                                                                                               gTGGCCAGTTTGCCAGGCGCATGACCAGAGACCGTCCAAGGGAAGCCAGCACGTCTTTAATCTCAGTCAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGCAGCTCCACCTCGGTACGCTGTACCAGATCACCGATATAGCGGATCGCtt                     >  2:469539/1‑150 (MQ=255)                                                                                                                  gCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAaga                                                                                                    >  4:290042/1‑68 (MQ=255)                                                                                                                  gCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAaga                                                                                                    <  3:290042/68‑1 (MQ=255)                                                                                                                  gCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGa                           >  3:198728/1‑141 (MQ=255)                                                                                                                  gCCAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGa                           <  4:198728/141‑1 (MQ=255)                                                                                                                   ccAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGc                 >  2:325919/1‑150 (MQ=255)                                                                                                                    cAGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCt                >  1:350416/1‑150 (MQ=255)                                                                                                                     aGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtct                                                                                                            <  4:521178/57‑1 (MQ=255)                                                                                                                     aGTTTTCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAAtct                                                                                                            >  3:521178/1‑57 (MQ=255)                                                                                                                       ttttCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa             <  2:4852/150‑1 (MQ=255)                                                                                                                       ttttCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa             >  1:186257/1‑150 (MQ=255)                                                                                                                       ttttCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTa             <  1:167457/150‑1 (MQ=255)                                                                                                                         ttCCAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAg           <  2:727615/150‑1 (MQ=255)                                                                                                                            cAGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCa        >  1:384035/1‑150 (MQ=255)                                                                                                                             aGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAg       >  4:460613/1‑150 (MQ=255)                                                                                                                             aGGCGCATGACCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCAATATAGTGGATAGCTTCTGCTTTAAGGCAg       >  1:553884/1‑150 (MQ=255)                                                                                                                                  cATGACAAGAGACCGCCCACGGGCCGCCAGTTCATCTTTGTCCTCGGGAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGc  <  3:492582/150‑1 (MQ=255)                                                                                                                                      |                                                                                                                                                  TGGCTGCTGTTGCGGTTCAGTTGACGACCACTCTTACGATGGCGCATGACCTTATCCTTCTCAGTAAAACCTTAACCTGTGATCCGGTTACTCGTCAGCGATGCTTGCCGGTGGCCAGTTTTCCAGGCGCATGCCCAGAGACAGTCCACGGGAAGCCAGCACGTCTTTAATCTCAGTAAGAGATTTTTTACCAAGGTTAGGCGTTTTAAGGAGCTCAACCTCGGTACGCTGTACCAGATCACCGATATAGTGGATAGCTTCTGCTTTAAGGCAGTTAGC  >  NC_000913_3_bme_pgi/3439953‑3440231

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑