Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	6,016,729	C→A	6.7%	intergenic (‑45/+182)	PP_5728 ← / ← PP_5267	rhs repeat‑associated core domain‑containing protein/cytochrome c5

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	6,016,729	0	C	A	6.7%	189.0 / 3.7	60	intergenic (‑45/+182)	PP_5728/PP_5267	rhs repeat‑associated core domain‑containing protein/cytochrome c5
Reads supporting (aligned to +/- strand):  ref base C (30/26);  new base A (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.88e-01

AATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGC  >  NC_002947/6016588‑6016875                                                                                                                                              |                                                                                                                                                   aaTAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGGGTTCTGAAGAACACTGacgc                                                                                                                                             <  5:34786/149‑4 (MQ=255)     ggCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATa                                                                                                                                          <  6:313438/148‑1 (MQ=255)      gCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAg                                                                                                                                         >  3:63465/1‑148 (MQ=255)       ccAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACTTATTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGaa                                                                                                                                       <  1:295888/149‑1 (MQ=255)          gATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGAGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTg                                                                                                                                    <  5:285320/149‑1 (MQ=255)           aTGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTgg                                                                                                                                   >  8:269407/1‑149 (MQ=255)                gAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTaaa                                                                                                                               >  8:83738/1‑148 (MQ=255)                      cgcTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTTGTAAAAATCCTa                                                                                                                        >  4:199401/1‑149 (MQ=255)                      cgcTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTa                                                                                                                        <  8:304151/149‑1 (MQ=255)                      cgcTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTa                                                                                                                        <  2:112988/149‑1 (MQ=255)                      cgcTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTa                                                                                                                        <  7:83738/149‑1 (MQ=255)                           ccGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGtt                                                                                                                     <  4:63465/147‑1 (MQ=255)                             ggggATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCa                                                                                                                 >  3:340547/1‑149 (MQ=255)                              gggATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCAt                                                                                                                <  7:137407/149‑1 (MQ=255)                                gATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCt                                                                                                                         >  8:186880/1‑138 (MQ=255)                                gATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCt                                                                                                                         <  7:186880/138‑1 (MQ=255)                                     gCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAAt                                                                                                                            >  7:301293/1‑130 (MQ=255)                                     gCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAAt                                                                                                                            <  8:301293/130‑1 (MQ=255)                                         ttGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACATTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAgg                                                                                                       <  3:76282/147‑1 (MQ=255)                                             ggCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCAACACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGc                                                                                                 >  8:287802/1‑149 (MQ=255)                                                tgtgTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCCCAGccc                                                                                               >  2:164918/1‑148 (MQ=255)                                                 gtgtTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGa                                                                                             >  2:23557/1‑149 (MQ=255)                                                  tgtTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGa                                                                                             <  3:129854/148‑1 (MQ=255)                                                   gtTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGAtt                                                                                           >  1:47037/1‑149 (MQ=255)                                                           ggCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTc                                                                                   <  4:340547/149‑1 (MQ=255)                                                             cTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAACCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTg                                                                                 >  5:52832/1‑149 (MQ=255)                                                                    ttGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACAt                                                                          <  2:47037/149‑1 (MQ=255)                                                                      gCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATc                                                                         >  3:58513/1‑148 (MQ=255)                                                                       cTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCa                                                                       >  2:339431/1‑149 (MQ=255)                                                                        tCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGgag                                                                              <  3:76230/141‑1 (MQ=255)                                                                        tCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGgag                                                                              >  4:76230/1‑141 (MQ=255)                                                                           ggACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGa                                                                   <  7:269407/149‑1 (MQ=255)                                                                           ggACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCCCAGCCCGATTGGAGCTTCTGGAGACATCCATGGa                                                                   <  1:164918/149‑1 (MQ=255)                                                                            gACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGa                                                                                        >  6:66167/1‑127 (MQ=255)                                                                            gACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGa                                                                                        <  5:66167/127‑1 (MQ=255)                                                                             aCAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACa                                                                 >  5:278850/1‑149 (MQ=255)                                                                               aGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATg                                                                                                               >  4:284639/1‑101 (MQ=255)                                                                               aGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATg                                                                                                               <  3:284639/101‑1 (MQ=255)                                                                               aGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATg                                                               >  5:277809/1‑149 (MQ=255)                                                                                      gCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAg                                                        >  2:62558/1‑149 (MQ=255)                                                                                       gccAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGc                                                       <  7:10391/148‑1 (MQ=255)                                                                                       cccAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGc                                                       >  6:86468/1‑149 (MQ=255)                                                                                              ccATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCt                                                                                  >  4:72163/1‑115 (MQ=255)                                                                                              ccATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCt                                                                                  <  3:72163/115‑1 (MQ=255)                                                                                               cATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCt                                                                                                                         <  4:212326/75‑1 (MQ=255)                                                                                               cATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCt                                                                                                                         >  3:212326/1‑75 (MQ=255)                                                                                               cATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGATATCCATGGACATGGCTCTAGCGggttggt                                               >  2:85786/1‑149 (MQ=255)                                                                                               cATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGggttggt                                               >  8:31851/1‑149 (MQ=255)                                                                                               cATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGggttgg                                                <  7:287802/148‑1 (MQ=255)                                                                                                aTTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTg                                                                                          >  4:193497/1‑105 (MQ=255)                                                                                                   gTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGc                                            <  7:31851/148‑1 (MQ=255)                                                                                                    ttAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCt                                          >  3:239454/1‑149 (MQ=255)                                                                                                                          ttttCAGGGTTCTGAAGAACACTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGcac                    >  5:316046/4‑149 (MQ=255)                                                                                                                           cGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGcacc                   <  1:23557/149‑1 (MQ=255)                                                                                                                               aGTGTTCTTCAGAACCCTGAAAAATAGAACTGGAAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACaa               <  1:62558/149‑1 (MQ=255)                                                                                                                                gtgtTCATCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGCCATCCATGGACATGGCTCTAGCGCGTTGGTTGCCTGTGCTGGACATTTCTCGGGCAcaacaa               >  7:169011/1‑148 (MQ=255)                                                                                                                                gggtTCTGAAGAACACTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACaa               >  4:325172/3‑148 (MQ=255)                                                                                                                                       tcAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAg       >  8:249214/1‑149 (MQ=255)                                                                                                                                        aaGAACACTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACaaa        <  5:43265/146‑1 (MQ=255)                                                                                                                                          gAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGccc    >  3:119096/1‑149 (MQ=255)                                                                                                                                             cccTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGc  >  2:192167/1‑148 (MQ=255)                                                                                                                                              |                                                                                                                                                   AATAGGCCAGATGCCGAAGGGCGCTGCCGGGGATCAGCTTTTGCGGCTGTGTTCGACCGGCTGCATGTTGCTCGGGACAGTTATCGCCCAGATCCATTGTTAGCCACCACTTCAATCAGAAGCGTCAGTGTTCTTCAGAACCCTGAAAAATAGAACTGGTAAAAATGCTAGTTCTCATGGAATCAGGCTCAGCCCGATTGGAGCTTCTGGAGACATCCATGGACATGGCTCTAGCGGGTTGGTTGCCTGTGCTGGACCTTTCTCGGGCACCACAAAAACAAAGCCCGC  >  NC_002947/6016588‑6016875

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