Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,226,707	T→A	6.6%	intergenic (‑217/‑114)	gltI ← / → PP_1072	glutamate/aspartate ABC transporter substrate‑binding protein/leucine‑rich repeat‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,226,707	0	T	A	6.6%	170.2 / 2.3	61	intergenic (‑217/‑114)	gltI/PP_1072	glutamate/aspartate ABC transporter substrate‑binding protein/leucine‑rich repeat‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (31/25);  new base A (2/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.85e-01

TGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTCCATC  >  NC_002947/1226567‑1226853                                                                                                                                               |                                                                                                                                                   tGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGaa                                                                                                                                            <  6:173224/149‑1 (MQ=255)    aTTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACg                                                                                                                                          >  1:274527/1‑148 (MQ=255)    aTTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGATAACCGAACga                                                                                                                                         <  5:46607/149‑2 (MQ=255)             ccAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTCTCGACCTTAAACCACTGGTATTGCTCGTGATTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCAAGGCCGAATAAACAGAGATTGCTCGTACGGTTATCCTAACGCATTGGCTGc                                                                                                                                >  5:223496/1‑149 (MQ=255)              cAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGc                                                                                                                                >  6:35119/1‑148 (MQ=255)              cAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGATAACCGAACga                                                                                                                                         >  1:151850/1‑138 (MQ=255)              cAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGATAACCGAACga                                                                                                                                         <  2:151850/139‑2 (MQ=255)               aaGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGGTATCCGAACGCATGGGCTGCCa                                                                                                                              >  4:78149/1‑149 (MQ=255)                   gtAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAAcc                                                                                                                           >  3:262476/1‑148 (MQ=255)                    tAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTa                                                                                                                         >  3:228942/1‑149 (MQ=255)                            gACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTa                                                                                                                 >  7:190175/1‑149 (MQ=255)                             aCTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGACTCCAGACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCgg                                                                                                                   >  1:97959/1‑148 (MQ=255)                                 gTGCCATGCCCTGCTGTCGACCTTATACCACTGGTATTGCTCGTGCTTGTCAGGCTC‑‑GACCTCTCTCCGGTGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGGCAACCTATCCCGGTATc                                                                                                               <  6:223496/146‑1 (MQ=255)                                     cAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCg                                                                                                        >  2:307406/1‑149 (MQ=255)                                      aaGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCgg                                                                                                       <  4:262476/149‑1 (MQ=255)                                              ctgTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTa                                                                                                 <  8:150151/147‑1 (MQ=255)                                                 tCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCa                                                                                             <  1:307406/148‑1 (MQ=255)                                                   gACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTg                                                                                          <  5:198878/149‑1 (MQ=255)                                                      cTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTgcgc                                                                                       >  2:270390/1‑149 (MQ=255)                                                      cTTAAACCACTGGTATTGCTCGCGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGt                                                                                                      >  2:27903/1‑134 (MQ=255)                                                       ttAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGcgcg                                                                                      >  6:74936/1‑149 (MQ=255)                                                             cACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAgg                                                                                <  8:190175/149‑1 (MQ=255)                                                              aCTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAggg                                                                               >  2:7719/1‑149 (MQ=255)                                                                       gCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACg                                                                       >  1:24342/1‑148 (MQ=255)                                                                        cTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACg                                                                       >  3:211959/1‑147 (MQ=255)                                                                        cTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTa                                                                     <  1:246953/149‑1 (MQ=255)                                                                         tCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAg                                                                    >  5:4940/1‑149 (MQ=255)                                                                             gatagACAGACTC‑‑GACCTCACTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCAGAACGCATGGGCTGCCAACCTATCCCGGTATCGCATGCTGTAACTATTCAGTGCGCGGCAACGGCTCAg                                                                          >  1:45750/5‑139 (MQ=255)                                                                             gCTTGACAGTCTC‑‑GACCTCTCTGCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGGCATCCTATCGGGGTATCGGTTGCGGTAACTATGCAGTGCGCTGCAAGGGCTCAg                                                                          <  2:45750/139‑1 (MQ=255)                                                                               ttGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCt                                                              >  8:169863/1‑149 (MQ=255)                                                                               ttGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGAATGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCt                                                              >  4:159553/1‑149 (MQ=255)                                                                                            gACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCttt                                                     <  7:273085/147‑1 (MQ=255)                                                                                            gACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCttt                                                     <  2:274527/147‑1 (MQ=255)                                                                                             aCCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGcc                                                  <  5:35119/149‑1 (MQ=255)                                                                                              ccTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCt                                                 <  2:24342/149‑1 (MQ=255)                                                                                               ctctctCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTg                                                >  4:220284/1‑149 (MQ=255)                                                                                                 ctctcCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCt                                                 >  6:155910/1‑146 (MQ=255)                                                                                                 ctctcCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCt                                                 <  5:155910/146‑1 (MQ=255)                                                                                                      cGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCa                                                                           >  8:271045/1‑115 (MQ=255)                                                                                                         ggTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCca                                      <  1:270390/149‑1 (MQ=255)                                                                                                           ttCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCcaca                                    <  4:211959/149‑1 (MQ=255)                                                                                                              gCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGAt                                 <  6:4940/149‑1 (MQ=255)                                                                                                                   gCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGga                                              <  8:151603/131‑1 (MQ=255)                                                                                                                   gCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGga                                              >  7:151603/1‑131 (MQ=255)                                                                                                                   gCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTcac                            >  7:197509/1‑149 (MQ=255)                                                                                                                   gCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTcac                            >  4:257268/1‑149 (MQ=255)                                                                                                                      gAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTcaccg                          >  4:23343/1‑148 (MQ=255)                                                                                                                      gAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTcaccg                          >  1:149845/1‑148 (MQ=255)                                                                                                                           aaCAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCgg                    >  6:75495/1‑149 (MQ=255)                                                                                                                            aCAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGc                   >  5:288182/1‑149 (MQ=255)                                                                                                                                  aTTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGa             <  1:7719/149‑1 (MQ=255)                                                                                                                                   ttGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGaa            <  5:74936/149‑1 (MQ=255)                                                                                                                                   ttGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGaa            <  3:257268/149‑1 (MQ=255)                                                                                                                                     gCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAAc           <  5:75495/148‑1 (MQ=255)                                                                                                                                       gcGTTCGGATAACCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGAt        >  4:272133/2‑149 (MQ=255)                                                                                                                                         gTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGAtt       >  3:23818/1‑148 (MQ=255)                                                                                                                                         gTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTc      <  3:78149/149‑1 (MQ=255)                                                                                                                                            cGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGc                                                          <  5:197403/94‑1 (MQ=255)                                                                                                                                            cGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGc                                                          >  6:197403/1‑94 (MQ=255)                                                                                                                                             ggTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTCCATc  <  3:159553/149‑1 (MQ=255)                                                                                                                                              gTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTCCATc  >  3:320676/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                   TGCATTCGGCTACCAAGTGTAGAGCATGACTCGTGCCAAGCCCTGCTGTCGACCTTAAACCACTGGTATTGCTCGTGCTTGACAGTCTC‑‑GACCTCTCTCCGGGGTTCGCCAGGCCGAATAAACAGCGATTGCTCGTTCGGTTATCCGAACGCATGGGCTGCCAACCTATCCCGGTATCGCTTGCGGTAACTATGCAGTGCGCGGCAAGGGCTCAGACGTAGTGAGCTTTGCTCTTTGCCTGGAGAGTTCCACAGATGTCACCGCACCGGCTCCCGAACGATTCCATC  >  NC_002947/1226567‑1226853

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