Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,591,576	T→A	16.1%	L258M (TTG→ATG) ‡	PP_3169 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,591,576	0	T	A	16.1%	143.1 / 19.6	62	L258M (TTG→ATG) ‡	PP_3169	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (28/24);  new base A (4/6);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.38e-01

TCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGG  >  NC_002947/3591430‑3591717                                                                                                                                                   |                                                                                                                                              tCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTg                                                                                                                                             >  5:406093/1‑149 (MQ=255)      cGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTg                                                                                                                                        <  4:11492/149‑1 (MQ=255)        ggTGCTGGCCGAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                      <  3:337920/149‑1 (MQ=255)        ggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                      <  1:268022/149‑1 (MQ=255)        ggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCAAGCGCTGGt                                                                                                                                      <  2:185531/149‑1 (MQ=255)           gCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCa                                                                                                                                   >  2:46060/1‑149 (MQ=255)               gCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCg                                                                                                                                <  6:406093/148‑1 (MQ=255)                ccAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGAACCAGCGCAAGCGCTGGTGCAACggg                                                                                                                              <  7:173302/149‑3 (MQ=255)                 cAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGGACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGAt                                                                                                                              >  2:390788/1‑148 (MQ=255)                  aaGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGa                                                                                                                            >  6:90302/1‑149 (MQ=255)                  aaGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGAACCAGCGCAAGCGCTGGTGCAACgggtc                                                                                                                            <  6:341285/149‑5 (MQ=255)                       cccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGAACCAGCGCAAGCGCTGGTGCaa                                                                                                                                  <  7:110633/138‑2 (MQ=255)                       cccGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGAACCAGCGCAAGCGCTGGTGCaa                                                                                                                                  >  8:110633/1‑137 (MQ=255)                         cGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTg                                                                                                                     >  2:272051/1‑149 (MQ=255)                             tGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCg                                                                                                                   >  2:188557/1‑147 (MQ=255)                             tGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCg                                                                                                                   <  1:188557/147‑1 (MQ=255)                                 gATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGa                                                                                                             <  1:46060/149‑1 (MQ=255)                                   ttACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGAc                                                                                                            >  5:259113/1‑148 (MQ=255)                                     actacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGTACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCa                                                                                                          <  1:390788/148‑1 (MQ=255)                                     actacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCa                                                                                                          <  3:310676/148‑1 (MQ=255)                                      ctacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGAc                                                                                                            >  1:142057/1‑145 (MQ=255)                                       tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGCCCAGAt                                                                                                       >  8:330875/1‑149 (MQ=255)                                       tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTACACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAt                                                                                                       >  5:161460/1‑149 (MQ=255)                                                aTCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACg                                                                                               >  5:339206/1‑148 (MQ=255)                                                aTCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACg                                                                                               >  7:211110/1‑148 (MQ=255)                                                aTCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACg                                                                                               >  4:291975/1‑148 (MQ=255)                                                aTCAGCAGCATGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACg                                                                                               >  7:69074/1‑148 (MQ=255)                                                   agcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCATCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCAc                                                                                            <  6:259113/148‑1 (MQ=255)                                                   agcagcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCAc                                                                                            <  4:173648/148‑1 (MQ=255)                                                        cGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggtgg                                                                                      <  6:339206/149‑1 (MQ=255)                                                          tGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTc                                                                                    <  3:291975/149‑1 (MQ=255)                                                                 ccAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGAcc                                                                             <  8:69074/149‑1 (MQ=255)                                                                          accCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                    <  1:31942/149‑1 (MQ=255)                                                                             cGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGa                                                                 >  2:236849/1‑149 (MQ=255)                                                                               aaaTGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATc                                                               <  2:235291/149‑1 (MQ=255)                                                                                    cAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                          <  5:90302/149‑1 (MQ=255)                                                                                    cAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                          <  1:272051/149‑1 (MQ=255)                                                                                    cAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                          <  6:93229/149‑1 (MQ=255)                                                                                     aaCGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                          >  7:76904/1‑148 (MQ=255)                                                                                     aaCGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                          >  3:123015/1‑148 (MQ=255)                                                                                                 gcAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCttt                                             >  4:234997/1‑149 (MQ=255)                                                                                                  cAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCACATGGGGATCCCGGGGAATTACACCTTTc                                            <  1:236849/149‑1 (MQ=255)                                                                                                     ccTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgcg                                                                         >  7:75752/1‑117 (MQ=255)                                                                                                     ccTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgcg                                                                         <  8:75752/117‑1 (MQ=255)                                                                                                      cTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCtt                                        >  3:235232/1‑149 (MQ=255)                                                                                                         gAACTGGCGCAAGCACTGAACGATGACCCGTTACACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTc                                       <  6:161460/147‑1 (MQ=255)                                                                                                               gcgcAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGc                               >  4:234653/1‑149 (MQ=255)                                                                                                                 gcAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGCTCCCGGGGAATTACACCTTTCGCTTCAAGATGGCTc                             >  6:307331/1‑149 (MQ=255)                                                                                                                   aaGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAg                           >  1:272625/1‑149 (MQ=255)                                                                                                                   aaGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTCCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGCGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCac                           >  6:286783/1‑148 (MQ=255)                                                                                                                    aGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAgc                          <  8:76904/149‑1 (MQ=255)                                                                                                                     gCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAgcg                         <  4:123015/149‑1 (MQ=255)                                                                                                                       aCTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAgcgc                        >  2:275977/1‑148 (MQ=255)                                                                                                                        cTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGa                      <  4:235232/149‑1 (MQ=255)                                                                                                                           aaCGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGatc                    <  3:234997/148‑1 (MQ=255)                                                                                                                                   cccGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAggg            >  5:7206/1‑148 (MQ=255)                                                                                                                                   cccGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAggg            <  3:234653/148‑1 (MQ=255)                                                                                                                                   cccGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAggg            >  1:195540/1‑148 (MQ=255)                                                                                                                                       gtgaaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTg       >  4:168944/5‑149 (MQ=255)                                                                                                                                         gaaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGa     >  7:112559/3‑149 (MQ=255)                                                                                                                                           aCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAATGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGTGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATg   >  1:300137/1‑149 (MQ=255)                                                                                                                                            ccAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATgg  >  1:167307/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                              TCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTTGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGG  >  NC_002947/3591430‑3591717

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