Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	498,050	T→A	8.0%	L640Q (CTG→CAG)	PP_0409 →	two‑component system sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	498,050	0	T	A	8.0%	184.5 / 4.5	62	L640Q (CTG→CAG)	PP_0409	two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base T (25/32);  new base A (2/3);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

CCTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCT  >  NC_002947/497904‑498196                                                                                                                                                   |                                                                                                                                                   ccTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGCGCTGGCCTAGGCCGAGCTGGACTCTTGCCTGAATGATCGCCTGGc                                                                                                                                                 <  2:897412/150‑1 (MQ=255) ccTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGc                                                                                                                                                 <  1:1126275/150‑1 (MQ=255) ccTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGc                                                                                                                                                 <  1:941180/150‑1 (MQ=255) ccTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCCGAATGATCGCCTGGc                                                                                                                                                 <  1:1058409/150‑1 (MQ=255)      agagCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCa                                                                                                                                            <  1:149938/150‑1 (MQ=255)          cgtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAttca                                                                                                                                        <  2:448339/150‑4 (MQ=255)           gtgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAttcag                                                                                                                                       >  2:537508/1‑146 (MQ=255)            tgcgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAttcagg                                                                                                                                      <  1:537508/150‑6 (MQ=255)              cgtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGc                                                                                                                                    <  2:872106/150‑1 (MQ=255)              cgtgAAGAGGAAAAAGCCCGTATCGACCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCGTGAGCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGc                                                                                                                                    <  1:210376/150‑1 (MQ=255)               gtgAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGcc                                                                                                                                   >  1:525011/1‑150 (MQ=255)                  aaGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCa                                                                                                                                            <  2:470332/138‑1 (MQ=255)                  aaGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCa                                                                                                                                            >  1:470332/1‑138 (MQ=255)                               ccGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAgttgtt                                                                                                                   <  1:776818/150‑1 (MQ=255)                                cGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTc                                                                                                                  <  1:810363/150‑1 (MQ=255)                                cGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTc                                                                                                                  <  1:75641/150‑1 (MQ=255)                                 gTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTcc                                                                                                                 >  1:356080/1‑150 (MQ=255)                                    tCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCAGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTCCACTTCATGAAACGCCTGATCGCCCAGTTGTTCCAGt                                                                                                              <  1:1001739/150‑1 (MQ=255)                                     cGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGtt                                                                                                             <  2:424301/150‑1 (MQ=255)                                       cccGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTgg                                                                                                           <  2:1173853/150‑1 (MQ=255)                                             aaGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAGGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCt                                                                                                                                  >  2:643666/1‑121 (MQ=255)                                             aaGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAGGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCt                                                                                                                                  <  1:643666/121‑1 (MQ=255)                                                tGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAtt                                                                                                                                          >  1:746759/1‑109 (MQ=255)                                                tGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGCTCGCCAGGCGATCAtt                                                                                                                                          <  2:746759/110‑2 (MQ=255)                                                    tgatgaACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTg                                                                                                            <  1:462563/136‑1 (MQ=255)                                                    tgatgaACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTg                                                                                                            >  2:462563/1‑136 (MQ=255)                                                     gatgaACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGcc                                                                                             >  1:730534/1‑150 (MQ=255)                                                          aCTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGc                                                                                        <  1:1015398/150‑1 (MQ=255)                                                          aCTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGc                                                                                        <  2:356080/150‑1 (MQ=255)                                                          aCTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGc                                                                                        >  2:819964/1‑150 (MQ=255)                                                           cTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGcc                                                                                       <  1:819964/150‑1 (MQ=255)                                                            tGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCt                                                                                      >  2:1186517/1‑150 (MQ=255)                                                                 cAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGt                                                                                 <  2:1003113/150‑1 (MQ=255)                                                                     tGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTccgc                                                                             <  1:1024994/150‑1 (MQ=255)                                                                      gCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTccgcc                                                                            >  1:661390/1‑150 (MQ=255)                                                                        tGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGa                                                                          >  1:16306/1‑150 (MQ=255)                                                                        tGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGa                                                                          >  2:241923/1‑150 (MQ=255)                                                                        tGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGa                                                                          <  1:1006681/150‑1 (MQ=255)                                                                                     gCTTGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCgg                                                             >  2:952951/1‑150 (MQ=255)                                                                                      cTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTc                                                                                                <  2:295754/114‑1 (MQ=255)                                                                                      cTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTc                                                                                                >  1:295754/1‑114 (MQ=255)                                                                                        ggAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATCCCGGCAt                                                          >  2:685578/1‑150 (MQ=255)                                                                                               tCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCg                                                   <  2:525011/150‑1 (MQ=255)                                                                                                             ggtcTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCa                                     <  1:952951/147‑1 (MQ=255)                                                                                                              gCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTg                                                                                    >  2:913550/1‑102 (MQ=255)                                                                                                              gCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTg                                                                                    <  1:913550/102‑1 (MQ=255)                                                                                                                    gCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTACCGGGCGTCCGCCGATCCCCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCgcgg                              >  1:317754/1‑150 (MQ=255)                                                                                                                            ggACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGa                      <  2:16306/150‑1 (MQ=255)                                                                                                                             gACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTccgc                                                                             <  2:280485/94‑1 (MQ=255)                                                                                                                             gACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCATGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTccgc                                                                             >  1:280485/1‑94 (MQ=255)                                                                                                                               cccTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGc                   <  2:730534/150‑1 (MQ=255)                                                                                                                                 cTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCATTGCGTCAGTGGCTCCGCCAGGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCAc                 >  1:142274/1‑150 (MQ=255)                                                                                                                                      cTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAg                     <  1:215172/141‑1 (MQ=255)                                                                                                                                      cTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAg                     >  2:215172/1‑141 (MQ=255)                                                                                                                                      cTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCAGCAGATCAACAAATCATTCACATCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACg            >  1:444969/1‑150 (MQ=255)                                                                                                                                       tGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGGCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGc           >  2:230475/1‑150 (MQ=255)                                                                                                                                         aaTGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGGACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCaa         >  1:936763/1‑150 (MQ=255)                                                                                                                                           tGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGc                                                                <  2:255189/93‑1 (MQ=255)                                                                                                                                           tGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGc                                                                >  1:255189/1‑93 (MQ=255)                                                                                                                                              tCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATccc    >  2:1134085/1‑150 (MQ=255)                                                                                                                                                gCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCt  <  2:661390/150‑1 (MQ=255)                                                                                                                                                gCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCAGAGCCTCGATGCCGGCATTGCCTAGCCCGGCGAACAGCAGGCTCGGCGTTTCGAGGCACGTACGCAAATCCCCt  >  1:270282/1‑150 (MQ=255)                                                                                                                                                   tGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCg                                                   >  2:248692/1‑98 (MQ=255)                                                                                                                                                   tGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCg                                                   <  1:248692/98‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCTGGAGAGCGTGCGTGAAGAGGAAAAAGCCCGTATCGCCCGGGAAGTGCATGATGAACTGGGGCAGATGCTGACCGTGCTCAAGCTGGAAGTGTCGATGTGCGAGCTGGCCTATGCCGAGCTGGACCCTGGCCTGAATGATCGCCTGGCGAGCATGAAACGCCTGATCGCCCAGTTGTTCCAGTTGGTACGCGATGTCGCCACTGCCTTGCGTCCGCCGATCCTCGATGCCGGCATTGCCTCGGCCATCGAATGGCAGGCGCGGCGTTTCGAGGCACGTACGCAAATCCCCT  >  NC_002947/497904‑498196

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