Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,157,564	A→G	6.7%	T768A (ACC→GCC)	mdtB →	multidrug efflux pump RND permease subunit MdtB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,157,564	0	A	G	6.7%	173.6 / 3.3	60	T768A (ACC→GCC)	mdtB	multidrug efflux pump RND permease subunit MdtB
Reads supporting (aligned to +/- strand):  ref base A (30/26);  new base G (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

CGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTGGATCAGTTCCCGGTAAC  >  NC_000913/2157415‑2157709                                                                                                                                                      |                                                                                                                                                  cGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCacaaca                                                                                                                                                   <  1:483515/150‑1 (MQ=255)     tGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGa                                                                                                                                               >  2:463290/1‑150 (MQ=255)                  cgcgACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGc                                                                                                                                  >  2:1241052/1‑150 (MQ=255)                  cgcgACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGc                                                                                                                                  >  1:1354949/1‑150 (MQ=255)                   gcgACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGcc                                                                                                                                 >  2:1097050/1‑150 (MQ=255)                     gACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTc                                                                                                                               >  1:239048/1‑150 (MQ=255)                       cAGCGCCAGCCGTCTTTGGATCAGCATTGAGGATGTAGATAACTCCCTGTACAAATCTTTTGTTCAGCGGCTGATTTCAACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGATCACAACACCGAAAATACCCCAGGACTcga                                                                                                                             <  1:925605/150‑2 (MQ=255)                             cAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCt                                                                                                                       <  2:703824/150‑1 (MQ=255)                                ccGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGa                                                                                                                    <  2:1064281/150‑1 (MQ=255)                                ccGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACAATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGa                                                                                                                    <  2:37293/150‑1 (MQ=255)                                    cTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATAcc                                                                                                                <  2:70767/150‑1 (MQ=255)                                       gggATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCAtt                                                                                                             >  1:1142417/1‑150 (MQ=255)                                        ggATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACGCCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTc                                                                                                            >  2:511917/1‑150 (MQ=255)                                        ggATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACGCCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTc                                                                                                            >  2:516866/1‑150 (MQ=255)                                         gATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCg                                                                                                           <  1:463290/150‑1 (MQ=255)                                              gCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGa                                                                                                      >  2:1125908/1‑150 (MQ=255)                                              gCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGa                                                                                                      >  2:1024105/1‑150 (MQ=255)                                                  ggCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGAt                                                                                                                   >  1:1165485/1‑133 (MQ=255)                                                  ggCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGAt                                                                                                                   <  2:1165485/133‑1 (MQ=255)                                                   gCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACcagc                                                                                                 >  2:1192589/1‑150 (MQ=255)                                                               aaCGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGcc                                                                                                         >  2:425980/1‑130 (MQ=255)                                                               aaCGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGcc                                                                                                         <  1:425980/130‑1 (MQ=255)                                                                  gCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCgtg                                                                                  >  1:621097/1‑150 (MQ=255)                                                                  gCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCgtg                                                                                  <  1:1241052/150‑1 (MQ=255)                                                                   cccTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCgtgg                                                                                 >  2:1087555/1‑150 (MQ=255)                                                                      tGTACAACGCGTTTGGTCTGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGc                                                                              <  1:613660/150‑1 (MQ=255)                                                                      tGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGc                                                                              <  1:516615/150‑1 (MQ=255)                                                                        tACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCg                                                                            <  2:82828/150‑1 (MQ=255)                                                                         aCAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACGCCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGc                                                                           <  1:516866/150‑1 (MQ=255)                                                                         aCAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACGCCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGc                                                                           <  1:511917/150‑1 (MQ=255)                                                                          cAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCt                                                                          >  1:304880/1‑150 (MQ=255)                                                                             cgcgTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGCCCAGCAGCGACGGCGGCGTGGTGCCGCTAAg                                                                       >  2:960938/1‑150 (MQ=255)                                                                             cgcgTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAg                                                                       <  2:606976/150‑1 (MQ=255)                                                                             cgcgTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAg                                                                       <  2:923571/150‑1 (MQ=255)                                                                                gTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTc                                                                    >  2:1146446/1‑150 (MQ=255)                                                                                gTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTc                                                                    <  2:1414626/150‑1 (MQ=255)                                                                                     gTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTg                                                               <  1:143284/150‑1 (MQ=255)                                                                                     gTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTg                                                               <  1:1146446/150‑1 (MQ=255)                                                                                      tCAGCGGCTGATTTCCACTATTTATACTCAGGTCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGc                                                              >  2:1202161/1‑150 (MQ=255)                                                                                       cAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGcc                                                             <  1:368157/150‑1 (MQ=255)                                                                                             cTGATTACCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAAtt                                                       >  1:1083660/1‑150 (MQ=255)                                                                                               gATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGa                                                     >  1:269048/1‑150 (MQ=255)                                                                                                   tCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGagcag                                                 <  2:1037975/150‑1 (MQ=255)                                                                                                      aCTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGt                                              >  2:414310/1‑150 (MQ=255)                                                                                                      aCTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGt                                              >  1:775584/1‑150 (MQ=255)                                                                                                           ttATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAGTACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTgc                                         >  1:1034634/1‑150 (MQ=255)                                                                                                               aCTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATANNANTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGcgccg                                     >  2:31878/1‑150 (MQ=255)                                                                                                               aCTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGcgccg                                     >  2:551506/1‑150 (MQ=255)                                                                                                                cTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGcgcccc                                    <  2:680908/150‑3 (MQ=255)                                                                                                                    ggCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCtctc                                >  1:202099/1‑150 (MQ=255)                                                                                                                      ccAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCa                              <  2:239048/150‑1 (MQ=255)                                                                                                                      ccAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCa                              >  1:209516/1‑150 (MQ=255)                                                                                                                      ccAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCa                              <  2:1354949/150‑1 (MQ=255)                                                                                                                            aGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAAcc                        >  1:725086/1‑150 (MQ=255)                                                                                                                             gTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCa                       <  2:764749/150‑1 (MQ=255)                                                                                                                                 cgcgTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTg                   >  1:1018608/1‑150 (MQ=255)                                                                                                                                  gcgTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTgg                  >  1:775513/1‑150 (MQ=255)                                                                                                                                  gcgTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTgg                  >  1:1370628/1‑150 (MQ=255)                                                                                                                                     tggtgCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTGGATc               <  1:690844/150‑1 (MQ=255)                                                                                                                                            ggAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTGGATCAGTTccc        <  2:1131795/150‑1 (MQ=255)                                                                                                                                             gAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTGGATCAGTTCCCg       <  2:304880/150‑1 (MQ=255)                                                                                                                                                  caacaCCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTGGATCAGTTCCCGGTAAc  >  1:1336604/1‑150 (MQ=255)                                                                                                                                                  aaacaCCGAAAATACCCCAGGCCTCGCGGCGCTGGAAACCATTCGCCTGACCAGCAGCGACGGCGGCGGGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTGGATCAGTTCCCGGTAAc  <  2:407138/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CGTATGTCAATGTTGATCGCGACAGCGCCAGCCGTCTGGGGATCAGCATGGCGGATGTCGATAACGCCCTGTACAACGCGTTTGGTCAGCGGCTGATTTCCACTATTTATACTCAGGCCAACCAGTATCGCGTGGTGCTGGAGCACAACACCGAAAATACCCCAGGCCTCGCGGCGCTGGATACCATTCGCCTGACCAGCAGCGACGGCGGCGTGGTGCCGCTAAGCTCAATTGCCAAAATTGAGCAGCGTTTTGCGCCGCTCTCCATCAACCATCTGGATCAGTTCCCGGTAAC  >  NC_000913/2157415‑2157709

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