Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,366,286	G→A	8.2%	V90V (GTC→GTT)	ftsH ←	integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,366,286	0	G	A	8.2%	174.8 / 4.3	61	V90V (GTC→GTT)	ftsH	integral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (29/27);  new base A (2/3);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAAGCTCATCGGCCCGCCCTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTCGACCTTGCCATCCTTGACCT  >  NC_002947/5366139‑5366431                                                                                                                                                    |                                                                                                                                                   gAAGCTCATCGGCCCGCCCTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTcgacga                                                                                                                                                   <  7:242970/149‑1 (MQ=255)  aaGCTCATCGGCCCGCCCTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCgacgac                                                                                                                                                  <  2:134255/149‑1 (MQ=255)   aGCTCATCGGCCCGCCCTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGAcc                                                                                                                                                 <  8:259253/149‑1 (MQ=255)    gCTCATCGGCCCGCCCTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCa                                                                                                                                                >  4:155918/1‑149 (MQ=255)                gcccTTGCCGCCCTCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCTCCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACAACCACGTGGTCGTCGa                                                                                                                                    <  4:119789/149‑1 (MQ=255)                  ccTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCAc                                                                                                                                               >  8:177831/1‑136 (MQ=255)                  ccTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCAc                                                                                                                                               <  7:177831/136‑1 (MQ=255)                  ccTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑c                                                                                                                                  >  8:132564/1‑148 (MQ=255)                  ccTTGCCGCCAGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑cc                                                                                                                                 <  5:109297/149‑1 (MQ=255)                        cgcccgcGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCg                                                                                                                           >  8:16841/1‑149 (MQ=255)                         gcccgcGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTTAGGCTGCTTGCCTTCGACAACCACGTGGTCGTCGaa                                                                                                                                   <  1:215139/141‑2 (MQ=255)                         gcccgcGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACAACCACGTGGTCGTCGaa                                                                                                                                   >  2:215139/1‑140 (MQ=255)                                 cACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAgg                                                                                                                  >  3:286723/1‑149 (MQ=255)                                        cATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGt                                                                                                           <  4:136055/149‑1 (MQ=255)                                         aTCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTc                                                                                                          >  2:298800/1‑149 (MQ=255)                                          tCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCAc                                                                                                                                               >  3:270746/1‑112 (MQ=255)                                          tCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCAc                                                                                                                                               <  4:270746/112‑1 (MQ=255)                                              gcgcATGAAGAACATGAACACGGCAATGATCACCAGGATCTGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGt                                                                                                       <  5:296507/147‑1 (MQ=255)                                              gcgcATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGt                                                                                                       <  3:155918/147‑1 (MQ=255)                                                    gaagaaCATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCg                                                                                               >  2:141530/1‑149 (MQ=255)                                                    gaagaaCATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCg                                                                                               >  5:313398/1‑149 (MQ=255)                                                    gaagaaCATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCg                                                                                               >  5:207493/1‑149 (MQ=255)                                                      agaaCATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCgg                                                                                              <  7:16841/148‑1 (MQ=255)                                                           aTGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGcac                                                                                        <  4:286723/149‑1 (MQ=255)                                                               acacGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTc                                                                                    <  1:298800/149‑1 (MQ=255)                                                                   ggCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTg                                                                                 >  3:272616/1‑148 (MQ=255)                                                                   ggCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTg                                                                                 >  5:317489/1‑148 (MQ=255)                                                                    gCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGaa                                                                               <  4:272616/149‑1 (MQ=255)                                                                              ccAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCa                                                                     >  7:90016/1‑149 (MQ=255)                                                                              ccAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCa                                                                     >  8:244055/1‑149 (MQ=255)                                                                                     cGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTgc                                                              >  3:311369/1‑149 (MQ=255)                                                                                       ggAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTg                                                                           >  7:305305/1‑134 (MQ=255)                                                                                       ggAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTg                                                                           <  8:305305/134‑1 (MQ=255)                                                                                       ggAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTgcgc                                                            >  6:213770/1‑149 (MQ=255)                                                                                       ggAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACAACCACGTGGTCGTCGaa                                                                                                                                   >  7:120164/1‑78 (MQ=255)                                                                                       ggAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACAACCACGTGGTCGTCGaa                                                                                                                                   <  8:120164/79‑2 (MQ=255)                                                                                           gCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGc                                                        <  7:132564/149‑1 (MQ=255)                                                                                           gCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGc                                                        <  6:317489/149‑1 (MQ=255)                                                                                           ccgcGCCCCCAGCAACGGCGCCCTGTTCCGCTGTGGCTCAGGCTGCTTGCCTTCCACGCCCACGTGGTTGTCTA‑CCACATCGCCCATCAGCCCATTGTCGGTCATGGCCGGGCGCACAGTCTTCAAGTTGCCGCCATCGGGGCGCTTGc                                                        <  4:311369/146‑1 (MQ=255)                                                                                             tCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGc                                                        >  3:100621/1‑147 (MQ=255)                                                                                               gcaaacaGAAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAg                                                     <  4:218371/143‑1 (MQ=255)                                                                                                 caccaGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGaaa                                                  >  2:109160/1‑149 (MQ=255)                                                                                                 caccaGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTaa                                                  <  7:244055/149‑1 (MQ=255)                                                                                                  accaGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGAAat                                                 <  8:90016/149‑1 (MQ=255)                                                                                                   ccaGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGcaca                                                                                       <  2:260366/110‑1 (MQ=255)                                                                                                   ccaGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGcaca                                                                                       >  1:260366/1‑110 (MQ=255)                                                                                                   ccaGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAatg                                                <  8:276478/149‑1 (MQ=255)                                                                                                   ccaGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAat                                                 >  7:276478/1‑148 (MQ=255)                                                                                                              gTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATc                                     >  7:88103/1‑149 (MQ=255)                                                                                                                   gATGCTCTGTTGGTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGcc                                                       <  2:294634/126‑1 (MQ=255)                                                                                                                   gATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGcc                                                       >  1:294634/1‑126 (MQ=255)                                                                                                                            ttGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGTAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCt                       <  6:313398/149‑1 (MQ=255)                                                                                                                              gCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTc                      <  6:207493/148‑1 (MQ=255)                                                                                                                              gCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTc                      <  1:109160/148‑1 (MQ=255)                                                                                                                               cTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTCGa                    >  4:147648/1‑149 (MQ=255)                                                                                                                                   ggCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTCGACCtt                <  5:213770/149‑1 (MQ=255)                                                                                                                                    gctgctTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCCCACGATCGACCGTg               >  8:214759/1‑149 (MQ=255)                                                                                                                                     ctgctTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTCGACCTTGc              <  8:88103/149‑1 (MQ=255)                                                                                                                                            ccTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGACAGTAATGATGTAGCCATCGACGGTCCCACGCTCGACGTTGCCATCCtt       >  2:74433/1‑149 (MQ=255)                                                                                                                                             cTTCGACGACCACGTGGTTGTCGAGCCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGACATCGACGGTCACACGCTCGACCTTGCCATCCtt       >  7:47185/1‑149 (MQ=255)                                                                                                                                              ttCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTc                                                                         <  8:126104/81‑1 (MQ=255)                                                                                                                                              ttCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTc                                                                         >  7:126104/1‑81 (MQ=255)                                                                                                                                               tCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTCGACCTTGCCATCCTTGAc    >  3:224664/1‑149 (MQ=255)                                                                                                                                                cgacgaCCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTCGACCTTGCCATCCTTGa     <  1:141530/147‑1 (MQ=255)                                                                                                                                                 gacgacCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACTGTCACACGCTCGACCTTGCCATCCTTGACct  <  7:214759/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GAAGCTCATCGGCCCGCCCTTGCCGCCCGCGCCACCTTGCATCTGGCGCATGAAGAACATGAACACGGCAATGATCACCAGGATCGGGAAGCTCGCCACCAGCAACTGCGTCCAGATGCTCTGTTGCTCAGGCTGCTTGCCTTCGACGACCACGTGGTTGTCGA‑CCAGATCGCCGATCAGGCCATTGTCGGTGATGGCCGGGCGCACAGTCTTGAAGTTGTCGCCATCGGTGCGCTTGCCAGTAATGATGTAGCCATCGACGGTCACACGCTCGACCTTGCCATCCTTGACCT  >  NC_002947/5366139‑5366431

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