Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	345,712	T→G	5.7%	K712Q (AAG→CAG)	PP_0287 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	345,712	0	T	G	5.7%	278.4 / 5.4	87	K712Q (AAG→CAG)	PP_0287	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (41/41);  new base G (3/2);  total (44/43)
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.57e-01

CAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTCGCC  >  NC_002947/345565‑345859                                                                                                                                                     |                                                                                                                                                    cAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTg                                                                                                                                                    >  1:356409/1‑149 (MQ=255)  aGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTGTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTgg                                                                                                                                                   >  8:328286/1‑149 (MQ=255)    tCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTCGAGCGCGTCTTGCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTggcg                                                                                                                                                 <  6:353628/149‑1 (MQ=255)      cAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTggcgg                                                                                                                                                >  8:165342/1‑148 (MQ=255)       aGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTggcggc                                                                                                                                               >  5:13896/1‑148 (MQ=255)        gCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGCCAAGCTGGCCGGCaa                                                                                                                                             <  6:45555/149‑1 (MQ=255)        gCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGCCAAGCTGGCCGGCaa                                                                                                                                             <  5:88489/149‑1 (MQ=255)         cACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAAc                                                                                                                                            >  2:197012/1‑149 (MQ=255)          aCTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAAcc                                                                                                                                           >  1:393764/1‑149 (MQ=255)               gcTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGcc                                                                                                                                      <  5:15450/149‑1 (MQ=255)                 tGGAGAACTGCTCGGGGATCATCGGTTGAACAGCCCCTTGAGCGCGGCTTTCAGTTCCGGGCTCTCTTTGTCTCCAAGTGTTTCT‑TCAAGTTTTTCATAGAGCATGTGTGTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTAGcaca                                                                                                                                    >  8:121204/1‑149 (MQ=255)                   gagaACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAgg                                                                                                                                  <  6:144077/149‑1 (MQ=255)                    agaACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGc                                                                                                                                 <  3:329484/149‑1 (MQ=255)                       aCTGCTCGGGGCTAATCGGTTTAACAGCCCCTTTAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCGCTCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGcc                                                                                                                                <  1:117463/148‑1 (MQ=255)                        cTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGtt                                                                                                                             >  5:261794/1‑149 (MQ=255)                        cTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGtt                                                                                                                             <  7:314910/149‑1 (MQ=255)                         tGCTCGGGGCTCATCGGTTGAACAGCCCCTTGACCGCGTCTTACACTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGAGCTTGTCTATCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTc                                                                                                                            >  5:41789/1‑149 (MQ=255)                           cTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTg                                                                                                                          >  1:49934/1‑149 (MQ=255)                           cTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCACGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTg                                                                                                                          >  1:52003/1‑149 (MQ=255)                              ggggCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTgg                                                                                                                         <  4:191843/147‑1 (MQ=255)                                     aTCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAggcgg                                                                                                                >  4:270778/1‑149 (MQ=255)                                      tCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAggcggc                                                                                                               <  7:165342/149‑1 (MQ=255)                                      tCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAggcggc                                                                                                               <  5:323071/149‑1 (MQ=255)                                      tCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAggcgg                                                                                                                >  2:11187/1‑148 (MQ=255)                                        ggTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAg                                                                                                             >  1:182939/1‑149 (MQ=255)                                          ttGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGc                                                                                                           <  2:393764/149‑1 (MQ=255)                                          ttGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGc                                                                                                           >  5:154606/1‑149 (MQ=255)                                              aCAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCttt                                                                                                       >  4:260408/1‑149 (MQ=255)                                                   cccTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGcc                                                                                                  <  2:356409/149‑1 (MQ=255)                                                     cTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCTGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGcc                                                                                                  >  7:155283/1‑147 (MQ=255)                                                     cTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGcc                                                                                                  >  3:186563/1‑147 (MQ=255)                                                     cTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGa                                                                                                <  2:49934/149‑1 (MQ=255)                                                     cTTGAGCGCGTCCTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGcc                                                                                                  >  6:281591/1‑147 (MQ=255)                                                        gAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGt                                                                                              <  8:312675/148‑1 (MQ=255)                                                         aGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTGCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTc                                                                                            <  6:314542/149‑1 (MQ=255)                                                          gcgcGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTc                                                                                            <  6:13896/148‑1 (MQ=255)                                                              gTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCg                                                                                       >  2:291790/1‑149 (MQ=255)                                                                   tCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAg                                                                                                             <  8:185752/122‑1 (MQ=255)                                                                   tCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAg                                                                                                             >  7:185752/1‑122 (MQ=255)                                                                     aGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACg                                                                                <  1:291790/149‑1 (MQ=255)                                                                      gTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACgg                                                                               <  8:31561/149‑1 (MQ=255)                                                                       ttCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACgg                                                                               >  6:139705/1‑148 (MQ=255)                                                                       ttCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACgg                                                                               >  6:311906/1‑148 (MQ=255)                                                                       ttCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACgg                                                                               >  1:42572/1‑148 (MQ=255)                                                                       ttCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGc                                                                              >  2:15433/1‑149 (MQ=255)                                                                         ccGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCa                                                                             <  2:182939/148‑1 (MQ=255)                                                                          cGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAgc                                                                           >  7:325446/1‑149 (MQ=255)                                                                          cGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAg                                                                            <  6:317836/148‑1 (MQ=255)                                                                            ggCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAgcgc                                                                         <  3:260408/149‑1 (MQ=255)                                                                                 aCTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAggg                                                                     <  1:11187/148‑1 (MQ=255)                                                                                 aCTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAgg                                                                      <  6:261794/147‑1 (MQ=255)                                                                                   tttGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTc                                                                                                                            >  5:225694/1‑91 (MQ=255)                                                                                   tttGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTc                                                                                                                            <  6:225694/91‑1 (MQ=255)                                                                                   tttGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGa                                                                  >  5:61127/1‑149 (MQ=255)                                                                                     tGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGAc                                                                 <  8:155283/148‑1 (MQ=255)                                                                                     tGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGAc                                                                 <  3:270778/148‑1 (MQ=255)                                                                                         tcCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCaa                                                            >  7:91145/1‑149 (MQ=255)                                                                                         tcCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCaa                                                            >  4:282705/1‑149 (MQ=255)                                                                                         tcCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCaa                                                            >  2:141601/1‑149 (MQ=255)                                                                                         tcCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCa                                                             >  3:185507/1‑148 (MQ=255)                                                                                           cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGGCGGGCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGGAGGGGGACTTCaa                                                            <  7:121204/147‑1 (MQ=255)                                                                                           cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTggc                                                                                                                                                  <  6:179334/61‑1 (MQ=255)                                                                                           cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAaca                                                          <  5:311906/149‑1 (MQ=255)                                                                                           cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAaca                                                          >  6:367132/1‑149 (MQ=255)                                                                                           cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAac                                                           <  6:41789/148‑1 (MQ=255)                                                                                           cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAac                                                           <  5:281591/148‑1 (MQ=255)                                                                                                 tttCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTGGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCcac                                                     <  6:293794/148‑1 (MQ=255)                                                                                                       tCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcg                                               >  6:259545/1‑149 (MQ=255)                                                                                                        cAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgt                                              >  3:333506/1‑149 (MQ=255)                                                                                                              tttCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAAcgttcgtt                                         <  2:52003/148‑1 (MQ=255)                                                                                                                 cATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGAc                                      <  7:328286/148‑1 (MQ=255)                                                                                                                  aTCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGAcc                                     <  4:185507/148‑1 (MQ=255)                                                                                                                  aTCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGAcc                                     <  2:42572/148‑1 (MQ=255)                                                                                                                  aTCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCt                                    >  3:306993/1‑149 (MQ=255)                                                                                                                   tCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTg                                   <  1:197012/149‑1 (MQ=255)                                                                                                                          tGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGcc                            >  7:227086/1‑149 (MQ=255)                                                                                                                                  aGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCg                     <  1:244713/148‑1 (MQ=255)                                                                                                                                    aaggTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCa                  >  4:26932/3‑149 (MQ=255)                                                                                                                                      ggTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATg                <  5:367132/149‑1 (MQ=255)                                                                                                                                      ggTTGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATg                >  2:249184/1‑149 (MQ=255)                                                                                                                                        ttGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATgg               >  2:184761/1‑148 (MQ=255)                                                                                                                                        ttGCCGCCAAGCTGGCCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGCACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGCGTCCGGCATgg               >  8:336271/1‑148 (MQ=255)                                                                                                                                             ggCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCt         >  3:221736/1‑149 (MQ=255)                                                                                                                                             ggCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCt         >  1:180933/1‑149 (MQ=255)                                                                                                                                               ccAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTgg       >  2:1043/1‑149 (MQ=255)                                                                                                                                                 aGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTc     <  6:222107/149‑1 (MQ=255)                                                                                                                                                 aGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTc     <  8:227086/149‑1 (MQ=255)                                                                                                                                                    ttGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTCGcc  >  5:390066/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTCGCC  >  NC_002947/345565‑345859

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