Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,970,884	A→C	8.7%	S228A (TCG→GCG)	PP_4379 ←	3‑oxoacyl‑ACP synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,970,884	0	A	C	8.7%	165.6 / 4.6	56	S228A (TCG→GCG)	PP_4379	3‑oxoacyl‑ACP synthase
Reads supporting (aligned to +/- strand):  ref base A (23/28);  new base C (3/2);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAACACTGGAGGACACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAGGTGCTCGGCACCGGAACCGTCAGTGCCGAAACGCGC  >  NC_002947/4970741‑4971031                                                                                                                                                 |                                                                                                                                                    cGGAACACTGGAGGACACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCa                                                                                                                                                <  1:160144/149‑1 (MQ=255)         tGGAGGACACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCTGAAGTCGTCGGATAT‑CGACGCCTTctgc                                                                                                                                              <  6:235986/143‑4 (MQ=255)         tGGAGGACACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCTGAAGTCGTCGGATAT‑CGACGCCTTctgc                                                                                                                                              >  5:235986/1‑140 (MQ=255)         tGGAGGACACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTc                                                                                                                                                 <  4:93010/140‑1 (MQ=255)         tGGAGGACACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTc                                                                                                                                                 >  3:93010/1‑140 (MQ=255)           gcggaCTCCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCGTCAAAGCGTCGCGCCATGGCGTCGACAATCGCCGCACTGCCCTGGCGGATGCAGCAGGCGTC‑GATATCAGACGCCTTCAGGCCCGTGGc                                                                                                                                      >  7:34645/3‑149 (MQ=255)           gaggaCACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGc                                                                                                                                      <  6:356964/149‑1 (MQ=255)              tacacCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATc                                                                                                                                   <  3:208652/148‑1 (MQ=255)              gaCACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATc                                                                                                                                   <  3:243976/149‑1 (MQ=255)                cacCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATc                                                                                                                                   >  4:208652/1‑147 (MQ=255)                cacCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATc                                                                                                                                   >  3:9873/1‑147 (MQ=255)                   cGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCcag                                                                                                                                >  5:264973/1‑147 (MQ=255)                       ttGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTc                                                                                                                          <  2:349421/149‑1 (MQ=255)                            ggTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCa                                                                                                                                                >  6:142085/1‑122 (MQ=255)                            ggTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCa                                                                                                                                                <  5:142085/122‑1 (MQ=255)                                         tCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCAc                                                                                                        <  3:318241/149‑1 (MQ=255)                                            ttGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCtt                                                                                                     >  1:2583/1‑149 (MQ=255)                                               aCGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGAt                                                                                                  >  3:5496/1‑149 (MQ=255)                                               aCGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGAt                                                                                                  >  3:45189/1‑149 (MQ=255)                                                 gAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGt                                                                                                                  <  4:220200/131‑1 (MQ=255)                                                   aCTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAgg                                                                                              <  2:2583/149‑1 (MQ=255)                                                        tcCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGaa                                                                                         <  4:9873/149‑1 (MQ=255)                                                          cGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGt                                                                                       <  6:264973/149‑1 (MQ=255)                                                                ttGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATa                                                                                 <  7:67923/149‑1 (MQ=255)                                                                   ctctTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTNCGCt                                                                                                             <  5:154359/118‑1 (MQ=255)                                                                   ctctTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCt                                                                                                             >  6:154359/1‑118 (MQ=255)                                                                    tctTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTg                                                                             >  6:23672/1‑149 (MQ=255)                                                                      ttCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGaa                                                                                   >  2:311840/1‑141 (MQ=255)                                                                      ttCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGaa                                                                                   <  1:311840/141‑1 (MQ=255)                                                                       tCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAAt                                                                                  >  3:86783/1‑141 (MQ=255)                                                                       tCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAAt                                                                                  <  4:86783/141‑1 (MQ=255)                                                                         aaaGCGTCGCGCCACGGCGTCGACAAACGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTggcggc                                                                        >  2:205625/1‑149 (MQ=255)                                                                            gCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGTGAAGTTGAATACCTGGCGGCCGt                                                                     <  4:45189/149‑1 (MQ=255)                                                                                   gcCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGGTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCa                                                                  >  7:185202/1‑145 (MQ=255)                                                                                   gcCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCa                                                                  <  8:185202/145‑1 (MQ=255)                                                                                   gcCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATgaa                                                              <  2:277593/149‑1 (MQ=255)                                                                                   gcCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATg                                                                >  1:277593/1‑147 (MQ=255)                                                                                          gcgcGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTg                                                       <  3:159960/146‑1 (MQ=255)                                                                                            tCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGcc                                                     <  4:5496/149‑1 (MQ=255)                                                                                              gACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGt                                                   >  2:283425/1‑149 (MQ=255)                                                                                              gACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGt                                                   >  6:202041/1‑149 (MQ=255)                                                                                                     gccgcACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGa                                             <  3:157419/148‑1 (MQ=255)                                                                                                       cgcACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACtt                                          >  4:333346/1‑149 (MQ=255)                                                                                                         cACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATgaagaa                                                           <  7:324223/130‑1 (MQ=255)                                                                                                         cACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATgaagaa                                                           >  8:324223/1‑130 (MQ=255)                                                                                                         cACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACttt                                         >  8:277057/1‑148 (MQ=255)                                                                                                           cTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACAGACTTTCAg                                      <  5:23672/149‑1 (MQ=255)                                                                                                              cccTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAg                                      >  5:45427/1‑146 (MQ=255)                                                                                                              cccTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAg                                      <  6:45427/146‑1 (MQ=255)                                                                                                                   gtggATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAGGTGCTc                                <  3:3463/147‑1 (MQ=255)                                                                                                                   gtggATGCAGAAGGCATC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAGGTGCTCgg                              >  2:144249/1‑149 (MQ=255)                                                                                                                        tGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAGGTGCTCGGCACCg                         <  3:333346/149‑1 (MQ=255)                                                                                                                           tgAAGTCGTCGGATAT‑CGACGCCTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCg                                                    <  5:260599/118‑1 (MQ=255)                                                                                                                           tgAAGTCGTCGGATAT‑CGACGCCTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCg                                                    >  6:260599/2‑119 (MQ=255)                                                                                                                                  gTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCg                                                    <  4:174716/112‑1 (MQ=255)                                                                                                                                  gTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCg                                                    >  3:174716/1‑112 (MQ=255)                                                                                                                                            cgacgaCTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAg                                      <  6:171364/117‑1 (MQ=255)                                                                                                                                            cgacgaCTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAg                                      >  5:171364/1‑117 (MQ=255)                                                                                                                                                gacTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAGGTGCTCGGCACCGGAAcc                    <  2:164508/131‑1 (MQ=255)                                                                                                                                                gacTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAGGTGCTCGGCACCGGAAcc                    >  1:164508/1‑131 (MQ=255)                                                                                                                                                gacTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAGGTGCTCGGCACCGGAACCGTCAGTGCCGAAACgcgc  >  4:181957/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                    CGGAACACTGGAGGACACCGTATTGCCGGTTTCCAGCATGTCCTTGACGAACTTCTCCGGGTGTTGCTCTTCAAAGCGTCGCGCCACGGCGTCGACAATCGCCGCACTGCCCTGGTGGATGCAGAAGGCGTC‑GATATCCGACGACTTCAGGCCGCTGGCATCCAGCAGCTCGTGCAGGTGCGCTGGCACCTTGATCAGGGCGAAGTTGAATACCTGGCGGCCGTTCATGAAGAACTTGCCGTCACTGACTTTCAGGTGCTCGGCACCGGAACCGTCAGTGCCGAAACGCGC  >  NC_002947/4970741‑4971031

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