Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,158,981	C→T	10.2%	intergenic (+318/+550)	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,981	0	C	T	10.2%	155.1 / 6.6	60	intergenic (+318/+550)	PP_5739/PP_5402	hypothetical protein/NAD‑dependent acetyl‑CoA synthetase deacetylase subunit
Reads supporting (aligned to +/- strand):  ref base C (26/26);  new base T (2/4);  total (30/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 = 3.46e-01

ACCAAGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCA  >  NC_002947/6158834‑6159124                                                                                                                                                     |                                                                                                                                                 aCCAAGCCGCAAGTGATGCATTTCGCTTTAGCGGCTAAATTTTTCGTTTTGTTGCGTCATTAA‑TTATTATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCct                                                                                                                                                 <  2:237136/149‑1 (MQ=255)     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                                                                                                                  <  7:33128/149‑1 (MQ=255)                               gCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTc                                                                                                                  <  5:21007/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                                                                                                             >  6:142955/1‑149 (MQ=255)                                    tAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATc                                                                                                             <  6:45200/149‑1 (MQ=255)                                    tAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATc                                                                                                             >  5:45200/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)                                                                                                                                                     |                                                                                                                                                 ACCAAGCCGCAAGTGATCCATTTCGCTTTAGCGGCTAAATTTTTCGTTTGGTTGCGTCATTAA‑TTATGATTTAGCTATAGTCCTCAACGTTCGACTGTTTCAGTACGTCGCCGACGAACCTCTCAATTCGTTGGCTAAAAAGCCGCCCTCTCAAGGCGGCTTTTTTTTGCCT‑CGGATTCGTATCGGTACCGTCGTGCTGCCAGTGTCGCCGAACTTGATCACTGTGAGGACCAGATTATCATGCGTACACTCAACCACCTGTGACTGTTCGGCCTGTTGGCAACCACGCCA  >  NC_002947/6158834‑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: ‑