Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	654,398	C→A	6.9%	L115L (CTG→CTT)	ctpL ←	methyl‑accepting chemotaxis protein CtpL

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	654,398	0	C	A	6.9%	214.1 / 4.1	67	L115L (CTG→CTT)	ctpL	methyl‑accepting chemotaxis protein CtpL
Reads supporting (aligned to +/- strand):  ref base C (32/30);  new base A (3/2);  total (35/32)
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.99e-01

CCAGATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGG  >  NC_002947/654253‑654545                                                                                                                                                   |                                                                                                                                                     ccAGATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGc                                                                                                                                                   >  3:234493/1‑149 (MQ=255)   aGATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCtt                                                                                                                                                 <  6:267266/149‑1 (MQ=255)               gCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCca                                                                                                                                      <  4:234493/148‑1 (MQ=255)                      gCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑g                                                                                                                             >  2:61471/1‑149 (MQ=255)                        aaGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑Gtt                                                                                                                           <  7:306851/149‑1 (MQ=255)                          gggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCACTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTgg                                                                                                                         <  1:39045/149‑1 (MQ=255)                          gggCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGAAGATGGCCAACTGCTGCAATTACGCCCCCAGAACCCCAACGGATTCCACCAGCA‑GGGATTGTGGTTCGCCGGCCACCTTCCCGCCGCCAAGCAGTGC‑GTAgg                                                                                                                         >  2:116804/1‑149 (MQ=255)                                 aaGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTc                                                                                                                               <  1:287609/137‑1 (MQ=255)                                 aaGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTc                                                                                                                               >  2:287609/1‑137 (MQ=255)                                 aaGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGaa                                                                                                                  <  1:130993/149‑1 (MQ=255)                                  aGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAg                                                                                                                 <  1:61471/149‑1 (MQ=255)                                    tAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGc                                                                                                               <  4:317943/149‑1 (MQ=255)                                          ttGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCGGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAg                                                                                                         <  1:238776/149‑1 (MQ=255)                                              ccTCGGCGTTGCCACTGTCGCGGGCGTCTCTTGCCATCTGCTCGAAGTTCGCACCCAGATCCCGCTCGGCTTGCAGAAGCA‑GGGCTGGTTGGTCGCCGGCACGCTTGCCGGCGGCCCGCAGTTC‑GGTGGCGGTGAAGGCGTGCAgga                                                                                                       >  6:360388/1‑146 (MQ=255)                                               cTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCgg                                                                                                                      <  7:123761/131‑1 (MQ=255)                                               cTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCgg                                                                                                                      >  8:123761/1‑131 (MQ=255)                                                tCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATc                                                                                                   >  1:290611/1‑149 (MQ=255)                                                 cGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGCCGCCGGCAAGCTGGCCGGCGAcc                                                                                                                                       >  2:201427/1‑113 (MQ=255)                                                 cGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGCCGCCGGCAAGCTGGCCGGCGAcc                                                                                                                                       <  1:201427/113‑1 (MQ=255)                                                   gCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGtt                                                                                                                                <  6:107404/118‑1 (MQ=255)                                                         ccACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTggg                                                                                          >  5:30910/1‑149 (MQ=255)                                                           aCTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgc                                                                                        >  3:88926/1‑149 (MQ=255)                                                            cTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGgcg                                                                                       >  1:140202/1‑149 (MQ=255)                                                                     gCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtctc                                                                              <  4:278327/149‑1 (MQ=255)                                                                       gTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCgg                                                                            <  2:140202/149‑1 (MQ=255)                                                                                 aTCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTAc                                                                  >  7:46357/1‑149 (MQ=255)                                                                                 aTCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTAc                                                                  >  3:271666/1‑149 (MQ=255)                                                                                      cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCAGGGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtctc                                                                              >  1:68393/1‑133 (MQ=255)                                                                                      cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgg                                                                                           <  6:159142/119‑1 (MQ=255)                                                                                      cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTgg                                                                                           >  5:159142/1‑119 (MQ=255)                                                                                      cTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTtctc                                                                              <  2:68393/132‑1 (MQ=255)                                                                                       tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCt                                                                                                             <  6:354528/100‑1 (MQ=255)                                                                                       tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCt                                                                                                             >  5:354528/1‑100 (MQ=255)                                                                                       tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                           <  7:40744/134‑1 (MQ=255)                                                                                       tCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGc                                                                           >  8:40744/1‑134 (MQ=255)                                                                                       tCGAAGTTAGCCCCCAGCTCCCCCTCGGCTTGCAGCAGCA‑GGGCGTGTGGGTCGCCGGCCACCTTGCCGGCGGCCAGCAGTTC‑GTTTGGGGTGTAGGGCGGCCGGCTATACACGGTGGGGGGCAACTTATAGGCCCGCTCTACCGGccc                                                            >  2:261068/1‑148 (MQ=255)                                                                                         gAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCa                                                            <  6:30910/147‑1 (MQ=255)                                                                                                gCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAg                                                   >  4:185825/1‑149 (MQ=255)                                                                                                   cccAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGc                                                >  1:62140/1‑149 (MQ=255)                                                                                                           ccGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                              >  2:49020/1‑127 (MQ=255)                                                                                                           ccGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                              <  1:49020/127‑1 (MQ=255)                                                                                                           ccGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTCGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                              >  4:136036/1‑127 (MQ=255)                                                                                                           ccGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTCGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGc                                                              <  3:136036/127‑1 (MQ=255)                                                                                                                 ggCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTcc                                                        >  2:204963/1‑127 (MQ=255)                                                                                                                    ttGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGt                               >  6:327469/1‑149 (MQ=255)                                                                                                                       cagcagca‑gGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGc                                         <  2:126199/136‑1 (MQ=255)                                                                                                                       cagcagca‑gGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGc                                         >  1:126199/1‑136 (MQ=255)                                                                                                                        agcagca‑gGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAgg                                                                                                        >  7:28494/1‑72 (MQ=255)                                                                                                                        agcagca‑gGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAgg                                                                                                        <  8:28494/72‑1 (MQ=255)                                                                                                                         gcagca‑gGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGcc                          >  4:268793/1‑149 (MQ=255)                                                                                                                              a‑gGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                     <  3:185825/149‑1 (MQ=255)                                                                                                                              a‑gGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                     >  8:66833/1‑149 (MQ=255)                                                                                                                                    ttGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGt                                                      >  6:218522/1‑111 (MQ=255)                                                                                                                                    ttGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGt                                                      <  5:218522/111‑1 (MQ=255)                                                                                                                                    ttGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggc                  >  1:151587/1‑147 (MQ=255)                                                                                                                                    ttGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggc                  <  2:151587/147‑1 (MQ=255)                                                                                                                                       tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                     <  7:69518/141‑1 (MQ=255)                                                                                                                                       tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCAt                     >  8:69518/1‑141 (MQ=255)                                                                                                                                       tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggcg              <  7:66833/148‑1 (MQ=255)                                                                                                                                       tgGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGGGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggcg              <  5:360388/148‑1 (MQ=255)                                                                                                                                       ctggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTcc                                                        >  7:191593/6‑106 (MQ=39)                                                                                                                                       ctggccGCCGGCAAGCTGGCCGGCGACCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTcc                                                        <  8:191593/101‑1 (MQ=39)                                                                                                                                        gggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATggcggcggc            <  2:62140/149‑1 (MQ=255)                                                                                                                                         ggTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGcc           >  4:132913/1‑149 (MQ=255)                                                                                                                                          gTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTCGGTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGcc           >  1:47205/1‑149 (MQ=255)                                                                                                                                           tCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGtt                                  >  7:29392/1‑124 (MQ=255)                                                                                                                                           tCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGtt                                  <  8:29392/124‑1 (MQ=255)                                                                                                                                                  ccAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAgg  <  5:327469/149‑1 (MQ=255)                                                                                                                                                   cAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCt                         <  4:196270/125‑1 (MQ=255)                                                                                                                                                   cAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCt                         >  3:196270/1‑125 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CCAGATGCAGGGTGGCGTCTAGCAAGGGCAGCAAGTAGCGATTGGCCTCGGCGTTGCCACTGTCGCGGGCGTAGCTGGCCATCTGCTCGAAGTTCGCCCCCAGCTCCCGCTCGGCTTGCAGCAGCA‑GGGCTTGTGGGTCGCCGGCCAGCTTGCCGGCGGCCAGCAGTTC‑GTTGGCGGTGAAGGCCTGCAGGCTATCCAGGCTGGGGCGCAACTTCTCGGCCAGCTCTACCGGCCAGTCCGTCAGTGCGGCCTGCAACTGTTGGTTCGCCTCCATGGCGGCGGCCTGGCGCAGG  >  NC_002947/654253‑654545

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