Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,324,257	A→G	25.3%	G258G (GGT→GGC)	yfdC ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,324,257	0	A	G	25.3%	118.3 / 44.1	75	G258G (GGT→GGC)	yfdC	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (29/27);  new base G (10/9);  total (39/36)
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.22e-01

CTCACCGTAGCGGCACGCTGAGTATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGAAATCGCCAAGGGCCATCAGGTAGGTGAT  >  NC_002947/1324110‑1324403                                                                                                                                                    |                                                                                                                                                   cTCACCGTAGCGGCACGCTGAGTATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCAc                                                                                                                                                   >  4:240393/1‑149 (MQ=255)  tCACCGTAGCGGCACGCTGAGTATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCAcc                                                                                                                                                  >  8:179233/1‑149 (MQ=255)        tAGCGGCAGGCTGAGTATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCTTTCAGAAGCCGCGAGGGTGGTTTGACACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCgatgat                                                                                                                                            <  6:272779/149‑1 (MQ=255)          gCGGCACGCTGAGTATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGatgatg                                                                                                                                           <  6:163955/148‑1 (MQ=255)            ggCACGCTGAGTATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTg                                                                                                                                        <  3:240393/149‑1 (MQ=255)             gCACGCTGAGTATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGc                                                                                                                                       >  2:169909/1‑149 (MQ=255)                       tATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGt                                                                                                                             >  6:216049/1‑149 (MQ=255)                       tATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGt                                                                                                                             >  7:249017/1‑149 (MQ=255)                       tATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCgatg                                                                                                                                              <  2:87783/132‑1 (MQ=255)                       tATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCgatg                                                                                                                                              >  1:87783/1‑132 (MQ=255)                       tATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGt                                                                                                                             >  8:113553/1‑149 (MQ=255)                            tttGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTGTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGAGCACTTGTGCATGGCTGATCAAGGAGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCgg                                                                                                                          >  1:33697/1‑147 (MQ=255)                               gtgCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGt                                                                                                                      <  8:324196/148‑1 (MQ=255)                                  cTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCa                                                                                                                   >  4:325310/1‑148 (MQ=255)                                  cTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGCGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCTAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCa                                                                                                                   >  8:307700/1‑148 (MQ=255)                                       gTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGTGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTc                                                                                                              >  8:32800/1‑148 (MQ=255)                                       gTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGTAGTGCAAATTc                                                                                                              >  1:266000/1‑148 (MQ=255)                                        tCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTcca                                                                                                            <  5:216049/149‑1 (MQ=255)                                        tCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTcca                                                                                                            <  7:179233/149‑1 (MQ=255)                                        tCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTcca                                                                                                            <  8:185120/149‑1 (MQ=255)                                          tgctgcGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTccagc                                                                                                          >  2:340244/1‑149 (MQ=255)                                              gcgcgcGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTccagccag                                                                                                       >  2:128601/1‑148 (MQ=255)                                                    ggACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGTGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGc                                                                                                <  7:32800/149‑1 (MQ=255)                                                           ggCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGt                                                                                         <  1:340244/149‑1 (MQ=255)                                                              tGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTtc                                                                                      >  7:315944/1‑149 (MQ=255)                                                                gtattCGTTCAGAAGCCGCGAGTGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGTCTGATCAAGGCGAAGATTAAACTGCCACCGTTGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGCTTtcgt                                                                                    <  2:33697/146‑1 (MQ=255)                                                                  cTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTccagc                                                                                                          >  7:91529/1‑125 (MQ=255)                                                                  cTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTccagc                                                                                                          <  8:91529/125‑1 (MQ=255)                                                                   ttCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGc                                                                                 >  3:18945/1‑149 (MQ=255)                                                                   ttCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGc                                                                                 >  5:174297/1‑149 (MQ=255)                                                                   ttCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTAACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGc                                                                                 >  1:237126/1‑149 (MQ=255)                                                                       ttCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGc                                                                             >  5:7752/1‑149 (MQ=255)                                                                         cAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGTAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGcc                                                                            <  2:266000/148‑1 (MQ=255)                                                                         cAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGcc                                                                            <  6:7752/148‑1 (MQ=255)                                                                          aGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCa                                                                               <  7:204971/144‑1 (MQ=255)                                                                          aGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCa                                                                               >  8:204971/1‑144 (MQ=255)                                                                          aGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCca                                                                          <  8:315944/149‑1 (MQ=255)                                                                            aaGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCAAGCCcaca                                                                        >  8:234136/1‑149 (MQ=255)                                                                               ccGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCacac                                                                       >  7:28128/1‑147 (MQ=255)                                                                                 gcgAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAgg                                                                   <  8:249017/149‑1 (MQ=255)                                                                                  cgAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAgg                                                                   >  4:338649/1‑148 (MQ=255)                                                                                   gAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAgg                                                                   >  6:35461/1‑147 (MQ=255)                                                                                    aGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCacac                                                                       <  2:49011/142‑1 (MQ=255)                                                                                    aGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCacac                                                                       >  1:49011/1‑142 (MQ=255)                                                                                    aGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAg                                                                <  7:113553/149‑1 (MQ=255)                                                                                     gggTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAg                                                                >  6:332558/1‑148 (MQ=255)                                                                                      ggtggtTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGa                                                              <  5:332558/149‑1 (MQ=255)                                                                                        tggtTTGACACTGTCACTGCGCACTTGTGCATCGATCATCAATGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCATATTCCAGACACAAGCGCCTCCACGTTTCGTCGCAAGCCCACACCACGTAGGac                                                             >  6:90737/1‑148 (MQ=255)                                                                                          gtTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGAAGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACt                                                          >  1:85123/1‑149 (MQ=255)                                                                                           tttGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTTCCACCGATGATGTTGCCCGCCAGCGTCGGCATTGCAAATTCCAGCCAGAAGCGGCTc                                                                                              <  5:66156/112‑1 (MQ=255)                                                                                           tttGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACtt                                                         <  8:28128/149‑1 (MQ=255)                                                                                           tttGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCTTCGCCAGTGCAAATTCCAGCCAGAACCGGCTc                                                                                              >  6:66156/1‑112 (MQ=255)                                                                                               ccACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGAc                                                     <  5:35461/149‑1 (MQ=255)                                                                                               ccACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGAc                                                     >  1:71752/1‑149 (MQ=255)                                                                                                cACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGt                                                                  <  6:334512/135‑1 (MQ=255)                                                                                                cACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGt                                                                  >  5:334512/1‑135 (MQ=255)                                                                                                 aCTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAg                                                                    <  3:247832/132‑1 (MQ=255)                                                                                                 aCTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAg                                                                    >  4:247832/1‑132 (MQ=255)                                                                                                   tGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGAc                                                 >  1:192399/1‑149 (MQ=255)                                                                                                     tCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCGCCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGAc                                                 <  4:18945/147‑1 (MQ=255)                                                                                                     tCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGAc                                                 <  3:325310/147‑1 (MQ=255)                                                                                                     tCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGAc                                                 <  1:169909/147‑1 (MQ=255)                                                                                                       aCTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCcca                                             >  1:73148/1‑149 (MQ=255)                                                                                                                ttGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTatat                                    >  1:205945/1‑149 (MQ=255)                                                                                                                     cATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGt                                <  2:210406/148‑1 (MQ=255)                                                                                                                         acTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGc                                                                                 >  3:198659/2‑95 (MQ=255)                                                                                                                         acTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGc                                                                                 <  4:198659/94‑1 (MQ=255)                                                                                                                             aTCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGAAATCGc                        <  2:205945/148‑1 (MQ=255)                                                                                                                              tCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTAACCCCACCACTATATGAGTGAAATCGcc                       >  7:341852/1‑148 (MQ=255)                                                                                                                                  ggCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGAAATCGCCAAggg                  <  3:338649/149‑1 (MQ=255)                                                                                                                                       aGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGTAATCGCCAAGGGCCATc             <  2:73148/149‑1 (MQ=255)                                                                                                                                            caccTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGAAATCGCCAAGGGCCATCaggt         <  1:128601/145‑1 (MQ=255)                                                                                                                                             aaCTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAATCGGCTCCAGCTTGCTTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGAAATCGCCAAGGGCCATCaggtag       <  7:234136/149‑1 (MQ=255)                                                                                                                                             aaCTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGAAATCGCCAAGGGCCATCaggtag       <  2:192399/149‑1 (MQ=255)                                                                                                                                                  ccACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGAAATCGCCAAGGGCCATCAGGTAGGTGat  >  2:1677/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CTCACCGTAGCGGCACGCTGAGTATCATTTGTGCTCATGTCTGCTGCGCGCGGACCGGGGCTGGTCTTCGTTCAGAAGCCGCGAGGGTGGTTTGCCACTGTCACTGCGCACTTGTGCATGGCTGATCAAGGCGAAGATGAAACTGCCACCGATGATGTTGCCCGCCAGCGTCGGCAGTGCAAATTCCAGCCAGAAGCGGCTCCAGGTTTCGTCGCCAGCCCACACCAGGTAGGACACTTCGACTGACCCCACCACTATATGAGTGAAATCGCCAAGGGCCATCAGGTAGGTGAT  >  NC_002947/1324110‑1324403

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