Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,211,904	C→T	7.9%	D12D (GAC→GAT)	PP_1955 →	cytochrome P450 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,211,904	0	C	T	7.9%	179.0 / 3.9	63	D12D (GAC→GAT)	PP_1955	cytochrome P450 family protein
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base T (3/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.07e-01

TGGCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGG  >  NC_002947/2211757‑2212051                                                                                                                                                    |                                                                                                                                                     tggCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACt                                                                                                                                                     >  6:296491/1‑149 (MQ=255)   gCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACtt                                                                                                                                                    <  4:131585/148‑1 (MQ=255)   gCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCATCAAGCGCCATCGGACtt                                                                                                                                                    >  7:291029/1‑148 (MQ=255)     ctctGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCaa                                                                                                                                                 <  8:72737/149‑1 (MQ=255)     ctctGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCa                                                                                                                                                  >  3:30898/1‑148 (MQ=255)          cGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATc                                                                                                                                               <  8:213911/146‑1 (MQ=255)          cGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGa                                                                                                                                            <  4:30898/149‑1 (MQ=255)                  tAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGATTGAAGTCCGAtggcgc                                                                                                                                      <  8:44387/147‑5 (MQ=255)                      ttAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCaga                                                                                                                                <  6:61301/149‑1 (MQ=255)                           tataGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGtttt                                                                                                                           <  5:209378/149‑1 (MQ=255)                                aTTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGa                                                                                                                      >  3:137153/1‑149 (MQ=255)                                     cTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCAGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGCCTTCATACCGATGAGTGAGGAGAGTTTTCGTGATCCAg                                                                                                                 >  8:131174/1‑149 (MQ=255)                                          cGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTc                                                                                                             >  6:136569/1‑148 (MQ=255)                                            gcATGTAGTGATACTTCCTTGGTGCGTAGCCGCGGTCCGTCTGCAATAATTGCTCGGAGAATACCAAAAATGTAACTCCTTGATCGTCATCACCCGCCATCGCACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAAt                                                                                                          <  4:21741/148‑1 (MQ=255)                                                gTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGc                                                                                                      >  8:133658/1‑149 (MQ=255)                                                   gtcaTACTTACTTGGTGTGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGCAAATACTTGATCGTCCTCAAGCGCCATCGGACATCAATCCGATGAGTGAGCAGAGTTTACGTGATCCAGCCGCAATCTGCCAg                                                                                                   >  6:16857/4‑149 (MQ=255)                                                       tacttactTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTg                                                                                               >  3:103566/1‑149 (MQ=255)                                                       tacttactTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTg                                                                                               >  8:133298/1‑149 (MQ=255)                                                         cttacttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTg‑c                                                                                             <  5:289809/148‑1 (MQ=255)                                                          ttacttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAtgagt                                                                                                                                       <  3:239680/106‑1 (MQ=255)                                                          ttacttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAtgagt                                                                                                                                       >  4:239680/1‑106 (MQ=255)                                                             cttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATCGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTg‑cgcgc                                                                                         <  2:308923/148‑1 (MQ=255)                                                                ggtggGTTGCCGCGTGCCTGAGGATATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATGGGTCTTCAATCCGATGAGTGAGCTGAGTTTTCGTGATCCAGGCTCAATCTGCCAGCGTG‑CGCGCgaa                                                                                      <  5:16857/148‑1 (MQ=255)                                                                      tAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCACCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCAAAGAAAccc                                                                               >  3:63462/1‑149 (MQ=255)                                                                                ccTGAGGAAATAATTGTTCGTAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCttttt                                                                     <  8:291029/149‑1 (MQ=255)                                                                                ccTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCACCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCAAAGAAACCCCTGTCttttt                                                                     <  4:63462/149‑1 (MQ=255)                                                                                      gAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATgc                                                               >  8:189075/1‑149 (MQ=255)                                                                                         ataatTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGcgcc                                                            >  8:198478/1‑149 (MQ=255)                                                                                          taatTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGcgccg                                                           <  4:137153/149‑1 (MQ=255)                                                                                           aatTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGcgccgc                                                          >  8:184190/1‑149 (MQ=255)                                                                                            atTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGcgccgc                                                          <  2:199343/148‑1 (MQ=255)                                                                                            atTGTTCGGAGAATAACAAAAATGGAAATACATGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGcgccgc                                                          >  1:199343/1‑148 (MQ=255)                                                                                                  cGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTg‑cg                                                                                            >  3:69844/1‑108 (MQ=255)                                                                                                  cGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTg‑cg                                                                                            <  4:69844/108‑1 (MQ=255)                                                                                                   ggAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTACGCGCCGCTAGGCgt                                                   <  4:102852/148‑1 (MQ=255)                                                                                                     agaATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTGCCGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCgtgt                                                 <  1:230949/149‑1 (MQ=255)                                                                                                        aTAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGAt                                             <  7:184190/149‑1 (MQ=255)                                                                                                             aaaaaTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTaa                                        >  2:8554/1‑149 (MQ=255)                                                                                                             aaaaaTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTaa                                        <  5:296491/149‑1 (MQ=255)                                                                                                             aaaaaTGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTaa                                        >  4:319265/1‑149 (MQ=255)                                                                                                                   ggAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGc                                                     <  1:200091/130‑1 (MQ=255)                                                                                                                   ggAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGc                                                     >  2:200091/1‑130 (MQ=255)                                                                                                                     aaaTACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGAGTc                                  >  7:174905/1‑147 (MQ=255)                                                                                                                      aaTACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTc                                  >  6:202188/1‑146 (MQ=255)                                                                                                                      aaTACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTc                                  >  3:230232/1‑146 (MQ=255)                                                                                                                          cTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGagaaga                           <  1:8554/149‑1 (MQ=255)                                                                                                                          cTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGagaag                            <  7:133658/148‑1 (MQ=255)                                                                                                                             gATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGc                        >  8:123416/1‑149 (MQ=255)                                                                                                                                 gTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGa                     <  4:230232/148‑1 (MQ=255)                                                                                                                                 gTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGa                     <  7:133298/148‑1 (MQ=255)                                                                                                                                    cTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAgcgc                 >  7:261391/1‑149 (MQ=255)                                                                                                                                    cTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAgcgc                 >  1:115556/1‑149 (MQ=255)                                                                                                                                      cAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGATGATGCCGAGCGCGt               >  8:140000/1‑149 (MQ=255)                                                                                                                                       aaGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtt              >  3:274602/1‑149 (MQ=255)                                                                                                                                        tcactcATCGGATTGAAGTCCGATGGGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTc                                  >  7:137822/6‑128 (MQ=255)                                                                                                                                        tcactcATCGGATTGAAGTCCGATGGGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTc                                  <  8:137822/123‑1 (MQ=255)                                                                                                                                         gcgcCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt            <  5:136569/149‑1 (MQ=255)                                                                                                                                         gcgcCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGtttt            <  4:103566/149‑1 (MQ=255)                                                                                                                                            tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAg         >  6:125464/2‑149 (MQ=255)                                                                                                                                            tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAg         >  1:5209/2‑149 (MQ=255)                                                                                                                                                cGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAg     >  5:41010/1‑149 (MQ=255)                                                                                                                                                cGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAg     >  4:260013/1‑149 (MQ=255)                                                                                                                                                   aCTTCAATCCGATGAGTGAGCAGAGTGTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGATAAGATGCCGAGCGCGTTTTGAGCGAGTgg  <  7:131174/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGGCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGG  >  NC_002947/2211757‑2212051

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