Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,984,042	A→G	9.7%	intergenic (+36/+154)	ilvE → / ← PP_3512	branched‑chain‑amino acid aminotransferase/membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,984,042	0	A	G	9.7%	167.1 / 8.0	62	intergenic (+36/+154)	ilvE/PP_3512	branched‑chain‑amino acid aminotransferase/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/25);  new base G (3/3);  total (34/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.46e-01

ACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGG  >  NC_002947/3983897‑3984186                                                                                                                                                    |                                                                                                                                                  aCGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATg                                                                                                                                                <  1:2033/149‑1 (MQ=255)    aCCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCa                                                                                                                                             <  4:233287/149‑1 (MQ=255)        ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCATTGCTGGCATGc                                                                                                                                          <  1:151738/148‑1 (MQ=255)        ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCATTGCTGGCATGc                                                                                                                                          >  2:151738/1‑148 (MQ=255)                  ggCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTc                                                                                                                               <  1:184960/149‑1 (MQ=255)                    cccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCATTGCTGGCATGCCCCGtt                                                                                                                                    <  1:17431/142‑1 (MQ=255)                    cccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCATTGCTGGCATGCCCCGtt                                                                                                                                    >  2:17431/1‑142 (MQ=255)                         ttACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGtt                                                                                                                        <  2:116224/149‑1 (MQ=255)                         ttACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGtt                                                                                                                        <  2:112443/149‑1 (MQ=255)                           aCCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGttt                                                                                                                       >  6:88615/1‑148 (MQ=255)                           aCCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGttt                                                                                                                       >  8:13662/1‑148 (MQ=255)                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTaca                                                                                                                    >  4:52734/1‑149 (MQ=255)                                   gCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGccc                                                                                                              >  3:76084/1‑149 (MQ=255)                                   gCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGccc                                                                                                              >  3:71702/1‑149 (MQ=255)                                          aaCGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                         >  5:170378/1‑147 (MQ=255)                                          aaCGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                         >  4:80182/1‑147 (MQ=255)                                                  gACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCggg                                                                                               >  5:306850/1‑149 (MQ=255)                                                   aCCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGt                                                                                              <  4:76084/149‑1 (MQ=255)                                                    ccGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGt                                                                                              <  4:71702/148‑1 (MQ=255)                                                     cGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGtt                                                                                             >  7:265834/1‑148 (MQ=255)                                                          tCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTAccc                                                                                        <  6:306850/148‑1 (MQ=255)                                                              gagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                   <  8:77560/149‑1 (MQ=255)                                                                gcggcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGATCGGGTTTACCCGCGAATa                                                                                 <  3:80182/149‑1 (MQ=255)                                                                  atctaTGTAGAGGCACGGGCAGGCTGGATGGTCATGGTCGCTTGAAACGAATCAGCTTTATT‑‑AAAAGGGGGAATGCCAGCAATGGATGGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAg                                                                                <  6:218376/144‑1 (MQ=255)                                                                   gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAgtc                                                                              >  2:217630/1‑147 (MQ=255)                                                                     gATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAg                                                                                >  5:218376/1‑145 (MQ=255)                                                                     gATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCca                                                                            >  2:10180/1‑149 (MQ=255)                                                                       tGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCcaca                                                                          <  6:146713/149‑1 (MQ=255)                                                                          agagGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCacacac                                                                       >  8:129818/1‑147 (MQ=255)                                                                              gCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgc                                                                     >  4:286414/1‑147 (MQ=255)                                                                               cACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgccag                                                                  <  3:286414/149‑1 (MQ=255)                                                                                aCCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgccagc                                                                 >  2:264806/1‑149 (MQ=255)                                                                                  cGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgccagc                                                                 <  5:94399/147‑1 (MQ=255)                                                                                  cGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgccagc                                                                 >  6:94399/1‑147 (MQ=255)                                                                                        gCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAAAACGGGGCATGCCATTGCTGGCATGcc                                                                                                                                         >  5:61953/1‑71 (MQ=17)                                                                                        gCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAAAACGGGGCATGCCATTGCTGGCATGcc                                                                                                                                         <  6:61953/71‑1 (MQ=17)                                                                                        gCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATg                                                         >  1:226843/1‑149 (MQ=255)                                                                                        gCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATg                                                         >  4:159036/1‑149 (MQ=255)                                                                                          tGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑cc                                                      >  4:309449/1‑149 (MQ=255)                                                                                            gATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑cc                                                      <  3:159036/147‑1 (MQ=255)                                                                                                 tCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACTGCCCCCTCGTAGGAGCgg                                                                                                <  6:218966/101‑1 (MQ=255)                                                                                                 tCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCgg                                                                                                >  5:218966/1‑101 (MQ=255)                                                                                                 tCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTgg                                                   >  6:41900/1‑145 (MQ=255)                                                                                                      gTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCAcc                                          >  3:126398/1‑149 (MQ=255)                                                                                                       tCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCg                                         >  2:309411/1‑149 (MQ=255)                                                                                                         gCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGtt                                       <  3:52734/149‑1 (MQ=255)                                                                                                         gCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGtt                                       <  3:232296/149‑1 (MQ=255)                                                                                                         gCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGtt                                       >  4:232296/1‑149 (MQ=255)                                                                                                            tGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGt                                     >  7:147463/1‑148 (MQ=255)                                                                                                                aCGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaat                                                             >  8:146570/1‑121 (MQ=255)                                                                                                                aCGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaat                                                             <  7:146570/121‑1 (MQ=255)                                                                                                                aCGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTg                                <  8:265834/149‑1 (MQ=255)                                                                                                                         cTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCa                        >  8:20139/1‑148 (MQ=255)                                                                                                                              tG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCg                  >  6:17695/1‑149 (MQ=255)                                                                                                                               g‑‑AAAACGGGGCATGCCAGGAATGGCATGCCCCGTTTTTTTATTTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGGCAATAATGCCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGc                   <  1:217630/148‑1 (MQ=255)                                                                                                                                  aaaaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCggg                <  5:17695/149‑1 (MQ=255)                                                                                                                                    aaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGc               >  1:114273/1‑148 (MQ=255)                                                                                                                                     aCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCaa             <  6:170378/149‑1 (MQ=255)                                                                                                                                       ggggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGc           <  3:309449/149‑1 (MQ=255)                                                                                                                                        gggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGcc          >  8:136515/1‑149 (MQ=255)                                                                                                                                            aTGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAAAG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcgcg      >  5:17075/1‑149 (MQ=255)                                                                                                                                                cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg  <  7:136515/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  ACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATG‑‑AAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATG‑CCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGG  >  NC_002947/3983897‑3984186

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