Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,158,984	T→A	10.1%	intergenic (+321/+547)	PP_5739 → / ← PP_5402	hypothetical protein/NAD‑dependent acetyl‑CoA synthetase deacetylase subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,158,984	0	T	A	10.1%	161.4 / 6.3	59	intergenic (+321/+547)	PP_5739/PP_5402	hypothetical protein/NAD‑dependent acetyl‑CoA synthetase deacetylase subunit
Reads supporting (aligned to +/- strand):  ref base T (27/26);  new base A (2/4);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCA  >  NC_002947/6158838‑6159124                                                                                                                                                    |                                                                                                                                              aGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCtctc                                                                                                                                              >  4:308609/1‑148 (MQ=255)   ccGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAg                                                                                                                                           >  8:50632/1‑149 (MQ=255)   ccGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAg                                                                                                                                           <  5:157089/149‑1 (MQ=255)      cAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAggcg                                                                                                                                        <  6:297929/149‑1 (MQ=255)             tCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCttttt                                                                                                                                 <  7:15152/149‑1 (MQ=255)             tCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGAATTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCtttt                                                                                                                                  <  7:268879/148‑1 (MQ=255)              ccATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCtttttt                                                                                                                                >  8:33128/1‑149 (MQ=255)              ccATTTCGCTTTAGCGGCAAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCttttta                                                                                                                                <  7:33964/149‑2 (MQ=255)               cATTTCGCTTTTGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCttttta                                                                                                                                <  1:129351/148‑2 (MQ=255)               cATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCTTGAGAGGGCGGCttttta                                                                                                                                <  5:286723/148‑2 (MQ=255)                       tttAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑Cg                                                                                                                       >  3:208036/1‑148 (MQ=255)                          aGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTTAAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCACTTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGAt                                                                                                                    <  5:72051/148‑1 (MQ=255)                           gCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTc                                                                                                                  <  5:21007/149‑1 (MQ=255)                           gCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTc                                                                                                                  <  7:33128/149‑1 (MQ=255)                           gCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAAATCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTc                                                                                                                  <  5:313155/149‑1 (MQ=255)                            cGGCTAAATTTTTCGTTTGGTTGCGTCATTAATTTATGATTTACCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGATAAAAAGCCGCCATCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTc                                                                                                                  >  1:44306/1‑149 (MQ=255)                               cTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTa                                                                                                               >  1:317322/1‑148 (MQ=255)                                tAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATc                                                                                                             >  5:45200/1‑149 (MQ=255)                                tAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATc                                                                                                             <  6:45200/149‑1 (MQ=255)                                tAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATc                                                                                                             >  6:142955/1‑149 (MQ=255)                                  aaTTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTGAGTACGTCGCCGACGAACCTCTAAATTCGTTGGCTAAAAAGCCTCCATCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCgg                                                                                                           >  3:192508/1‑149 (MQ=255)                                  aaTTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCTCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGAGGCTTTTTTTTGCCT‑CGGATTCGTATCgg                                                                                                           <  7:3607/149‑1 (MQ=255)                                        tCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACcgt                                                                                                     >  7:267027/1‑149 (MQ=255)                                           tttGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACcgtcgt                                                                                                  <  2:99016/149‑1 (MQ=255)                                                ttGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGtgctgc                                                                                             <  4:208036/149‑1 (MQ=255)                                                           a‑ttaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGcc                                                                                   <  7:50632/148‑1 (MQ=255)                                                             ttaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGcc                                                                                            >  5:28444/1‑138 (MQ=255)                                                             ttaTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGcc                                                                                            <  6:28444/138‑1 (MQ=255)                                                              taTGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGaa                                                                                >  5:260932/1‑149 (MQ=255)                                                                   tttAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTg                                                                            <  3:60599/148‑1 (MQ=255)                                                                     tAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGAt                                                                          <  2:317322/148‑1 (MQ=255)                                                                        cTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACt                                                                      <  3:147200/149‑1 (MQ=255)                                                                        cTATAGTCCTCAACGTTCGACTGTTTCACTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACt                                                                      >  1:181881/1‑149 (MQ=255)                                                                         tataGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTATTCCA‑CGGATTAGTAACGGTACCGTCGTGCTGCCAGTCTCGCCGAACTTGATCACtt                                                                     >  7:104286/1‑148 (MQ=255)                                                                                  cAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAc                                                              >  6:94796/1‑147 (MQ=255)                                                                                   aaCGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTTGCCGAACTTGATCACTGTGAGGAc                                                              >  7:30877/1‑146 (MQ=255)                                                                                   aaCGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGGGAGGAc                                                              >  4:197741/1‑146 (MQ=255)                                                                                           aCTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCAt                                                   >  8:44663/1‑149 (MQ=255)                                                                                             tGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGc                                                 >  5:16172/1‑149 (MQ=255)                                                                                               tttCATTACGTTGCCGACCAACCTCTCAATTCGTTGGATAAAAAGCCGCCCTCTCAAGGCGGGTTTTTTTTGCCG‑CGGATTCGTATCGGTACCGTCGGGCGGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCg                                                <  2:166745/148‑1 (MQ=255)                                                                                                   aGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTaca                                            >  8:111027/1‑148 (MQ=255)                                                                                                       cGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCaa                                       <  3:219287/149‑1 (MQ=255)                                                                                                          cGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTTAAATCCGTACACACAaccc                                    >  8:277648/1‑148 (MQ=255)                                                                                                               acgaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCtgt                               <  7:111027/149‑1 (MQ=255)                                                                                                                 gaACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGa                             >  5:127746/1‑149 (MQ=255)                                                                                                                        tcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACtgt                                                                    >  8:138664/1‑103 (MQ=255)                                                                                                                        tcAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCg                      >  5:139858/1‑149 (MQ=255)                                                                                                                          aaTTCGTTGTCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCg                      <  2:44306/147‑1 (MQ=255)                                                                                                                          aaTTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTTGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCg                      <  8:30877/147‑1 (MQ=255)                                                                                                                          aaTTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGTACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGc                    <  6:260932/149‑1 (MQ=255)                                                                                                                           aTTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGc                                                 <  8:208259/119‑1 (MQ=255)                                                                                                                           aTTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGACGAACTTGATCACTGTGAGGACCAGATTATCATGc                                                 >  7:208259/1‑119 (MQ=255)                                                                                                                               gTTGGCTAAAAAGCCGCCCTCTCAAGCCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACAGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGtt               >  8:133418/1‑149 (MQ=255)                                                                                                                                  ggCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTgg             <  8:267027/148‑1 (MQ=255)                                                                                                                                   gCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGc            <  7:44663/148‑1 (MQ=255)                                                                                                                                   aaaaaaaaGCCGCCTTGAGAGGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGAACAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCa           >  5:149256/4‑149 (MQ=255)                                                                                                                                     tAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAAc         >  3:179299/1‑149 (MQ=255)                                                                                                                                      aaaaaGCCGCCTTGAGAGGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTCTGACTGTTCGGCCTGTTGGCAAc         >  6:286080/1‑148 (MQ=255)                                                                                                                                         aaGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCTGCGGATTCGTATCGGTACCGTCGTGGTGCCAGTGTCGCCGAACTTGATCACTGTGGGGCCCAGATTTTCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCAc      >  1:203488/1‑149 (MQ=255)                                                                                                                                            ccgccCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCa  <  5:94796/149‑1 (MQ=255)                                                                                                                                                   ctcAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTc                       >  8:125577/1‑121 (MQ=255)                                                                                                                                                   ctcAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTc                       <  7:125577/121‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              AGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCA  >  NC_002947/6158838‑6159124

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