Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,189,751	A→C	6.4%	L206L (CTA→CTC)	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,189,751	0	A	C	6.4%	230.9 / 6.6	78	L206L (CTA→CTC)	PP_1937	helicase
Reads supporting (aligned to +/- strand):  ref base A (36/36);  new base C (3/2);  total (39/39)
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.24e-01

GACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGCAACGCCAAGGCA  >  NC_002947/2189612‑2189896                                                                                                                                             |                                                                                                                                                  gACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCtt                                                                                                                                          <  5:52758/149‑1 (MQ=255)  aCCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTCTAGCAGCtca                                                                                                                                         <  1:48549/149‑3 (MQ=255)    cAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTCTAGCAGCtcaat                                                                                                                                       <  5:349361/149‑5 (MQ=255)           aCATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTc                                                                                                                                 <  5:101656/148‑1 (MQ=255)           aCATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTc                                                                                                                                 <  6:215067/148‑1 (MQ=255)           aCATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCg                                                                                                                                >  2:330654/1‑149 (MQ=255)            cATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCTCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTAGAAAGTTCg                                                                                                                                >  7:127045/1‑148 (MQ=255)             aTACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGt                                                                                                                               >  3:344295/1‑148 (MQ=255)             aTACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTc                                                                                                                              <  2:74235/149‑1 (MQ=255)                cTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCAc                                                                                                                            <  5:132628/148‑1 (MQ=255)                 tCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCAcca                                                                                                                          >  3:164471/1‑148 (MQ=255)                 tCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCt                                                                                                                          >  2:304329/1‑149 (MQ=255)                 tCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCt                                                                                                                          >  7:99786/1‑149 (MQ=255)                   aGCTTGATGCGAACGTGCGCTTCACGTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCGTTTGGAGCTTGGGTTGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTg                                                                                                                         <  4:37071/148‑1 (MQ=255)                   aGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGCACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCAcc                                                                                                                           <  1:78601/147‑1 (MQ=255)                   aGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCAc                                                                                                                            <  4:164471/145‑1 (MQ=255)                        gAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCTAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCaa                                                                                                                     >  4:340214/1‑147 (MQ=255)                         aaGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAagag                                                                                                                  <  6:143819/149‑1 (MQ=255)                           gCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGc                                                                                                                       <  2:325471/142‑1 (MQ=255)                           gCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAagagca                                                                                                                >  1:325471/1‑147 (MQ=255)                           gCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGc                                                                                                                >  6:150581/1‑149 (MQ=255)                               aCGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCAcc                                                                                                                           <  4:291507/134‑1 (MQ=255)                               aCGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCAcc                                                                                                                           >  3:291507/1‑134 (MQ=255)                                 gTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTg                                                                                                            <  7:306509/147‑1 (MQ=255)                                 gTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTg                                                                                                            >  8:306509/1‑147 (MQ=255)                                         cACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCttt                                                                                                             >  7:309425/1‑138 (MQ=255)                                         cACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCttt                                                                                                             <  8:309425/138‑1 (MQ=255)                                            ctctGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCg                                                                                               >  3:226831/1‑149 (MQ=255)                                             tctGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCg                                                                                               >  1:115786/1‑148 (MQ=255)                                                 gggAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGgaga                                                                                           <  6:249103/148‑1 (MQ=255)                                                 gggAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGagag                                                                                          <  1:330654/149‑1 (MQ=255)                                                    aaaGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGTGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGATGCTTTGTCTAGGGCTTGCTGAGAgatg                                                                                       >  2:119531/1‑149 (MQ=255)                                                    aaaGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAgctg                                                                                       >  8:305472/1‑149 (MQ=255)                                                        agCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTc                                                                                    >  3:303618/1‑148 (MQ=255)                                                           gcgcTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGaaa                                                                                <  1:304329/149‑1 (MQ=255)                                                               tttGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTCGCGGAGAGCTGCTCGAAAGCGa                                                                            >  3:231630/1‑149 (MQ=255)                                                               tttGTTGAGACCAATCTCCTAGATGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGATGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCTAAAGCGa                                                                            >  7:349137/1‑149 (MQ=255)                                                                  gttgAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGAtt                                                                          <  8:99786/148‑1 (MQ=255)                                                                  gttgAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCAGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGAATg                                                                         >  5:27552/1‑149 (MQ=255)                                                                  gttgAGACCAATCTCCTAGCTGAGTGCTATTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCTTCACCTGGCAAGAGGCTTTGCCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTg                                                                         >  4:105268/1‑149 (MQ=255)                                                                   ttgAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTgg                                                                        <  1:315232/149‑1 (MQ=255)                                                                   ttgAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTg                                                                         <  4:226831/148‑1 (MQ=255)                                                                   ttgAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGCGCTTGCGGAGAGCTGCTCGAAAGCGATTgg                                                                        >  8:122215/1‑149 (MQ=255)                                                                          cAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGt                                                                 <  5:131438/149‑1 (MQ=255)                                                                           aaTCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTc                                                                >  6:303416/1‑149 (MQ=255)                                                                            aTCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATAGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCt                                                               >  6:35317/1‑149 (MQ=255)                                                                                ccTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTAtt                                                           >  3:79258/1‑149 (MQ=255)                                                                                     cTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGAt                                                      <  7:122215/149‑1 (MQ=255)                                                                                        aGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTc                                                    <  5:73798/148‑1 (MQ=255)                                                                                        aGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTc                                                    <  6:159695/148‑1 (MQ=255)                                                                                            cTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTa                                                <  2:115786/148‑1 (MQ=255)                                                                                             ttttGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTAtt                                              >  3:181940/1‑149 (MQ=255)                                                                                             ttttGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTAtt                                              <  7:266168/149‑1 (MQ=255)                                                                                             ttttGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTAtt                                              >  5:102884/1‑149 (MQ=255)                                                                                               ttGGAGCTTGTGTAGAAACGCCGAAT‑ACGCTTTCATTGAGATTCTAGAGCAGCTTCTCAAGTTAGTCACCTGGCAAGTGGCTTTTTCTAGGGCTTGAGGAGAGGTGCTCGAAAGCGATTGTGCCAAGTCTTATTTGGATTCGTTATTc                                             <  5:309060/148‑1 (MQ=255)                                                                                               ttGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCg                                            >  7:115446/1‑149 (MQ=255)                                                                                                 ggAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGa                                           <  4:303618/148‑1 (MQ=255)                                                                                                  gAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGATAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAgg                                         <  3:340214/149‑1 (MQ=255)                                                                                                   aGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGc                                        <  7:305472/149‑1 (MQ=255)                                                                                                     cTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCg                                      <  5:303416/149‑1 (MQ=255)                                                                                                        ggggAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGc                                                                                                          <  8:354471/78‑1 (MQ=255)                                                                                                        ggggAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGc                                                                                                          >  7:354471/1‑78 (MQ=255)                                                                                                        ggggAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGAt                                    <  5:35317/148‑1 (MQ=255)                                                                                                           gagaAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGc                                                                                                          <  4:232485/75‑1 (MQ=255)                                                                                                           gagaAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGc                                                                                                          >  3:232485/1‑75 (MQ=255)                                                                                                           gagaAACGCCCAAG‑ACGCTTTCATTCAGCTGCTAGAGCAGCTACTCAAGTTCGTCACCAGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTCCTCGAAAGCGATTGGGCCAAGTCTTATTTGCATTCCTTATTCGAGGCCGATCCGt                                >  2:147017/1‑149 (MQ=255)                                                                                                                   ccGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCtg                        <  4:79258/149‑1 (MQ=255)                                                                                                                             cgttCATTGCGCTGCTAGAGCAGCGTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCg               <  4:344295/147‑1 (MQ=255)                                                                                                                              tttCATTGATCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCtgt                       <  5:112437/140‑1 (MQ=255)                                                                                                                              tttCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGAAAGAGGCTTTGGAGAGGGGTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTAGTATTTGGATTCGTTATTCGATGCCGATACGTTAAATCtgt                       >  6:112437/1‑140 (MQ=255)                                                                                                                               cttgagaaGCTGCTCTAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGc              >  2:320503/8‑148 (MQ=255)                                                                                                                                tCATTGAGCTGCTAGAGCAGCTTCTCAAGTTAGTCACATGGCAAGAGGCTTTGGCGAGGGATTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTAGATGCCGATCCGTTAAATCTCTGGCGATCGCaa            >  6:356918/1‑149 (MQ=255)                                                                                                                                  gagaaGCTGCTCTAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGCAACg          >  8:108142/5‑149 (MQ=255)                                                                                                                                  gagaaGCTGCTCTAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGCAACg          >  8:313639/5‑149 (MQ=255)                                                                                                                                  aTTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGCAACg          >  3:272474/1‑149 (MQ=255)                                                                                                                                      aGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGCAACGcc        >  2:180989/1‑147 (MQ=255)                                                                                                                                       gctgctTGAGCAGCTTCTCAAGTTCGTCACCTTGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGCAACGCCAAg     <  8:349137/149‑1 (MQ=255)                                                                                                                                          gctAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTCGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGCAACGCCAAGGCa  >  6:272946/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                  GACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAG‑ACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAGGCTTTGGCGAGGGCTTGCGGAGAGCTGCTCGAAAGCGATTGGGCCAAGTCTTATTTGGATTCGTTATTCGAGGCCGATCCGTTAAATCTGTGGCCATCGCAACGCCAAGGCA  >  NC_002947/2189612‑2189896

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑