Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,839,227	A→C	5.1%	F99C (TTC→TGC)	ecd ←	flavoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,839,227	0	A	C	5.1%	253.4 / 2.9	79	F99C (TTC→TGC)	ecd	flavoprotein
Reads supporting (aligned to +/- strand):  ref base A (44/31);  new base C (2/2);  total (46/33)
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.36e-01

CGAAAGCTAGTGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTGGGCCGCCGCGATCGCCCAGCTC  >  NC_002947/2839081‑2839373                                                                                                                                                   |                                                                                                                                                   cGAAAGCTAGTGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGCAc                                                                                                                                                  >  4:18924/1‑149 (MQ=255) cGAAAGCTAGTGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGCAc                                                                                                                                                  >  4:122286/1‑149 (MQ=255)  gAAAGCTAGTGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGaa                                                                                                                                                   <  6:311116/147‑1 (MQ=255)  gAAAGCTAGTGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACa                                                                                                                                                 >  4:8299/1‑149 (MQ=255)     aGCTAGTGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCg                                                                                                                                              <  5:117921/149‑1 (MQ=255)          gTGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGc                                                                                                                                          <  7:117087/148‑1 (MQ=255)           tGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCt                                                                                                                                        <  5:291332/149‑1 (MQ=255)            gcggcgTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTg                                                                                                                                       >  3:75258/1‑149 (MQ=255)             cggcgTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGc                                                                                                                                      <  4:75258/149‑1 (MQ=255)               gcgTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTc                                                                                                                                    >  6:4843/1‑149 (MQ=255)               gcgTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTc                                                                                                                                    >  4:49903/1‑149 (MQ=255)                    cTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGCACCGCTGTTCCTGGTctggc                                                                                                                                <  3:18924/148‑5 (MQ=255)                    cTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGCACCGCTGTTCCTGGTctggc                                                                                                                                <  3:122286/148‑5 (MQ=255)                    cTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACa                                                                                                                                <  3:116256/148‑1 (MQ=255)                      gAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATgg                                                                                                                             >  4:81732/1‑149 (MQ=255)                       aGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGgc                                                                                                                            >  3:84985/1‑149 (MQ=255)                                agCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTAc                                                                                                                   >  2:352494/1‑149 (MQ=255)                                  cTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCg                                                                                                                 >  5:20670/1‑149 (MQ=255)                                   tGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCg                                                                                                                 >  5:256586/1‑148 (MQ=255)                                         tACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTc                                                                                                          >  5:79232/1‑149 (MQ=255)                                             acacCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGc                                                                                                      <  3:49903/149‑1 (MQ=255)                                               acCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGTGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCaa                                                                                                    <  1:73962/149‑1 (MQ=255)                                                 cGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGc                                                                                                  >  7:138456/1‑149 (MQ=255)                                                    cGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTc                                                                                               >  6:314554/1‑149 (MQ=255)                                                      gtgtAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTcgc                                                                                             <  3:81732/149‑1 (MQ=255)                                                       tgtAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTcgcg                                                                                            >  6:319990/1‑149 (MQ=255)                                                        gtAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCga                                                                                           <  5:272011/149‑1 (MQ=255)                                                         tAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCgag                                                                                          <  6:206319/149‑1 (MQ=255)                                                             tttCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGc                                                                                      >  6:40563/1‑149 (MQ=255)                                                             tttCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGc                                                                                      >  6:378182/1‑149 (MQ=255)                                                              ttCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACGGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCt                                                                                     >  2:50348/1‑149 (MQ=255)                                                                cTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTc                                                                                    <  3:8299/148‑1 (MQ=255)                                                                  aaaTAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTcgc                                                                                  >  4:52272/1‑148 (MQ=255)                                                                     tAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATc                                                                              >  5:139986/1‑149 (MQ=255)                                                                     tAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATc                                                                              >  3:38810/1‑149 (MQ=255)                                                                      aGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCGCGCGATcc                                                                             >  6:211857/1‑149 (MQ=255)                                                                          gATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCa                                                                         >  7:275200/1‑149 (MQ=255)                                                                          gATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTAGGCGTAGGCGTGACCAGTCCACGCGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCa                                                                         >  1:265387/1‑149 (MQ=255)                                                                               cTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGc                                                                     >  7:121562/1‑148 (MQ=255)                                                                                tGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGc                                                                     >  3:362355/1‑147 (MQ=255)                                                                                tGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGc                                                                     >  6:50714/1‑147 (MQ=255)                                                                                tGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGc                                                                     >  5:245708/1‑147 (MQ=255)                                                                                  cAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGa                                                                   <  4:38810/147‑1 (MQ=255)                                                                                   aGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCa                                                                <  5:40563/149‑1 (MQ=255)                                                                                     tCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCa                                                                <  8:138456/147‑1 (MQ=255)                                                                                       ggTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCt                                                            >  7:121783/1‑149 (MQ=255)                                                                                       ggTGCCTGAAGGCTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCt                                                            >  8:215136/1‑149 (MQ=255)                                                                                           ccTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGACGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTcc                                                          >  4:99399/1‑147 (MQ=255)                                                                                             tGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGc                                                      >  3:207198/1‑149 (MQ=255)                                                                                                   gTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAg                                                <  6:256586/149‑1 (MQ=255)                                                                                                    ttCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAgg                                               >  6:102388/1‑149 (MQ=255)                                                                                                       tGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTAGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTc                                            >  5:2198/1‑149 (MQ=255)                                                                                                         gCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGTAAGAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACTGCGAGCACGCTCCAAGCTTGCAGGTTCGa                                          <  7:258337/149‑1 (MQ=255)                                                                                                         gCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGc                                                  >  1:124246/1‑141 (MQ=255)                                                                                                         gCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGc                                                  <  2:124246/141‑1 (MQ=255)                                                                                                         gCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTc                                            <  3:99399/147‑1 (MQ=255)                                                                                                         gCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGa                                          <  5:333438/149‑1 (MQ=255)                                                                                                          cTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAg                                         <  8:121783/149‑1 (MQ=255)                                                                                                          cTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAg                                         <  4:207198/149‑1 (MQ=255)                                                                                                          cTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAg                                         >  2:257636/1‑149 (MQ=255)                                                                                                           tAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAg                                         >  4:80418/1‑148 (MQ=255)                                                                                                             gTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGaa                                      <  5:4843/149‑1 (MQ=255)                                                                                                              tCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTAGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGaa                                      <  6:2198/148‑1 (MQ=255)                                                                                                                ggCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACGGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTc                                   <  1:50348/149‑1 (MQ=255)                                                                                                                ggCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGt                                    <  3:52272/148‑1 (MQ=255)                                                                                                                 gCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTc                                   >  8:332067/1‑148 (MQ=255)                                                                                                                 gCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCg                                  >  3:328452/1‑149 (MQ=255)                                                                                                                  cGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGa                                 <  5:314554/149‑1 (MQ=255)                                                                                                                   gTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAg                                <  6:83358/149‑1 (MQ=255)                                                                                                                        cGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGa                           >  4:294737/1‑149 (MQ=255)                                                                                                                         gTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGAc                          >  8:137809/1‑149 (MQ=255)                                                                                                                           gACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTg                        >  6:281756/1‑149 (MQ=255)                                                                                                                            aCCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTgg                       <  8:121562/149‑1 (MQ=255)                                                                                                                              cAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTGGgc                     >  3:136479/1‑149 (MQ=255)                                                                                                                                   cACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTTGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTGGGccgccg                >  1:103982/1‑149 (MQ=255)                                                                                                                                   cACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTgg                       <  2:187325/142‑1 (MQ=255)                                                                                                                                   cACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTgg                       >  1:187325/1‑142 (MQ=255)                                                                                                                                         ccAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTGGGCCGCCGCGATCg          >  8:266552/1‑149 (MQ=255)                                                                                                                                                aGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTGGGCCGCCGCGATCGCCCAGCt   <  4:136479/149‑1 (MQ=255)                                                                                                                                                 ggAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGTCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTGGGCCGCCGCGATCGCCCAGCTc  <  7:70765/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CGAAAGCTAGTGCGGCGTTCTGAGCGGCCAGAGCTGCATCTACAACACCGACGGTGTAGCTTTCTAAATAGTCGATACCTGCAGTCGGTGCCTGAAGGGTTCTGGCTAGTCGGCGTAGGCGTGACCAGTCCACGAGCCAGACCAGGAACAGCGGTGCCTGCTCGACATGGCGCTGGTTACCGGACAGTCGGGCAAGCCTCGCGAGACGCTCGCGATCCCGCACGGCGAGCACGCTCCAAGCTTGCAGGTTCGAGGAAGTCGAGGCTGACTGGGCCGCCGCGATCGCCCAGCTC  >  NC_002947/2839081‑2839373

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