Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,659,467	A→C	7.3%	E79A (GAA→GCA) ‡	mdtE →	multidrug efflux pump membrane fusion protein MdtE

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,659,467	0	A	C	7.3%	202.8 / 8.4	69	E79A (GAA→GCA) ‡	mdtE	multidrug efflux pump membrane fusion protein MdtE
Reads supporting (aligned to +/- strand):  ref base A (34/30);  new base C (3/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

CGCCGCCATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTG  >  NC_000913/3659318‑3659615                                                                                                                                                      |                                                                                                                                                     cgccgccATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGa                                                                                                                                                      <  1:913585/150‑1 (MQ=255)   ccgccATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAg                                                                                                                                                    <  2:314549/150‑1 (MQ=255)     gccATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGc                                                                                                                                                  <  1:1398369/150‑1 (MQ=255)            cGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAg                                                                                                                                           >  2:128253/1‑150 (MQ=255)            cGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAg                                                                                                                                           <  2:1387532/150‑1 (MQ=255)             gCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGt                                                                                                                                          <  1:1183294/150‑1 (MQ=255)                tGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGaa                                                                                                                                       >  1:416235/1‑150 (MQ=255)                 gAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAGttgcg                                                                                                                                      <  1:937546/150‑5 (MQ=255)                     tCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAgg                                                                                                                                  >  2:444853/1‑150 (MQ=255)                           tcgtcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAg                                                                                                                                           >  1:595425/1‑135 (MQ=255)                           tcgtcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGCCTTCGATAAAg                                                                                                                                           <  2:595425/135‑1 (MQ=255)                           tcgtcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGAtt                                                                                                                            <  2:1142013/150‑1 (MQ=255)                              tcACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGc                                                                                                                         <  2:982442/150‑1 (MQ=255)                                acacTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTg                                                                                                                       <  2:1144389/150‑1 (MQ=255)                                 cacTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGt                                                                                                                      <  1:372273/150‑1 (MQ=255)                                    tctcCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATc                                                                                                                   <  1:937146/150‑1 (MQ=255)                                       cccccGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGa                                                                                                                >  2:1126206/1‑150 (MQ=255)                                              ttCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATcc                                                                                                         <  1:1110157/150‑1 (MQ=255)                                              ttCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATcc                                                                                                         >  2:607226/1‑150 (MQ=255)                                                cGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGATTTACGTCCCCAGGTGGGCGGTATTATCATTAAACGAAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTg                                                                                                       <  2:400372/150‑1 (MQ=255)                                                cGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTg                                                                                                       >  2:819886/1‑150 (MQ=255)                                                 ggTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGc                                                                                                      >  1:1187916/1‑150 (MQ=255)                                                    cAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCCAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCATGCAcc                                                                                                   >  1:1356489/1‑150 (MQ=255)                                                          gtTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACa                                                                                             >  1:211213/1‑150 (MQ=255)                                                            tGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAgg                                                                                           >  2:1552947/1‑150 (MQ=255)                                                            tGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAgg                                                                                           >  2:457137/1‑150 (MQ=255)                                                            tGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAgg                                                                                           >  1:1072585/1‑150 (MQ=255)                                                             gAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGc                                                                                          >  1:548134/1‑150 (MQ=255)                                                             gAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGc                                                                                          >  2:754468/1‑150 (MQ=255)                                                              aGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGAGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGcc                                                                                         >  2:1348343/1‑150 (MQ=255)                                                               gCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCg                                                                                        <  1:128253/150‑1 (MQ=255)                                                                  aaTTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGc                                                                                     >  1:1375170/1‑150 (MQ=255)                                                                   aTTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCt                                                                                    <  2:416235/150‑1 (MQ=255)                                                                     tGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTaa                                                                                  <  1:244670/150‑1 (MQ=255)                                                                       cccGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAAc                                                                                <  1:1105010/150‑1 (MQ=255)                                                                       cccGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAAc                                                                                >  2:314534/1‑150 (MQ=255)                                                                         cGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTc                                                                              >  1:442196/1‑150 (MQ=255)                                                                             aGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGcc                                                                          <  1:559898/150‑1 (MQ=255)                                                                             aGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGcc                                                                          <  1:882211/150‑1 (MQ=255)                                                                             aGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGcc                                                                          <  1:444853/150‑1 (MQ=255)                                                                             aGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGcc                                                                          <  1:1126206/150‑1 (MQ=255)                                                                                       cTTATGAAGTTGCCGAGATACGTCCCGAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGc                                                                >  2:1387978/1‑150 (MQ=255)                                                                                       cTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGc                                                                >  2:1393701/1‑150 (MQ=255)                                                                                        ttATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCt                                                               >  2:438238/1‑150 (MQ=255)                                                                                        ttATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCt                                                               <  1:1652229/150‑1 (MQ=255)                                                                                         tATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTg                                                              >  2:221438/1‑150 (MQ=255)                                                                                           tGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGc                                                            >  1:1511611/1‑150 (MQ=255)                                                                                            gAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCg                                                           <  2:1375170/150‑1 (MQ=255)                                                                                                 tGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCggtaga                                                      >  2:889621/1‑145 (MQ=255)                                                                                                 tGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAgc                                                      >  2:315881/1‑150 (MQ=255)                                                                                                  gCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAgcg                                                     >  1:1035729/1‑150 (MQ=255)                                                                                                    cGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGct                                                   <  1:967524/150‑1 (MQ=255)                                                                                                         tACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTa                                              >  2:1300303/1‑150 (MQ=255)                                                                                                                 cAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAACGCGCTCTCTACCGCCAGc                                      >  2:1548130/1‑150 (MQ=255)                                                                                                                  aGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCa                                     <  2:1072585/150‑1 (MQ=255)                                                                                                                   ggTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCaa                                    <  1:1393701/150‑1 (MQ=255)                                                                                                                   ggTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCaa                                    <  1:1387978/150‑1 (MQ=255)                                                                                                                           gTATTATCATTAAACGCTAGTTTATCGAAGGCGATATAGTGAACGAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGACCTAGCTAAACTCCGCCAAAGGCTCGCTGGCTTAAGCGCTCTCTACCGCCAGCAATGCCCGCa                            <  2:204124/150‑1 (MQ=255)                                                                                                                           gTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCa                            <  2:109745/150‑1 (MQ=255)                                                                                                                               tatCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCAc                        >  1:549639/1‑150 (MQ=255)                                                                                                                                   aTTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCttt                    >  1:248523/1‑150 (MQ=255)                                                                                                                                     tAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTaa                  <  2:1187916/150‑1 (MQ=255)                                                                                                                                       aaCGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAAcc                >  2:1275503/1‑150 (MQ=255)                                                                                                                                           tcaCTTTATCGCCTTCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCa            >  1:31790/3‑150 (MQ=255)                                                                                                                                           cAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCa            <  1:315881/150‑1 (MQ=255)                                                                                                                                            caCTTTATCGCCTTCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAg           >  2:255802/2‑150 (MQ=255)                                                                                                                                                     gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGc              >  2:345726/1‑138 (MQ=255)                                                                                                                                                     gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGc              <  1:345726/138‑1 (MQ=255)                                                                                                                                                     gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTg  >  1:271175/1‑150 (MQ=255)                                                                                                                                                     gAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATAAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTg  >  2:324291/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CGCCGCCATGACGCCTGAGGTCGGTGTCGTCACACTCTCCCCCGGTTCGGTCAATGTGTTGAGCGAATTGCCCGGTAGAACCGTTCCTTATGAAGTTGCCGAGATACGTCCCCAGGTGGGCGGTATTATCATTAAACGCAACTTTATCGAAGGCGATAAAGTGAACCAGGGCGATTCGCTGTATCAGATTGATCCTGCACCTTTACAGGCCGAGCTAAACTCCGCCAAAGGCTCGCTGGCGAAAGCGCTCTCTACCGCCAGCAATGCCCGCATCACCTTTAACCGCCAGGCATCGTTG  >  NC_000913/3659318‑3659615

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