Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,579,595	C→A	10.7%	intergenic (+54/+22)	psd → / ← serB	phosphatidylserine decarboxylase/phosphoserine phosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,579,595	0	C	A	10.7%	182.7 / 7.3	60	intergenic (+54/+22)	psd/serB	phosphatidylserine decarboxylase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base C (26/27);  new base A (3/4);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCA  >  NC_002947/5579447‑5579743                                                                                                                                                     |                                                                                                                                                      tGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGgccgcc                                                                                                                                                       >  2:78578/1‑149 (MQ=255)   ttCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCtt                                                                                                                                                     >  1:68577/1‑149 (MQ=255)       ggCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTAGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCgg                                                                                                                                                 >  3:124031/1‑149 (MQ=255)        gCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGc                                                                                                                                                <  7:303871/149‑1 (MQ=255)         ccAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCTGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTGGGGCCGCAAGGCGGcc                                                                                                                                               >  4:289212/1‑149 (MQ=255)         ccAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCTGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTGGGGCCGCAAGGCGGcc                                                                                                                                               >  4:294022/1‑149 (MQ=255)            gagCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCAAGGCGGcccaa                                                                                                                                            <  8:296245/149‑3 (MQ=255)             agCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCTGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTGGGGCCGCAAGGCGGcccaac                                                                                                                                           <  3:289212/149‑4 (MQ=255)             agCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCTGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTGGGGCCGCAAGGCGGcccaac                                                                                                                                           <  3:294022/149‑4 (MQ=255)                   tGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCAc                                                                                                                                     >  6:221403/1‑149 (MQ=255)                   tGAAAGGGGCTGAGCGCCTGGGGGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCAc                                                                                                                                     <  2:236037/149‑1 (MQ=255)                               agagCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGcccc                                                                                                                                             <  8:256764/129‑1 (MQ=255)                               agagCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGcccc                                                                                                                                             >  7:256764/1‑129 (MQ=255)                                 agCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCggg                                                                                                                       <  5:244380/149‑1 (MQ=255)                                    cTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGt                                                                                                                    >  5:29561/1‑149 (MQ=255)                                      gggTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTcg                                                                                                                  >  3:217416/1‑149 (MQ=255)                                           cAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATc                                                                                                             <  2:68577/149‑1 (MQ=255)                                             ggCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCAc                                                                                                           >  5:151956/1‑149 (MQ=255)                                              gCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGCTACGCCGCGGGCGTCGCGATCACg                                                                                                          >  4:16190/1‑149 (MQ=255)                                                tctGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGa                                                                                                               >  7:65275/1‑142 (MQ=255)                                                tctGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGa                                                                                                               <  8:65275/142‑1 (MQ=255)                                                  tGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGt                                                                                                                    >  5:7200/1‑135 (MQ=255)                                                  tGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGt                                                                                                                    <  6:7200/135‑1 (MQ=255)                                                  tGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑gg                                                                                                     >  4:245892/1‑149 (MQ=255)                                                   gCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGc                                                                                                    >  3:179608/1‑149 (MQ=255)                                                         cGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGa                                                                                                                                          >  1:258381/1‑106 (MQ=255)                                                         cGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGa                                                                                                                                          <  2:258381/106‑1 (MQ=255)                                                          gCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCcagca                                                                                              <  4:99620/148‑1 (MQ=255)                                                               gggCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTAc                                                                                         <  5:221403/148‑1 (MQ=255)                                                               gggCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTAc                                                                                         <  5:181906/148‑1 (MQ=255)                                                                    gcagcTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCAcc                                                                                   <  4:293725/149‑1 (MQ=255)                                                                           gcggcgCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCa                                                                            <  4:217416/149‑1 (MQ=255)                                                                            cggcgCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCAc                                                                                                                                     >  6:87991/1‑92 (MQ=255)                                                                            cggcgCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCAc                                                                                                                                     <  5:87991/92‑1 (MQ=255)                                                                              gcggCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGGCCTGCATATGGAGGGCCTTTTTCGTTGGGGCGCCTTGCGGCCCCGGTTTCACGGCTCAGCCGCGGGCGTCGCGATCACGCATGCCCCAGCAGGTACAGCACCCCGTCCAGc                                                                          <  1:4526/149‑1 (MQ=255)                                                                                          gCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGaaa                                                             >  5:309174/1‑149 (MQ=255)                                                                                            ttGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATc                                                           >  3:78051/1‑149 (MQ=255)                                                                                               aTTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGc                                                         <  6:29561/148‑1 (MQ=255)                                                                                                 tCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCt                                                       <  4:124031/148‑1 (MQ=255)                                                                                                  cGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGc                                                     >  1:319816/1‑149 (MQ=255)                                                                                                   gAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCt                                                    >  4:86644/1‑149 (MQ=255)                                                                                                      ccATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTgg                                                 <  4:78051/149‑1 (MQ=255)                                                                                                          atGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGa                                             >  3:77100/1‑149 (MQ=255)                                                                                                            gAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGCCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACt                                           >  8:127620/1‑149 (MQ=255)                                                                                                             aaaaGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACt                                           <  1:78578/148‑1 (MQ=255)                                                                                                               aaGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTg                                          <  3:265767/147‑1 (MQ=255)                                                                                                               aaGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTg                                          <  3:261229/147‑1 (MQ=255)                                                                                                               aaGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTg                                          >  4:265767/1‑147 (MQ=255)                                                                                                               aaGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTg                                          >  4:261229/1‑147 (MQ=255)                                                                                                                    ccTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCAcc                                   <  4:77100/149‑1 (MQ=255)                                                                                                                    ccTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCAc                                    >  2:67711/1‑148 (MQ=255)                                                                                                                            tGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCtt                           <  6:309174/149‑1 (MQ=255)                                                                                                                              cAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTgg                         <  4:179608/149‑1 (MQ=255)                                                                                                                              cAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTg                          <  3:86644/148‑1 (MQ=255)                                                                                                                                     tttttCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAg                                                                    >  3:203909/1‑99 (MQ=255)                                                                                                                                     tttttCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAg                                                                    <  4:203909/99‑1 (MQ=255)                                                                                                                                         atccggggCCGCAAGGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTg                          <  3:47832/132‑1 (MQ=255)                                                                                                                                         atccggggCCGCAAGGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTg                          >  4:47832/6‑137 (MQ=255)                                                                                                                                          cGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAacgcca             >  7:99942/1‑149 (MQ=255)                                                                                                                                           gTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTgg                         >  7:213295/1‑136 (MQ=255)                                                                                                                                           gTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTgg                         <  8:213295/136‑1 (MQ=255)                                                                                                                                                 ccgccTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTgg                         <  2:193860/130‑1 (MQ=255)                                                                                                                                                 ccgccTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTgg                         >  1:193860/1‑130 (MQ=255)                                                                                                                                                     cTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACcca  <  2:319816/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TGTTCGGGCCAGAGCAGGTGAAATGGGCTGAGAGCCTGGGTGCAGGCTCTGCAGTTCGCATGGGGCAGCAGCTGGCGGCGCCTGCGCAGGCTTGATTCGAGCCATATGAAAAGGCCCTGCATATGCAGGGCCTTTTTCGTTGGGCCGCCTTGCGGCCCCGGATTCACGGATCAGCCGCGGGCGTCGCGATCACGCA‑GGCCCAGCAGGTACAGCACCCCGTCCAGCCCCAAGGTGGAAATCGCCTGCTTGGCCGACTGGCGCACCAGCGGCTTGGCACGAAACGCCACGCCCAACCCA  >  NC_002947/5579447‑5579743

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