Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,837	T→C	S35S (AGT→AGC)	slt →	soluble lytic murein transglycosylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,630,837	0	T	C	95.3%	112.1 / ‑2.1	44	S35S (AGT→AGC)	slt	soluble lytic murein transglycosylase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (20/21);  total (22/22)
Fisher's exact test for biased strand distribution p-value = 4.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.00e-01

CCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGC  >  NC_000913/4630588‑4631083                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       ccGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAGGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTAATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTc                                                                                                                                                                                                                                                       >  1:131987/1‑251 (MQ=255)        aTCCCAGTTGGTCGATCACCTTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTAAGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTAcgc                                                                                                                                                                                                                                                >  1:342701/1‑251 (MQ=255)                             tCGCAGTAAAAGTAAAAAAGTGTCCTTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTggg                                                                                                                                                                                                                            <  1:237142/250‑1 (MQ=255)                                  gtaaaagtaaaaAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATc                                                                                                                                                                                                                      <  1:297168/251‑1 (MQ=255)                                   taaaagtaaaaAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCg                                                                                                                                                                                                                     >  1:372683/1‑251 (MQ=255)                                                 tgtCCGCAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCCTTGATGTATTTACACTTAGGGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGAGAATCGACAAATGGCTGTgg                                                                                                                                                                                                       >  2:146205/1‑251 (MQ=255)                                                  gtCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGt                                                                                                                                                                                                      <  1:13306/251‑1 (MQ=255)                                                    cctgaaCGTGGCGGAATCGGCAATTATTGGTTTGCATAAATCTTTTCCGCGGCATTATGCGTCATTGATGTATTTACACTTAGAGGATGCCCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGGGTCTGTCTGCTGACGGTCAGCAGCTTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCg                                                                                                                                                                                                    <  2:118417/247‑1 (MQ=255)                                                               cGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCCTGACGCTGCATTGATGTATTTACACTTAGAGGATTCTCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCTCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAatgatg                                                                                                                                                                                         <  1:256989/251‑1 (MQ=255)                                                                 tAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGcc                                                                                                                                                                                       >  1:156863/1‑251 (MQ=255)                                                                  aaaCGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCt                                                                                                                                                                                      <  1:149516/251‑1 (MQ=255)                                                                          aTGACTGGTTAGCATAAAGCATTTACGGGGCATGACGTGGCATTGATTTTGTTCCACTTGGAGGATGCTCTGGTTGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGTCTGCCGGTGTCTGTCGTCAGACGGTCAGCAGCGTGGCGCGAGCCGATTTACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATTATGCCTGGACTGaa                                                                                                                                                                              <  2:117077/251‑1 (MQ=255)                                                                            gACTGGTTAGCATAACTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGaa                                                                                                                                                                              >  1:82991/1‑249 (MQ=255)                                                                                 ggtagcactaaTCTAATGCGCGGCATTACGCTGCATTGATGTATTTACACTTAGAGTATGCGCCTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTTCTGACTGTAAGCGGCGTGGCCCTAGCCGACTCACTGGATGAGCAGCGTATCCGGTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAAGGATGCCTGGACTGaaatatta                                                                                                                                                                        <  2:16872/241‑7 (MQ=255)                                                                                       aaaaatttaatacgcggatagACGCTGCATGGATGTATTTCAACTTAGAGGATGCGCTTGGGGAAAAAGCCAAAAAAGTTCCCTGTCGGGTGGTGGCTGCCGGTGTCTTTCTGCTGTCGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGACGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGGGGGCGAACAAATGATGCCTGGACTGAAAGAGTATCCACtt                                                                                                                                                                 <  2:156863/231‑1 (MQ=255)                                                                                          aaTCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTAt                                                                                                                                                              >  1:18737/1‑251 (MQ=255)                                                                                            tCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGATGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGGCAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATcc                                                                                                                                                            >  1:362976/1‑251 (MQ=255)                                                                                                   cgcgGCATGACGCTGCATTGATGTATTTACACTTAGAGGATTCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTg                                                                                                                                                     >  1:152035/1‑251 (MQ=255)                                                                                                     cgGCATGACGCTTCATTGATGTATTTACACGTAGAGGAGGCGCTTGTGGAAAAAGCCAAACAAGTTACCGGGCGGCTGTTGGCTGCCGGTGTCTGTCGTCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGa                                                                                                                                                   <  1:383609/251‑1 (MQ=255)                                                                                                          tGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAATCCAAACAAGTTACCTGACGGCTGTTGGCTGCCGGTGTCAGTCTGCTTAAGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAACAGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGGTGTGGTCGAACAACTGATGCCTGGAATGAAAGATTATCCACTTTATCCCTAACTGGAATAcc                                                                                                                                              >  2:287179/1‑251 (MQ=255)                                                                                                                      gATGTATTTACACTTAGAGGATGCGCTTGTTGAAAAAGCCAAACAAGTTACCAGGCGGCTGTTGGCTGCCGGTGTCTGGCTGCTGACGGTCAGCAGCTTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCg                                                                                                                                  <  2:61042/251‑1 (MQ=255)                                                                                                                         gTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAATTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGAGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAAACAAGCAGGCCTGGGATAATCGACAAATGGCTGCGGTCGCACAAATGATGCCTGGACTGAAGGATTATCCGCTTTACCCCTACCTGGAATACCGCCAGATCACCgatg                                                                                                                               >  2:166629/1‑251 (MQ=255)                                                                                                                         gTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCgat                                                                                                                                <  1:303670/250‑1 (MQ=255)                                                                                                                              taaacTTTGAGTATGCGCTGGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGTTGTCTGTTTGCTGACGGTCAGCAGCGGGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCgatgacct                                                                                                                           <  2:338386/247‑4 (MQ=255)                                                                                                                                aaattAGATGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCtgat                                                                                                                        <  1:257291/248‑1 (MQ=255)                                                                                                                                     agagGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATc                                                                                                                   >  1:76005/1‑251 (MQ=255)                                                                                                                                             cgcTTGTGGAAAACGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACcggcg                                                                                                           <  1:232798/251‑1 (MQ=255)                                                                                                                                                 tggggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCTACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTCCCTGGAATACCGTCAGAGCACCGATGACTTGATGAATCAACCGGCGGt                                                                                                         <  1:41860/246‑1 (MQ=255)                                                                                                                                                    ggAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGActg                                                                                                    >  1:93071/1‑249 (MQ=255)                                                                                                                                                             cAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGAAGCGTAGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGATACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGCGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACATGATGAATCAACCGGCGGTGACTGTCTCTAACt                                                                                           >  1:182575/1‑251 (MQ=255)                                                                                                                                                                  aaGTTTCCTGGCGGCTGTTGTCGGCCGGTGTCTGGCTGCTGACGGTCAGCAGCGTGGCGCGAGCCTACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACGTTATCCCTACCTGGAATACCGTCAGATCACCGTTGGCCTGATGAATCAACCGGCGGTGATTGTCACGAACTTTGtt                                                                                      <  2:192554/251‑1 (MQ=255)                                                                                                                                                                                 tGTTGGCTGCCGGTGACTGTCTGCTGACGGTCAGAAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGGGACTGTCACTAACTTTGTTCGCGCTAAACCCACg                                                                       >  1:109147/1‑251 (MQ=255)                                                                                                                                                                                        tGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTtcctc                                                                >  1:291028/1‑251 (MQ=255)                                                                                                                                                                                                tctgtctgCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGc                                                        >  1:33422/1‑251 (MQ=255)                                                                                                                                                                                                          gACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCa                                                <  1:344879/249‑1 (MQ=255)                                                                                                                                                                                                               tCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCACTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGGCTGTCACTAACTTTGTTCGCGCAAACCCCACGCTTCCTCCCACTCGCACGCTGCAATCTCGt                                         >  1:114909/1‑251 (MQ=255)                                                                                                                                                                                                                              gAGCCGACGCACTGGATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACtt                          <  1:235210/251‑2 (MQ=255)                                                                                                                                                                                                                                           ggATGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                           >  1:148927/1‑221 (MQ=255)                                                                                                                                                                                                                                             aTGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGAc           <  1:146205/251‑1 (MQ=255)                                                                                                                                                                                                                                              tGAGCAGCGTAGCCGTTCCGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCa                                                                                                                  >  2:171816/1‑147 (MQ=255)                                                                                                                                                                                                                                              tGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCa                                                                                                                  <  1:171816/147‑1 (MQ=255)                                                                                                                                                                                                                                              tGAGCAGCGTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCACTTTATCCCTACCTGGAATACCGTCAGATCACCGATGACCTGATGAATCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                           >  1:347843/1‑218 (MQ=255)                                                                                                                                                                                                                                              ggAGCAGCGTAGCCGTGACGCGCAAATCAAGCAGGCCTGGTATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCCTCTTTATCCCTACCTGGAATACCGTCAGATCTCCGATGACCTGATGAAGCAACCGGCGGTGACTGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAAtctc                                           <  2:347843/217‑1 (MQ=255)                                                                                                                                                                                                                                                      gTAGCCGTTACGCGCAAATCAAGCAGGCCTGGGTTAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAAGATTATCAACTTTATCCCTACCTGGAATACCATCAGATCACAGATGACCTGATGAATCAACCGGCGGTGACTGTCCCTAACTTTTTTCGCGCTAAACCCACGCTTCCCCACACCCGCCCTCCGAAATCCCCTTTCGTCACTGAACTGGCGCGGCGTGAAGaacagcgtggg  >  2:160982/1‑241 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       CCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGTAAAAGTAAAAAAGTGTCCGTAACGTGGCGTAAACGGCAATGACTGGTTAGCATAAATCTATTACGCGGCATGACGCTGCATTGATGTATTTACACTTAGAGGATGCGCTTGTGGAAAAAGCCAAACAAGTTACCTGGCGGCTGTTGGCTGCCGGTGTCTGTCTGCTGACGGTCAGCAGCGTGGCGCGAGCCGACTCACTGGATGAGCAGCGTAGTCGTTACGCGCAAATCAAGCAGGCCTGGGATAATCGACAAATGGATGTGGTCGAACAAATGATGCCTGGACTGAAGGATTATCCGCTTTATCCCTACCTGGAATACCGCCAGATCACCGATGATCTGATGAATCAACCGGCGGTGACGGTCACTAACTTTGTTCGCGCTAACCCCACGCTTCCTCCCGCTCGCACGCTGCAATCTCGTTTCGTCAATGAACTGGCGCGGCGTGAAGACTGGCGTGGC  >  NC_000913/4630588‑4631083

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