Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	508,221	T→C	6.9%	N264D (AAC→GAC) ‡	estP ←	esterase EstP

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	508,221	0	T	C	6.9%	163.6 / 2.9	58	N264D (AAC→GAC) ‡	estP	esterase EstP
Reads supporting (aligned to +/- strand):  ref base T (34/20);  new base C (2/2);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.34e-01

TAATCGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGTGG  >  NC_002947/508076‑508368                                                                                                                                                   |                                                                                                                                                     tAATCGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTg                                                                                                                                                   >  6:202858/1‑149 (MQ=255)     cGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCg                                                                                                                                               >  5:216476/1‑149 (MQ=255)        cAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGACCTGGTCggca                                                                                                                                            <  5:260778/149‑4 (MQ=255)         aaTCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTgg                                                                                                                                           >  6:310055/1‑149 (MQ=255)          aTCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGc                                                                                                                                          <  5:13386/149‑1 (MQ=255)              gctgctGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAgg                                                                                                                                      <  6:216476/149‑1 (MQ=255)                      cccGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAact                                                                                                                               >  3:228526/1‑148 (MQ=255)                                   ggTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGc                                                                                                                  >  1:6483/1‑148 (MQ=255)                                   ggTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGc                                                                                                                  >  3:240019/1‑148 (MQ=255)                                       gACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAgc                                                                                                              <  5:202858/148‑1 (MQ=255)                                        aCCGAGTCGTTGAACAGCAGTGTGGTCGGGTCTGGCGTTGTGCCGTTGATCCAGTACACCGGGTTTTACACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAgcgc                                                                                                            <  4:188501/149‑1 (MQ=255)                                          cGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCt                                                                                                          >  2:239584/1‑149 (MQ=255)                                            aGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCg                                                                                                        >  2:116080/1‑149 (MQ=255)                                                    aaCAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACa                                                                                                 <  5:310055/148‑1 (MQ=255)                                                    aaCAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACa                                                                                                 <  2:6483/148‑1 (MQ=255)                                                     aCAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTgggg                                                                                                                       <  7:103504/125‑1 (MQ=255)                                                     aCAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTgggg                                                                                                                       >  8:103504/1‑125 (MQ=255)                                                           gTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAgg                                                                                          >  7:79811/1‑148 (MQ=255)                                                           gTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAgg                                                                                          >  2:20978/1‑148 (MQ=255)                                                            tttGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCATCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTg                                                                                                                                            >  4:174799/1‑97 (MQ=255)                                                            tttGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGa                                                                                        >  5:8675/1‑149 (MQ=255)                                                            tttGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGa                                                                                        >  1:228447/1‑149 (MQ=255)                                                                 tCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGtt                                                                                    >  5:182973/1‑148 (MQ=255)                                                                      tCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTg                                                                               >  8:151948/1‑148 (MQ=255)                                                                       cTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTCGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGa                                                                             <  1:239584/149‑1 (MQ=255)                                                                          gCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGa                                                                             <  4:240019/146‑1 (MQ=255)                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGc                                                                                                                                          >  6:84243/1‑70 (MQ=255)                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGc                                                                                                                                          <  5:84243/70‑1 (MQ=255)                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATc                                                          >  8:190845/1‑149 (MQ=255)                                                                                         tCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATc                                                          >  8:269975/1‑149 (MQ=255)                                                                                           ccGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGTGCGCTGGCCAGGCCGAACTGACGGGGGCGGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCt                                                         <  8:305412/148‑1 (MQ=255)                                                                                           ccGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTg                                                        >  4:293377/1‑149 (MQ=255)                                                                                                cacCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGAt                                                                                       >  4:291906/1‑114 (MQ=255)                                                                                                cacCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGAt                                                                                       <  3:291906/114‑1 (MQ=255)                                                                                                   cGGGTTTTCCACGCAGCTATCGCCGCTATGAGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAgc                                                 >  1:198880/1‑149 (MQ=255)                                                                                                   cGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAgct                                                >  6:159171/1‑149 (MQ=255)                                                                                                         ttCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGa                                          >  6:267860/1‑149 (MQ=255)                                                                                                            cACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCa                                       <  7:151948/149‑1 (MQ=255)                                                                                                              cacAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGt                                     <  7:20054/147‑1 (MQ=255)                                                                                                                 aGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATTTGCCCTAGCTGGCTGATCAGTGa                                   <  5:293766/148‑1 (MQ=255)                                                                                                                 aGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATTTGCCCTAGCTGGCTGATCAGTGa                                   >  6:293766/1‑148 (MQ=255)                                                                                                                        gccgcTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTg                            >  3:323303/1‑148 (MQ=255)                                                                                                                         ccgcTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGTCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACATCAAAGGGATGTTCAGt                                                                                <  3:312552/96‑1 (MQ=255)                                                                                                                         ccgcTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCAAGGTCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACATCAAAGGGATGTTCAGt                                                                                >  4:312552/1‑96 (MQ=255)                                                                                                                                    aGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAACTGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGga                 <  1:97740/149‑1 (MQ=255)                                                                                                                                        agcGACCAGAACCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGa                                          <  4:288446/117‑1 (MQ=255)                                                                                                                                        agcGACCAGAACCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGa                                          >  3:288446/3‑119 (MQ=255)                                                                                                                                        agcGACCAGAACCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGagag              >  1:33779/3‑147 (MQ=255)                                                                                                                                         ccGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCa            <  1:116080/148‑1 (MQ=255)                                                                                                                                         ccGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCa            <  7:190845/148‑1 (MQ=255)                                                                                                                                         ccGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAg           >  8:40047/1‑149 (MQ=255)                                                                                                                                           gACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGt          >  7:214503/1‑148 (MQ=255)                                                                                                                                            aCCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGa        >  7:98855/1‑149 (MQ=255)                                                                                                                                            aCCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGa        <  4:228526/149‑1 (MQ=255)                                                                                                                                              cAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGa        >  6:342589/1‑147 (MQ=255)                                                                                                                                              cAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGa        >  1:299764/1‑147 (MQ=255)                                                                                                                                                 gTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGCGAGCTGTGACAggg   >  4:203611/1‑147 (MQ=255)                                                                                                                                                 gTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGTg   >  5:128551/1‑149 (MQ=255)                                                                                                                                                  ttCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGTgg  <  6:8675/149‑1 (MQ=255)                                                                                                                                                  ttCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGTgg  <  2:198880/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     TAATCGGCAATCAGCTGCTGACCCGCGATGGTCGGGTGGACCGAGTCGTTGAACAGCAGTTTGGTCGGGTCTGGCGTTGTGCCGTTGATCCCGTACACCGGGTTTTCCACGCAGCTATCGCCGCTAT‑AGCAGGTGCCGACCAGGTTCTGGTCGCTGGCCAGGCCGAACTGACTGGGGCTGGCCAGCGCCTCGCTCAACAGCAAAGGGATGTTCAGTGGAA‑TGATCTGGGCATCGATCTGCCCTAGCTGGCTGATCAGTGACTGGTTGAACGCAGCGGAGAGCAGTGACAGTGG  >  NC_002947/508076‑508368

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