Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,211,906	T→G	8.3%	F13C (TTC→TGC) ‡	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,906	0	T	G	8.3%	179.7 / 7.2	60	F13C (TTC→TGC) ‡	PP_1955	cytochrome P450 family protein
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base G (3/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.93e-01

GCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGGG  >  NC_002947/2211759‑2212052                                                                                                                                                    |                                                                                                                                                    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                                                                                                                                        >  4:239680/1‑106 (MQ=255)                                                        ttacttGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGAtgagt                                                                                                                                        <  3:239680/106‑1 (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                                                                                             <  4:69844/108‑1 (MQ=255)                                                                                                cGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTg‑cg                                                                                             >  3:69844/1‑108 (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                  >  1:115556/1‑149 (MQ=255)                                                                                                                                  cTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAgcgc                  >  7:261391/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          >  1:5209/2‑149 (MQ=255)                                                                                                                                          tcATCGGATTGAAGTCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAg          >  6:125464/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)                                                                                                                                                   ttCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTggg  >  7:243906/1‑148 (MQ=255)                                                                                                                                                   ttCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTggg  >  7:160589/1‑148 (MQ=255)                                                                                                                                                   ttCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTggg  >  6:313324/1‑148 (MQ=255)                                                                                                                                                   ttCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTATTCTATGCGCCGCTAGGCGAGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTggg  >  1:103489/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GCCTCTGCGGCACTCTAAGTTAAATATAGATTTCCTCAGCGTCATGTAGTGATACTTACTTGGTGGGTAGCCGCGTGCCTGAGGAAATAATTGTTCGGAGAATAACAAAAATGGAAATACTTGATCGTCCTCAAGCGCCATCGGACTTCAATCCGATGAGTGAGCAGAGTTTTCGTGATCCAGCCTCAATCTGCCAGCGTG‑CGCGCGAAGAAACCCCTGTCTTTTTCTATGCGCCGCTAGGCGTGTGGATGGTAACGCGTCGAGAAGATGCCGAGCGCGTTTTGAGCGAGTGGG  >  NC_002947/2211759‑2212052

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