Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,860,813	G→C	5.7%	L79L (CTC→CTG)	ligA ←	DNA ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,860,813	0	G	C	5.7%	261.3 / 4.5	81	L79L (CTC→CTG)	ligA	DNA ligase
Reads supporting (aligned to +/- strand):  ref base G (34/42);  new base C (3/2);  total (37/44)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATC  >  NC_002947/4860669‑4860959                                                                                                                                                 |                                                                                                                                                    gCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCa                                                                                                                                                 >  5:13416/1‑149 (MQ=255)  ccGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCAt                                                                                                                                                <  5:69978/149‑1 (MQ=255)   cGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATg                                                                                                                                               <  3:309442/149‑1 (MQ=255)   cGTCGCGGTACAGCAGGCGCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCTCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATg                                                                                                                                               <  8:105647/149‑1 (MQ=255)    gTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATgg                                                                                                                                              <  6:34813/149‑1 (MQ=255)    gTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATgg                                                                                                                                              <  5:163316/149‑1 (MQ=255)       gcgGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCATGCTCAGCCTGGGc                                                                                                                                            <  3:116190/148‑1 (MQ=255)              agcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTg                                                                                                                                    <  3:347710/149‑1 (MQ=255)              agcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTg                                                                                                                                    >  1:212790/1‑149 (MQ=255)              agcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTg                                                                                                                                    <  2:217765/149‑1 (MQ=255)              agcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTg                                                                                                                                    <  4:245060/149‑1 (MQ=255)              agcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTg                                                                                                                                    <  3:147202/149‑1 (MQ=255)               gcagGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTgg                                                                                                                                   >  6:96506/1‑149 (MQ=255)                   gCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTa                                                                                                                               >  6:26458/1‑149 (MQ=255)                     tCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTAcc                                                                                                                             >  3:260925/1‑149 (MQ=255)                       aCGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTCCCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTg                                                                                                                           >  2:108047/1‑149 (MQ=255)                       aCGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTg                                                                                                                           >  8:57841/1‑149 (MQ=255)                        cGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTgg                                                                                                                          <  6:26211/149‑1 (MQ=255)                          gccaggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCt                                                                                                                        >  3:76292/1‑149 (MQ=255)                              ggccaTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAc                                                                                                                    <  5:299611/149‑1 (MQ=255)                                 caTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAcgcc                                                                                                                 >  6:57015/1‑149 (MQ=255)                                  aTCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAAcgccg                                                                                                                <  6:121744/149‑1 (MQ=255)                                        cTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAgg                                                                                                          <  3:103869/149‑1 (MQ=255)                                                  cAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccacc                                                                                                <  5:26458/149‑1 (MQ=255)                                                       gTAGTCGCCCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTccacccaccc                                                                                           >  3:232807/1‑147 (MQ=255)                                                       gTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCacac                                                                                           >  5:216469/1‑149 (MQ=255)                                                        tAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACCCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACg                                                                                          >  7:208138/1‑149 (MQ=255)                                                          gTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGc                                                                                         <  4:20080/148‑1 (MQ=255)                                                             gACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgg                                                                                      <  6:13416/148‑1 (MQ=255)                                                              aCCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggt                                                                                    <  3:194050/149‑1 (MQ=255)                                                                 gcgcCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtggg                                                                                 <  2:70590/149‑1 (MQ=255)                                                                  cgcCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGc                                                                                <  2:212790/149‑1 (MQ=255)                                                                  cgcCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGc                                                                                <  4:260925/149‑1 (MQ=255)                                                                    cccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGa                                                                              <  6:221974/149‑1 (MQ=255)                                                                     ccGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAg                                                                             >  2:313806/1‑149 (MQ=255)                                                                        gCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑g                                                                         >  3:151328/1‑149 (MQ=255)                                                                         cTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑gg                                                                        >  5:122101/1‑149 (MQ=255)                                                                         cTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑gg                                                                        <  5:57015/149‑1 (MQ=255)                                                                                aGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGAGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCGGGGTGGGCGCGTC‑GGGGGTTAc                                                                 >  6:194813/1‑149 (MQ=255)                                                                                          ccACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCg                                                         <  3:58170/147‑1 (MQ=255)                                                                                           cACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGgt                                                      <  4:151328/149‑1 (MQ=255)                                                                                           cACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGgt                                                      <  1:18300/149‑1 (MQ=255)                                                                                             cgcgACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTgg                                                                                                                                   >  5:383880/1‑71 (MQ=255)                                                                                             cgcgACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTgg                                                                                                                                   <  6:383880/71‑1 (MQ=255)                                                                                               cgACGGCCGAACTCACGGAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCGGATCATGGGCAGTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACGAGGTGCGGGTGTTc                                                  <  4:87725/149‑1 (MQ=255)                                                                                                 aCGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACCCGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGGTCAg                                                <  4:232807/149‑1 (MQ=255)                                                                                                 aCGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAg                                                <  6:216469/149‑1 (MQ=255)                                                                                                  cGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGc                                               <  5:266745/149‑1 (MQ=255)                                                                                                       gAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTTCCAGGTGCGGGTGTTCAGCTTCCa                                          <  7:57841/149‑1 (MQ=255)                                                                                                           tCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAgcgc                                      >  8:51903/1‑149 (MQ=255)                                                                                                             aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCTAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCt                                    <  8:14767/149‑1 (MQ=255)                                                                                                             aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtg                                                                                   <  2:250615/103‑1 (MQ=255)                                                                                                             aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTAcc                                                                >  1:26338/1‑121 (MQ=255)                                                                                                             aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTAcc                                                                <  2:26338/121‑1 (MQ=255)                                                                                                             aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCt                                    >  4:152891/1‑149 (MQ=255)                                                                                                             aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCt                                    >  3:281640/1‑149 (MQ=255)                                                                                                             aCGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGCCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCtgggtg                                                                                   >  1:250615/1‑103 (MQ=255)                                                                                                              cGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCtt                                   <  6:122101/149‑1 (MQ=255)                                                                                                                cAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGa                                 <  8:208138/149‑1 (MQ=255)                                                                                                                    tCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTc                             >  1:147296/1‑149 (MQ=255)                                                                                                                      ggCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGt                           >  8:296649/1‑149 (MQ=255)                                                                                                                       gCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGgtgt                                                    <  6:22563/123‑1 (MQ=255)                                                                                                                       gCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGgtgt                                                    >  5:22563/1‑123 (MQ=255)                                                                                                                         ttcttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGGTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTg                         >  2:224987/1‑148 (MQ=255)                                                                                                                           cttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGAGGGTGTTCAGCTTCCAGCGACTTGAGCTCGTTGaa                       >  8:371731/1‑148 (MQ=255)                                                                                                                            ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                       >  2:372345/1‑147 (MQ=255)                                                                                                                            ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGaa                       >  7:302424/1‑147 (MQ=255)                                                                                                                            ttcGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAAc                      >  3:277142/1‑148 (MQ=255)                                                                                                                              cGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAAc                      <  4:277142/146‑1 (MQ=255)                                                                                                                                cacgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCg                 >  3:95985/8‑149 (MQ=255)                                                                                                                                  cgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGa                                 <  4:117715/126‑1 (MQ=255)                                                                                                                                  cgaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGa                                 >  3:117715/6‑131 (MQ=255)                                                                                                                                   gaagtGCCCATGCTCAGCCTGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAAGAGGCGGTc              >  4:106314/5‑149 (MQ=255)                                                                                                                                    cggtGCCCAGGCTGTGCATGGGGACTTCGTGGGGTACCTGGCTGTACGCCGCCAAGGCAGGTCCACCCACAAGCTGGGTGGGCGAGTCGGGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTc              <  1:108047/146‑1 (MQ=255)                                                                                                                                     gTTGCTCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCg             >  6:366221/1‑148 (MQ=255)                                                                                                                                     gTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCg             >  1:220387/1‑148 (MQ=255)                                                                                                                                             ggCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                          <  1:27179/127‑1 (MQ=255)                                                                                                                                             ggCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGtt                          >  2:27179/1‑127 (MQ=255)                                                                                                                                             ggCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCa    <  6:356114/149‑1 (MQ=255)                                                                                                                                               cTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATc  <  2:220387/149‑1 (MQ=255)                                                                                                                                                tGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCAt   >  2:205481/1‑147 (MQ=255)                                                                                                                                                 |                                                                                                                                                    GCCGTCGCGGTACAGCAGGCTCACGGCCAGGCCATCGAGCTTGGGTTCGCAGCTGTAGTCGACCGCGCCCGGCTGATCGAGGCCTTCTACCACGCGACGGCCGAACTCACGCAGGTCGGCTTCTTCGAAGGCGTTGCCCAGGCTGAGCATGGGCACTTCGTGGCGTACCTGGCTGAACGCCGCCAAGGCAGCTCCACCCACACGCTGGGTGGGCGAGTC‑GGGGGTTACCAGGTGCGGGTGTTCAGCTTCCAGCGCCTTGAGCTCGTTGAACAGGCGGTCGTATTCGGCATC  >  NC_002947/4860669‑4860959

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