Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,368,196	A→G	8.1%	intergenic (+25/‑158)	PP_1193 → / → PP_1194	membrane protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,368,196	0	A	G	8.1%	165.9 / 3.5	61	intergenic (+25/‑158)	PP_1193/PP_1194	membrane protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (3/2);  total (32/29)
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.98e-01

AACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGTGCCTC  >  NC_002947/1368061‑1368344                                                                                                                                         |                                                                                                                                                     aaCCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGtttt                                                                                                                                         <  6:276819/149‑1 (MQ=255)       gTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGa                                                                                                                                   <  8:190804/149‑1 (MQ=255)        tGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAg                                                                                                                                  <  8:286814/149‑1 (MQ=255)        tGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAg                                                                                                                                  >  6:147944/1‑149 (MQ=255)          aCGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCa                                                                                                                                <  8:98852/149‑1 (MQ=255)            gATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAgg                                                                                                                              >  7:264973/1‑149 (MQ=255)              tCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTa                                                                                                                                      >  8:259217/1‑139 (MQ=255)              tCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTa                                                                                                                                      <  7:259217/139‑1 (MQ=255)              tCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGtt                                                                                                                            <  2:177298/149‑1 (MQ=255)                     aaCCATGTGGGTGGCGAACCCGGTAAGTGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTATgg                                                                                                                                   >  1:3748/1‑135 (MQ=255)                     aaCCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCGTCCGCTGGGTTTTGTATgg                                                                                                                                   <  2:3748/135‑2 (MQ=255)                         aTGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCt                                                                                                                 <  8:264973/149‑1 (MQ=255)                          tgtgGGTTGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCtt                                                                                                                <  4:58795/149‑1 (MQ=255)                          tgtgGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGCCCCGAAGTGGCCCATACAAACCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTGtt                                                                                                                >  6:113420/1‑149 (MQ=255)                           gtgggtggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGtt                                                                                                                            <  6:160821/136‑1 (MQ=255)                           gtgggtggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGtt                                                                                                                            >  5:160821/1‑136 (MQ=255)                               gtggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCtt                                                                                                           <  2:293423/149‑1 (MQ=255)                               gtggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCtt                                                                                                           >  1:293423/1‑149 (MQ=255)                                tggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCttt                                                                                                          <  1:60775/149‑1 (MQ=255)                                tggCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCtt                                                                                                           >  2:60775/1‑148 (MQ=255)                                          ggTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGCCGGGGTTTGTATGAGCAGGGTCTGGCCTGGCTGGCCTTTATCACGGCag                                                                                                >  2:128766/1‑148 (MQ=255)                                          ggTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCAt                                                                                                >  7:78749/1‑149 (MQ=255)                                            tAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACTTCAttt                                                                                              >  4:246413/1‑149 (MQ=255)                                                    cACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAAt                                                                                         >  1:150158/1‑147 (MQ=255)                                                    cACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAAt                                                                                         >  8:175180/1‑147 (MQ=255)                                                    cACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAAt                                                                                         >  3:133928/1‑147 (MQ=255)                                                      cGGTGTTCATGCCGGCGGCGAGGCGGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACCAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAgc                                                                                     <  1:128766/149‑1 (MQ=255)                                                      cGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTc                                                                                                                  >  2:299025/1‑120 (MQ=255)                                                      cGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTc                                                                                                                  <  1:299025/120‑1 (MQ=255)                                                              aTGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAg                                                                             >  5:313503/1‑149 (MQ=255)                                                                gCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATc                                                                                                       <  8:126173/121‑1 (MQ=255)                                                                           ggCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTc                                                                   <  3:85585/146‑1 (MQ=255)                                                                           ggCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTc                                                                   >  4:85585/1‑146 (MQ=255)                                                                             ccGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAgcg                                                                                    <  5:115375/127‑1 (MQ=255)                                                                             ccGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAgcg                                                                                    >  6:115375/1‑127 (MQ=255)                                                                              cgtggATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGTTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCAt                                                             <  5:249363/146‑1 (MQ=255)                                                                                    tGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTc                                                                      <  6:63639/134‑1 (MQ=255)                                                                                    tGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTc                                                                      >  5:63639/1‑134 (MQ=255)                                                                                     gACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGc                                                                                                                       >  4:33055/1‑84 (MQ=255)                                                                                     gACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGc                                                                                                                       <  3:33055/84‑1 (MQ=255)                                                                                     gACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCggg                                                       <  7:175180/148‑1 (MQ=255)                                                                                              agAATCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGc                                             <  2:78779/149‑1 (MQ=255)                                                                                                  gCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGCCAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGtt                                         >  5:177337/1‑149 (MQ=255)                                                                                                      aTTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCg                                     <  8:78749/149‑1 (MQ=255)                                                                                                            gcTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGAt                               <  6:313503/149‑1 (MQ=255)                                                                                                            gcTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGAt                               <  5:147944/149‑1 (MQ=255)                                                                                                                  ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTc                                                                      <  2:169968/104‑1 (MQ=255)                                                                                                                  ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTc                                                                      >  1:169968/1‑104 (MQ=255)                                                                                                                  ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTgg                                  >  3:146470/1‑140 (MQ=255)                                                                                                                  ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTgg                                  <  4:146470/140‑1 (MQ=255)                                                                                                                  ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTc                         >  1:261553/1‑149 (MQ=255)                                                                                                                  ccGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTc                         >  4:257460/1‑149 (MQ=255)                                                                                                                           cccATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCAGTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATCATTGCGGTTTGCGTGGGATTAATTCCCGCGTCt                 >  5:225049/1‑148 (MQ=255)                                                                                                                            ccATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTta               >  1:329125/1‑149 (MQ=255)                                                                                                                              aTACAAAACCCGTCCGCTGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTtata             >  1:193510/1‑149 (MQ=255)                                                                                                                              aTACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTta               >  3:196735/1‑147 (MQ=255)                                                                                                                                  aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg         <  2:261553/149‑1 (MQ=255)                                                                                                                                  aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg         >  5:133311/1‑149 (MQ=255)                                                                                                                                  aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg         <  4:196735/149‑1 (MQ=255)                                                                                                                                  aaaaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGAGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTgg         >  6:281608/1‑149 (MQ=255)                                                                                                                                    aaCCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGt       >  7:26897/1‑149 (MQ=255)                                                                                                                                         aGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGTGCctc  >  4:167127/1‑149 (MQ=255)                                                                                                                                         |                                                                                                                                                     AACCTGGTGACGATCACGGCAACCATGTGGGTGGCGAACCCGGTAA‑TGACCACGGTGTTCATGCCGGCGGCGAGGCCGGGGATGACAACAGCAGAAGCTGATTGGCGCTGACCCGAAGTGGCCCATACAAAACCCAGCGGACGGGTTTTGTATGAGCAGGTTCTGGCCTGTCTTGCCTTTATCACGTCATTTCCAATGAGCGCTGACAGCTTTTTCCTCGTCCATGCCGGGATATGATTGCGGTTTGCGTGGGATTAATTCCCGCGTCTTATACTGGGTGCCTC  >  NC_002947/1368061‑1368344

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