Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,646,746	C→T	5.8%	A194T (GCA→ACA)	pbpC ←	penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,646,746	0	C	T	5.8%	224.4 / 2.2	68	A194T (GCA→ACA)	pbpC	penicillin‑insensitive murein repair transglycosylase; inactive transpeptidase domain protein
Reads supporting (aligned to +/- strand):  ref base C (29/35);  new base T (2/2);  total (31/37)
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.98e-01

CCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGCCATTCCAGTT  >  NC_000913/2646603‑2646889                                                                                                                                                |                                                                                                                                                ccTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCt                                                                                                                                            >  1:502895/1‑149 (MQ=255)  cTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGGAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTc                                                                                                                                           >  4:661135/1‑149 (MQ=255)  cTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTc                                                                                                                                           >  4:98835/1‑149 (MQ=255)    tGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCgg                                                                                                                                         >  4:760167/1‑149 (MQ=255)        tctTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATa                                                                                                                                     >  4:627476/1‑149 (MQ=255)          ttttACCTGCTCACGGGACCACACACCTTGCACTGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGc                                                                                                                                   >  4:579025/1‑149 (MQ=255)          ttttACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGc                                                                                                                                   >  1:231352/1‑149 (MQ=255)              aCCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTaa                                                                                                                               >  1:747322/1‑149 (MQ=255)               ccTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTaaa                                                                                                                              >  3:797828/1‑149 (MQ=255)                  gCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTc                                                                                                                           >  1:456671/1‑149 (MQ=255)                     cACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCa                                                                                                                        <  2:747322/149‑1 (MQ=255)                        gggACCACACACCTTGCACGGCCATCAGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGc                                                                                                                     >  1:94483/1‑149 (MQ=255)                         ggACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCg                                                                                                                    <  2:758321/149‑1 (MQ=255)                             cacacaCCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGAttt                                                                                                                <  2:239314/149‑1 (MQ=255)                              acacacCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGAtttt                                                                                                               >  4:697587/1‑149 (MQ=255)                               cacacCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTc                                                                                                              <  3:162787/149‑1 (MQ=255)                               cacacCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTc                                                                                                              <  2:699851/149‑1 (MQ=255)                                  acCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCga                                                                                                           <  4:701585/149‑1 (MQ=255)                                    cTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCgaga                                                                                                         <  2:358442/149‑1 (MQ=255)                                       gCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATaa                                                                                                      <  4:797828/149‑1 (MQ=255)                                         aCGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGc                                                                                                    >  1:293315/1‑149 (MQ=255)                                             ccATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCaa                                                                                                <  1:750309/149‑1 (MQ=255)                                                tCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACtt                                                                                             >  4:801613/1‑149 (MQ=255)                                                tCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACtt                                                                                             >  4:842707/1‑149 (MQ=255)                                                         gCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCg                                                                                    >  3:445012/1‑149 (MQ=255)                                                             tttATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATcc                                                                                >  2:468764/1‑149 (MQ=255)                                                             tttATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATcc                                                                                <  3:272966/149‑1 (MQ=255)                                                             tttATTACGCGCGGCTTCGGCACGCTCCGGCCAAAGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCTATTTTCCGAGATAAGCCCAACTTGCCGCATCGATcc                                                                                >  3:577626/1‑149 (MQ=255)                                                                attaCGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCACCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCt                                                                             <  3:733936/149‑1 (MQ=255)                                                                   aCGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCa                                                                          >  4:316304/1‑149 (MQ=255)                                                                     gcgcgGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGTCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAAc                                                                        >  1:564474/1‑149 (MQ=255)                                                                     gcgcgGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAAc                                                                        >  4:177495/1‑149 (MQ=255)                                                                     gcgcgGCTTCGGCACGCCCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGTCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAAc                                                                        >  1:207554/1‑149 (MQ=255)                                                                          gCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTAcc                                                                   <  3:801613/149‑1 (MQ=255)                                                                           cTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCg                                                                  <  3:697587/149‑1 (MQ=255)                                                                            ttCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGc                                                                 >  1:167089/1‑149 (MQ=255)                                                                                    gCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGa                                                         <  3:98835/149‑1 (MQ=255)                                                                                     cTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAg                                                        <  2:167089/149‑1 (MQ=255)                                                                                      tCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAgc                                                       <  3:835644/149‑1 (MQ=255)                                                                                             aaCGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTa                                                >  1:863398/1‑149 (MQ=255)                                                                                              aCGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTaa                                               >  2:618245/1‑149 (MQ=255)                                                                                              aCGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTaa                                               <  3:591848/149‑1 (MQ=255)                                                                                               cGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAg                                              <  4:577626/149‑1 (MQ=255)                                                                                                  tCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATa                                           >  3:487065/1‑149 (MQ=255)                                                                                                  tCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATa                                           >  1:58610/1‑149 (MQ=255)                                                                                                     gggCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACaa                                        >  1:679326/1‑149 (MQ=255)                                                                                                     gggCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACaa                                        <  3:760167/149‑1 (MQ=255)                                                                                                      ggCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAg                                       >  4:446690/1‑149 (MQ=255)                                                                                                         ggaGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTc                                    <  3:661135/147‑1 (MQ=255)                                                                                                          aaGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCa                                   <  3:297218/149‑1 (MQ=255)                                                                                                             ccGGCTGGGCGCTTGCGGCAAAACCGCCAGCATAGCCGCCTCGGAATAGCTTAAATGCGCAGGCGATTTTCCGAAATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGaa                                <  4:720235/148‑1 (MQ=255)                                                                                                                 cTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTc                            <  3:400055/149‑1 (MQ=255)                                                                                                                   gggCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCAc                          <  2:863398/149‑1 (MQ=255)                                                                                                                   gggCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCAc                          >  2:253428/1‑149 (MQ=255)                                                                                                                       gcTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAGCGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCtt                      <  2:293315/149‑1 (MQ=255)                                                                                                                        cTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTa                     <  1:2524/149‑1 (MQ=255)                                                                                                                              gCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGa               >  2:875910/1‑149 (MQ=255)                                                                                                                               cAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGAt              <  2:502895/149‑1 (MQ=255)                                                                                                                               cAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGAt              <  3:316304/149‑1 (MQ=255)                                                                                                                                 aaaCCGCCAGCATTGTCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGc            <  2:564474/149‑1 (MQ=255)                                                                                                                                 aaaCCGCCAGCATTGTCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGc            <  2:207554/149‑1 (MQ=255)                                                                                                                                  aaCCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGcc           <  2:90611/149‑1 (MQ=255)                                                                                                                                      gCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGCCATTc       <  4:485655/149‑1 (MQ=255)                                                                                                                                      gCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCACCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGCCATTc       <  3:107077/149‑1 (MQ=255)                                                                                                                                        cAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGCCATTCCa     <  2:537455/149‑1 (MQ=255)                                                                                                                                           cATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGCCATTCCAGtt  <  1:721987/149‑1 (MQ=255)                                                                                                                                           cATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGCCATTCCAGtt  <  2:456671/149‑1 (MQ=255)                                                                                                                                           cATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGCCATTCCAGtt  <  2:58610/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                CCTTGACTCTTTTACCTGCTCACGGGACCACACACCTTGCACGGCCATCCGTTCGAGCACTTTATTACGCGCGGCTTCGGCACGCTCCGGCCAACGATCCGGGCGAAGACGGCTGGGCGCTTGCGGCAAAACCGCCAGCATTGCCGCCTCGGAATAGCTTAAATTCGCAGGCGATTTTCCGAGATAAGCCCAACTTGCCGCACCGATCCCCTGCAACGTACCGCCAAACGGAGCGCGGTTAAGATACAAGGTCAGAATTTCACGCTTAGACAGATGCCATTCCAGTT  >  NC_000913/2646603‑2646889

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑