Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	318,938	C→T	7.4%	S269N (AGC→AAC)	PP_0264 ←	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	318,938	0	C	T	7.4%	188.6 / 4.9	67	S269N (AGC→AAC)	PP_0264	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base T (3/2);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CACGGGCTTCGAAGTCGATGCGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGT  >  NC_002947/318789‑319080                                                                                                                                                      |                                                                                                                                               cACGGGCTTCGAAGTCGATGCGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGc                                                                                                                                                >  2:543634/1‑150 (MQ=255)  aCGGGCTTCGAAGTCGATGCGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCt                                                                                                                                               <  2:953468/150‑1 (MQ=255)  aCGGGCTTCGAAGTCGATGCGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCt                                                                                                                                               <  1:747205/150‑1 (MQ=255)     ggCTTCGAAGTCGATGCGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGAc                                                                                                                                            <  1:58783/150‑1 (MQ=255)      gCTTCGAAGTCGATGCGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACt                                                                                                                                           <  2:414401/150‑1 (MQ=255)                 gATGCGGTCGTTGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGt                                                                                                                                >  2:593219/1‑150 (MQ=255)                 gATGCGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGt                                                                                                                                >  2:833602/1‑150 (MQ=255)                     cGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGTCCTGCTTGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTcc                                                                                                                            <  1:1096524/150‑1 (MQ=255)                         cGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGa                                                                                                                        >  2:221442/1‑150 (MQ=255)                               aaTGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTcgc                                                                                                                  >  2:192420/1‑150 (MQ=255)                                aTGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTcgcg                                                                                                                 >  2:181804/1‑150 (MQ=255)                                 tGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCGACATCAGGGCCAGCAGCGTACCTGCACCGATGGCCATTACCGTCTGCCCCGAGCGATCAACAAGCGCGCCCGGGGCGCCGACGCTGACGTGCCACCCAGTTTCCTTGATGTCGCgg                                                                                                                >  2:346208/1‑150 (MQ=255)                                 tGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCgg                                                                                                                <  2:72398/150‑1 (MQ=255)                                        gccgTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTggg                                                                                                         <  1:181804/150‑1 (MQ=255)                                                aTCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGTTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAg                                                                                                 >  2:194460/1‑150 (MQ=255)                                                aTCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAg                                                                                                 <  1:833602/150‑1 (MQ=255)                                                   aCCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGcc                                                                                              >  1:186813/1‑150 (MQ=255)                                                        gCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTATCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCgg                                                                                         <  2:583057/150‑1 (MQ=255)                                                        gCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCgg                                                                                         >  2:353088/1‑150 (MQ=255)                                                          cTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGt                                                                                       >  2:549550/1‑150 (MQ=255)                                                             ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGa                                                                                                      <  1:935443/132‑1 (MQ=255)                                                             ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGa                                                                                                      >  2:935443/1‑132 (MQ=255)                                                             ccAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTga                                                                                    <  1:593219/150‑1 (MQ=255)                                                              cAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTgag                                                                                   >  1:440152/1‑150 (MQ=255)                                                                  aaCATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAAcc                                                                               >  2:1113654/1‑150 (MQ=255)                                                                   aCATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGTTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCa                                                                              >  1:328332/1‑150 (MQ=255)                                                                     aTCAGTACCAGCAGCGTACCGGCACCGATGGCCAGTACCGTCTGCACCGAGCGATCAACCAGCGTTAGGGTATTCAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGc                                                                            <  1:192420/150‑1 (MQ=255)                                                                      tCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCg                                                                           >  1:112378/1‑150 (MQ=255)                                                                         gTACCAGCAGCGTACCCGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCt                                                                        <  2:567916/150‑1 (MQ=255)                                                                             cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTgg                                                                    <  1:549550/150‑1 (MQ=255)                                                                             cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTgg                                                                    >  1:966030/1‑150 (MQ=255)                                                                             cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTgg                                                                    >  1:957928/1‑150 (MQ=255)                                                                             cagcagCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTgg                                                                    >  1:822457/1‑150 (MQ=255)                                                                                  gcGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCt                                                               <  2:822457/150‑1 (MQ=255)                                                                                    gTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTgg                                                             <  2:80464/150‑1 (MQ=255)                                                                                                   gCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTgg                                              >  1:276035/1‑150 (MQ=255)                                                                                                    ccATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGa                                             >  1:594743/1‑150 (MQ=255)                                                                                                     cATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGTTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGAt                                            >  2:134362/1‑150 (MQ=255)                                                                                                      aTTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCtt                                                                                                                          <  1:299619/71‑1 (MQ=255)                                                                                                      aTTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCtt                                                                                                                          >  2:299619/1‑71 (MQ=255)                                                                                                       ttACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATgg                                          <  1:496944/150‑1 (MQ=255)                                                                                                          ccGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGTTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCga                                       <  1:134362/150‑1 (MQ=255)                                                                                                              cTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatgat                                   <  2:440152/150‑1 (MQ=255)                                                                                                                 cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCt                                                               >  2:875632/1‑119 (MQ=255)                                                                                                                 cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCt                                                               <  1:875632/119‑1 (MQ=255)                                                                                                                 cACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGc                                <  1:1113654/150‑1 (MQ=255)                                                                                                                   ccGAGCGATCAACCAGTGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCt                              <  2:961138/150‑1 (MQ=255)                                                                                                                    cGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTg                             <  2:186813/150‑1 (MQ=255)                                                                                                                    cGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTg                             <  2:135807/150‑1 (MQ=255)                                                                                                                     gAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTgg                            >  1:24297/1‑150 (MQ=255)                                                                                                                        cGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGgcgc                         <  2:878395/150‑1 (MQ=255)                                                                                                                         gATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCt                        >  2:809334/1‑150 (MQ=255)                                                                                                                          aTCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCtt                                                                                                                          <  1:638796/51‑1 (MQ=255)                                                                                                                          aTCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCtt                                                                                                                          >  2:638796/1‑51 (MQ=255)                                                                                                                              aCCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTgg                                                                    >  2:517113/1‑101 (MQ=255)                                                                                                                              aCCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTgg                                                                    <  1:517113/101‑1 (MQ=255)                                                                                                                               ccAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatga                                    <  1:1038987/132‑1 (MQ=255)                                                                                                                               ccAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCgatga                                    >  2:1038987/1‑132 (MQ=255)                                                                                                                               ccAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGt                  >  1:698370/1‑150 (MQ=255)                                                                                                                                cAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCCCGGCGTc                 >  2:58997/1‑150 (MQ=255)                                                                                                                                 aGCGTGCGCGGGGCGAGGATGTTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCg                <  2:328332/150‑1 (MQ=255)                                                                                                                                 aGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCg                >  2:1039054/1‑150 (MQ=255)                                                                                                                                     ggcgcGGGGCGAGGATGCTGACTTGCCAATCTGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCa            <  1:504810/149‑1 (MQ=255)                                                                                                                                           gggCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGc      <  1:924659/150‑1 (MQ=255)                                                                                                                                              cGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCggg   <  1:221442/150‑1 (MQ=255)                                                                                                                                               gAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTGTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGACATGTTGGCATGGCGCTCGGCGTCGGTCAGCGTGCGGGt  >  1:608109/1‑150 (MQ=255)                                                                                                                                                  gATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGTgg    >  1:401265/1‑145 (MQ=255)                                                                                                                                                  gATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGTgg    <  2:401265/145‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                               CACGGGCTTCGAAGTCGATGCGGTCGATGTAATGGCGCCGCCGTTGCATCACCAGGCCTGCCAGCAACATCAGTACCAGCAGCGTACCGGCACCGATGGCCATTACCGTCTGCACCGAGCGATCAACCAGCGTGCGCGGGGCGAGGATGCTGACTTGCCAACCAGTTTCCTTGATGTCGCGGGTCTGGGTGATCCAGGCCTGCGGGTTGAGAACCAGCGGCTGTGGGGCCTGGGTCGGATAGGGCTGGATGGCGATGATGGCCTGGCGCTCGGCGTCGGTCAGCGTGCGGGT  >  NC_002947/318789‑319080

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