Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,767,352	T→G	5.3%	M1174L (ATG→CTG)	PP_5060 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,767,352	0	T	G	5.3%	212.6 / 2.0	69	M1174L (ATG→CTG)	PP_5060	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base G (2/2);  total (35/34)
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.14e-01

GGCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGA  >  NC_002947/5767206‑5767495                                                                                                                                                   |                                                                                                                                                 ggCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATgg                                                                                                                                                <  4:15697/149‑1 (MQ=255)    gACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCt                                                                                                                                             >  1:211875/1‑149 (MQ=255)    gACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCt                                                                                                                                             <  7:199629/149‑1 (MQ=255)       aGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTgg                                                                                                                                           >  2:51311/1‑148 (MQ=255)        ggCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggtc                                                                                                                                         <  2:353433/149‑3 (MQ=255)         gCAAGTCACAGAACATCAGCAGGTGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGaa                                                                                                                                         >  7:56523/1‑148 (MQ=255)         gCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGaa                                                                                                                                         >  4:79451/1‑148 (MQ=255)         gCAAGTCACAGAACATCAGCAGATGAAGCGATTGTGCGGCCACCCGCACACGCTGACTGGCGAGCGCATGGCGGAAGATCATGGAGTGCTGGGTAGGAGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCAAGGCCTCGaa                                                                                                                                         >  4:59179/1‑148 (MQ=255)           aaGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGaaaac                                                                                                                                      >  4:196536/1‑148 (MQ=255)                acaGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑a                                                                                                                                >  7:342229/1‑149 (MQ=255)                acaGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑a                                                                                                                                >  1:75350/1‑149 (MQ=255)                  aGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATg                                                                                                                              >  5:113245/1‑149 (MQ=255)                    aaCATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGtt                                                                                                                            <  5:377859/149‑1 (MQ=255)                       aTCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGAt                                                                                                                         >  7:55743/1‑149 (MQ=255)                            cagATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGttg                                                                                                                    <  7:249556/149‑3 (MQ=255)                               atgatgCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGgtcg                                                                                                                 >  7:9387/1‑149 (MQ=255)                                     cGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTAc                                                                                                           <  8:51271/149‑1 (MQ=255)                                               gCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGa                                                                                                  >  7:234493/1‑148 (MQ=255)                                                  aCCCGCAGAGGCTGACTGTCGAGAGCATGGCGGAAGATGATGTAGTGCTGGGTAGGCGTTCCTTTGTGCCACTGCAGGGCCAGCTCTACCAGGCCCATTGCCTGGATAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAAt                                                                                              <  3:241017/149‑1 (MQ=255)                                                  aCCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAAt                                                                                              <  8:9387/149‑1 (MQ=255)                                                  aCCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAAt                                                                                              <  3:79451/149‑1 (MQ=255)                                                      gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTTTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTc                                                                                          <  6:122321/149‑1 (MQ=255)                                                      gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTc                                                                                          <  4:230377/149‑1 (MQ=255)                                                      gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTc                                                                                          <  8:342229/149‑1 (MQ=255)                                                      gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggtcga                                                                                                                                       <  1:51816/105‑5 (MQ=255)                                                      gCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCATGGGCCTggtcga                                                                                                                                       >  2:51816/1‑101 (MQ=255)                                                        aGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGCGTAGGCGTGCCTTTGTGCCTCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCga                                                                                        <  3:59179/149‑1 (MQ=255)                                                             cTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACgg                                                                                   <  3:196536/149‑1 (MQ=255)                                                               gACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGt                                                                                  >  4:85898/1‑148 (MQ=255)                                                                aCTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGt                                                                                  <  6:363151/147‑1 (MQ=255)                                                                  tGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCa                                                                                >  6:342423/1‑147 (MQ=255)                                                                  tGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCa                                                                                >  6:118170/1‑147 (MQ=255)                                                                    gCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGt                                                                            >  6:307528/1‑149 (MQ=255)                                                                       aGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGt                                                                         <  2:211875/149‑1 (MQ=255)                                                                          gcATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCGCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTAtt                                                                      >  7:30920/1‑149 (MQ=255)                                                                          gcATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCGCTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTAtt                                                                      >  7:36447/1‑149 (MQ=255)                                                                            aTGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGa                                                                    >  5:65207/1‑149 (MQ=255)                                                                            aTGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGa                                                                    >  5:98641/1‑149 (MQ=255)                                                                                       tCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACt                                                         >  3:42969/1‑149 (MQ=255)                                                                                         aTGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTg                                                       <  1:338184/149‑1 (MQ=255)                                                                                          tGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGc                                                      <  3:85898/149‑1 (MQ=255)                                                                                                  tGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgat                                              >  7:101162/1‑149 (MQ=255)                                                                                                        ggCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCgg                                         <  2:75350/148‑1 (MQ=255)                                                                                                        ggCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGt                                        <  6:113245/149‑1 (MQ=255)                                                                                                           gTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGAt                                     >  5:139923/1‑149 (MQ=255)                                                                                                           gTGCCTTTGAGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGAt                                     >  1:192552/1‑149 (MQ=255)                                                                                                            tGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATg                                    <  6:98641/149‑1 (MQ=255)                                                                                                               cTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgatgtcg                                          <  6:170020/140‑1 (MQ=255)                                                                                                               cTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgtcgatgtcg                                          >  5:170020/1‑140 (MQ=255)                                                                                                               cTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGc                                  >  4:282324/1‑148 (MQ=255)                                                                                                                    tgCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAg                            >  7:323177/1‑149 (MQ=255)                                                                                                                       cACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGc                          >  8:198168/1‑148 (MQ=255)                                                                                                                               ggCCAGCGCGACCAGGCCCAGGGCCGGGAAAACACGG‑AGGTTGAGCGTGGTCGGGTTACGGCTCGCGATAATGCTCGAGCCGGGCAGAGTGGTATTGAGGTATGGCACTTGCATGTCGATGTCGGGGATGGCCACCCGGGGTGCCAgct                 >  8:150706/1‑149 (MQ=255)                                                                                                                                gCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGGAATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATg                                    >  1:60419/1‑130 (MQ=255)                                                                                                                                gCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATg                                    <  2:60419/129‑1 (MQ=255)                                                                                                                                   aGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgt                                                  >  6:301877/1‑112 (MQ=255)                                                                                                                                   aGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCatgt                                                  <  5:301877/112‑1 (MQ=255)                                                                                                                                   aGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAgctggct             <  4:330896/149‑1 (MQ=255)                                                                                                                                    gCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAgctggctg            <  8:56523/149‑1 (MQ=255)                                                                                                                                     cTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGctggctgg           <  6:139923/149‑1 (MQ=255)                                                                                                                                     cTCGACCAGGCCCATGGCCAGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCAAACAGGCGATCCAgctggctg            >  2:170428/1‑148 (MQ=255)                                                                                                                                        gACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGgag         <  6:65207/148‑1 (MQ=255)                                                                                                                                          cgaGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACa      <  4:42969/147‑1 (MQ=255)                                                                                                                                           cAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGCCGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACAt     <  7:150706/149‑1 (MQ=255)                                                                                                                                             ggCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACAtt    >  6:158106/1‑148 (MQ=255)                                                                                                                                             ggCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtg   <  7:198168/149‑1 (MQ=255)                                                                                                                                              gCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtga  >  7:286972/1‑149 (MQ=255)                                                                                                                                              gCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtga  <  3:282324/149‑1 (MQ=255)                                                                                                                                              gCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATtga  <  5:342423/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GGCGACAGGCAAGTCACAGAACATCAGCAGATGATGCGTTTGTGCGGCCACCCGCAGACGCTGACTGGCGAGCGCATGGCGGAAGATCATGTAGTGCTGGGTAGGCGTGCCTTTGTGCCACTGCAGGGCCAGCTCGACCAGGCCCATGGCCTGGAAAAACCGG‑ATGTTGATCGTGGTCGTCTTACGGCTCGCGATAATGCTCGAGACGGTCAGAGTGGTATTGACGTAAGGCACTTGCATGTCGATGTCGGTGATGGCCAACAGGCGATCCAGCTGGCTGGAGACATTGA  >  NC_002947/5767206‑5767495

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