Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	345,714	G→C	6.1%	A711G (GCC→GGC)	PP_0287 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	345,714	0	G	C	6.1%	267.7 / 6.5	82	A711G (GCC→GGC)	PP_0287	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (38/39);  new base C (3/2);  total (41/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.15e-01

AGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTCGCCTTC  >  NC_002947/345566‑345862                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                <  7:314910/149‑1 (MQ=255)                       cTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGtt                                                                                                                                >  5:261794/1‑149 (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                                                                                                                >  7:185752/1‑122 (MQ=255)                                                                  tCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAg                                                                                                                <  8:185752/122‑1 (MQ=255)                                                                    aGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACg                                                                                   <  1:291790/149‑1 (MQ=255)                                                                     gTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACgg                                                                                  <  8:31561/149‑1 (MQ=255)                                                                      ttCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACgg                                                                                  >  6:311906/1‑148 (MQ=255)                                                                      ttCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACgg                                                                                  >  1:42572/1‑148 (MQ=255)                                                                      ttCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACgg                                                                                  >  6:139705/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                                                               >  4:282705/1‑149 (MQ=255)                                                                                        tcCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCaa                                                               >  2:141601/1‑149 (MQ=255)                                                                                        tcCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCaa                                                               >  7:91145/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                                                             >  6:367132/1‑149 (MQ=255)                                                                                          cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAaca                                                             <  5:311906/149‑1 (MQ=255)                                                                                          cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAac                                                              <  5:281591/148‑1 (MQ=255)                                                                                          cAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAac                                                              <  6:41789/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        <  8:227086/149‑1 (MQ=255)                                                                                                                                                aGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTc        <  6:222107/149‑1 (MQ=255)                                                                                                                                                   ttGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTCGcc     >  5:390066/1‑149 (MQ=255)                                                                                                                                                      gcggcAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTCGCCTTc  >  7:104352/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     AGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGGTCGCCTTC  >  NC_002947/345566‑345862

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