Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	345,707	C→A	5.8%	L713L (CTG→CTT)	PP_0287 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	345,707	0	C	A	5.8%	279.9 / 5.9	86	L713L (CTG→CTT)	PP_0287	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (42/39);  new base A (3/2);  total (45/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 = 2.99e-01

TCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGG  >  NC_002947/345559‑345854                                                                                                                                                      |                                                                                                                                                    tCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGcc                                                                                                                                                     >  4:372603/1‑149 (MQ=255)     aCCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCt                                                                                                                                                 >  5:222107/1‑149 (MQ=255)     aCCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCt                                                                                                                                                 >  2:244713/1‑149 (MQ=255)       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                                                                                                          <  5:323071/149‑1 (MQ=255)                                            tCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAggcggc                                                                                                          <  7:165342/149‑1 (MQ=255)                                            tCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAggcgg                                                                                                           >  2:11187/1‑148 (MQ=255)                                              ggTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAg                                                                                                        >  1:182939/1‑149 (MQ=255)                                                ttGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGc                                                                                                      >  5:154606/1‑149 (MQ=255)                                                ttGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGc                                                                                                      <  2:393764/149‑1 (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                                                       >  7:91145/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                                                     >  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)                                                                                                                                                      |                                                                                                                                                    TCATACCAGTCCAGCACTGCGCTGGAGAACTGCTCGGGGCTCATCGGTTGAACAGCCCCTTGAGCGCGTCTTTCAGTTCCGGGCTCACTTTGTCTCCAAGTTTTTCG‑TCAAGTTTTTCATCGATCTTGTCTTTCAGGCGGTTGCCGGCCAGCTTGGCGGCAACCTTGCCCAGGCCGTTCTGGTCGAGGCGGCAGGCCTTTGCGCCGAGTTCCAGCGGACCACGGCAGCGCAGGGGGACTTCAACACCCACATAACGTTCGTTGACCTGGCAGGCCGGGTCCGGCATGGCACGCTGG  >  NC_002947/345559‑345854

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