Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,984,790:1	+T	8.8%	intergenic (‑43/+10)	cheR ← / ← PP_4393	chemotaxis protein methyltransferase/chemotaxis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,984,790	1	.	T	8.8%	248.4 / 18.6	68	intergenic (‑43/+10)	cheR/PP_4393	chemotaxis protein methyltransferase/chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base . (29/33);  new base T (4/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

ACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAGTAT  >  NC_002947/4984646‑4984927                                                                                                                                                    |                                                                                                                                          actgCTGGAAACCAGGTACCGCTCATCCTCGTTTAGTAGGATGCCACAGGCTTGCAACAGAACGACCAGGACCTGTTCGAGATCCAAATGACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGG‑GGGCATGCCCTCCGgc                                                                                                                                      <  2:1348764/147‑1 (MQ=255) aCGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGgc                                                                                                                                      >  1:1410241/1‑150 (MQ=255) aCGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGG‑GGGCATGCCCTCCGgc                                                                                                                                      >  1:801590/1‑150 (MQ=255)      tGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCtt                                                                                                                                 <  2:1869313/150‑1 (MQ=255)       ggAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTc                                                                                                                                <  2:1497500/150‑1 (MQ=255)        gAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCa                                                                                                                               >  2:568569/1‑150 (MQ=255)             cAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATgcg                                                                                                                          >  2:1481328/1‑150 (MQ=255)              aGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATgcgc                                                                                                                         <  2:779091/150‑1 (MQ=255)                   cTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCt                                                                                                                    <  1:997031/150‑1 (MQ=255)                   cTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGG‑GGGCATGCCCTCCGgca                                                                                                                                     <  1:1500843/133‑2 (MQ=255)                   cTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGG‑GGGCATGCCCTCCGgca                                                                                                                                     >  2:1500843/1‑132 (MQ=255)                        tatccgcGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACGGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATc                                                                                                               <  2:1178314/144‑1 (MQ=255)                            ctcGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCCATCCCAGTTTCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGt                                                                                                           <  1:1572986/150‑1 (MQ=255)                              cGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTc                                                                                                          >  1:222027/1‑149 (MQ=255)                              cGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTc                                                                                                          <  2:222027/149‑1 (MQ=255)                                         aTGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATG‑CC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTgg                                                                                             >  1:1198993/1‑150 (MQ=255)                                          tGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTgg                                                                                             <  1:1481328/150‑1 (MQ=255)                                              acaGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCa                                                                                         >  2:715936/1‑150 (MQ=255)                                               caGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAg                                                                                        <  2:1105688/150‑1 (MQ=255)                                               caGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAg                                                                                        <  1:1887937/150‑1 (MQ=255)                                                         ccAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGc                                                                              >  1:1857572/1‑150 (MQ=255)                                                         ccAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGc                                                                              <  1:625433/150‑1 (MQ=255)                                                          cAGGAAGACCCGGAACTGTTCGAAATGCAACTTACCGGTAGGCACT‑CT‑ACGCTTCGTACCGCAATCACTGGCTCCTGAGGGCATGACC‑CAGGCTGCTTCAATGGGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCGTCCGGCt                                                                             >  1:514020/1‑150 (MQ=255)                                                          cAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCt                                                                             <  1:996825/150‑1 (MQ=255)                                                                 aCCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTgg                                                                                             <  1:607785/127‑1 (MQ=255)                                                                 aCCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTgg                                                                                             >  2:607785/1‑127 (MQ=255)                                                                   ccGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTgg                                                                    >  1:1312469/1‑150 (MQ=255)                                                                   ccGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTgg                                                                    >  1:171979/1‑150 (MQ=255)                                                                       aaCTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAg                                                                <  1:715936/150‑1 (MQ=255)                                                                        aCTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAgg                                                               >  2:508666/1‑150 (MQ=255)                                                                          tGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGaa                                                             <  1:1865251/150‑1 (MQ=255)                                                                            ttCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGt                                                           >  2:454278/1‑150 (MQ=255)                                                                               gAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGAt                                                                                                                >  1:475124/1‑94 (MQ=255)                                                                               gAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGAt                                                                                                                <  2:475124/94‑1 (MQ=255)                                                                                aaaTCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAgtcgtc                                                       >  1:1508160/1‑150 (MQ=255)                                                                                    ccAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGCAGTCATCGGcc                                                   >  1:717042/1‑149 (MQ=255)                                                                                         ttACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAAGGTGGCCAGGAAGTCATCGGCACCGAc                                              >  1:292588/1‑150 (MQ=255)                                                                                             ccGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACcttc                                          <  1:213393/150‑1 (MQ=255)                                                                                             ccGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACcttc                                          <  2:1742414/150‑1 (MQ=255)                                                                                             ccGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCGGTGCTGGGTCATCCGGCTGGAGCTTGGCCAGGAAGTCATCGGCACCGACcttc                                          >  2:990750/1‑150 (MQ=255)                                                                                              cGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACcttct                                         >  2:855209/1‑150 (MQ=255)                                                                                                  gACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGc                                                                   <  2:982395/120‑1 (MQ=255)                                                                                                  gACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGc                                                                   >  1:982395/1‑120 (MQ=255)                                                                                                  gACACTACTGCCGCCTCTT‑‑TTCAAGCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTTATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAGCTTGCCAAGGAAGTTATCGGTACCGACCTTCTTGAc                                     >  1:822255/1‑150 (MQ=255)                                                                                                       tactGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACctt                                           <  2:1009771/139‑1 (MQ=255)                                                                                                       tactGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACctt                                           >  1:1009771/1‑139 (MQ=255)                                                                                                            ccgncTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTgg                           >  2:1197351/1‑150 (MQ=255)                                                                                                             cgccTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGt                          <  2:1312469/150‑1 (MQ=255)                                                                                                              gccTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGtt                         <  1:38310/150‑1 (MQ=255)                                                                                                              gccTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGtt                         <  2:292588/150‑1 (MQ=255)                                                                                                              gccTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATCCTGTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGtt                         >  1:1078454/1‑150 (MQ=255)                                                                                                               ccTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCAtt                                 <  2:388238/141‑1 (MQ=255)                                                                                                               ccTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCAtt                                 >  1:388238/1‑141 (MQ=255)                                                                                                               ccTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTg                        >  2:1702154/1‑150 (MQ=255)                                                                                                                 tctT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGaa                      <  1:1197351/150‑1 (MQ=255)                                                                                                                 tctT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGaa                      <  1:1702154/150‑1 (MQ=255)                                                                                                                   tt‑‑ttCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAg                                                                                        >  2:243652/1‑82 (MQ=255)                                                                                                                   tt‑‑ttCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAg                                                                                        <  1:243652/82‑1 (MQ=255)                                                                                                                              aCTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGAt                                                                                                                >  1:651422/1‑49 (MQ=255)                                                                                                                              aCTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGAt                                                                                                                <  2:651422/49‑1 (MQ=255)                                                                                                                                tgttGCCGGAGGGCATGCCC‑CCGGCTGCTTGAAGCCCCTTCCTTGATTCTGGCGTCTATCCGCTGGGCCAGCTCATCCGTCTTTAGCTTGGCCAGGAATTCATCCTGACCTGCCTTCTTATCCATTGCCTTTTAGTACACCGCGGAAAGc         <  2:822255/147‑1 (MQ=255)                                                                                                                                  gTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGa       <  1:724863/150‑1 (MQ=255)                                                                                                                                   tGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGaa      <  1:508666/150‑1 (MQ=255)                                                                                                                                   tGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGaa      <  1:846288/150‑1 (MQ=255)                                                                                                                                   cagccgggggCATGCCCTCCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGa       >  2:767045/8‑150 (MQ=255)                                                                                                                                    gCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAAGACCCCCGGCGTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAg     <  2:1894639/150‑1 (MQ=255)                                                                                                                                      ccgggggCATGCCCTCCGGCTGCTTCGATACGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAGt    >  1:1792213/5‑150 (MQ=255)                                                                                                                                       ggAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGGACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAGTAt  >  2:1503957/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                          ACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAGTAT  >  NC_002947/4984646‑4984927

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