Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,984,780	Δ1 bp	5.6%	intergenic (‑33/+20)	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,780	0	A	.	5.6%	211.8 / 6.8	72	intergenic (‑33/+20)	cheR/PP_4393	chemotaxis protein methyltransferase/chemotaxis protein
Reads supporting (aligned to +/- strand):  ref base A (31/37);  new base . (2/2);  total (33/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

ATCAGCTTGTTGAGACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAGTAT  >  NC_002947/4984632‑4984927                                                                                                                                                       |                                                                                                                                                     aTCAGCTTGTTGAGACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAg                                                                                                                                                     >  2:996825/1‑150 (MQ=255) aTCAGCTTGTTGAGACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAg                                                                                                                                                     <  1:948273/150‑1 (MQ=255)  tCAGCTTGTTGAGACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAgg                                                                                                                                                    <  2:1839398/150‑1 (MQ=255)  tCAGCTTGTTGAGACGGCTCGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAgg                                                                                                                                                    <  1:6585/150‑1 (MQ=255)      cTTGTTGAGACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTTCCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCAt                                                                                                                                                <  1:800065/150‑1 (MQ=255)      cTTGTTGAGACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGCCACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCAt                                                                                                                                                >  2:724863/1‑150 (MQ=255)               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                                                                                                                                     >  2:1500843/1‑132 (MQ=255)                                 cTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGG‑GGGCATGCCCTCCGgca                                                                                                                                     <  1:1500843/133‑2 (MQ=255)                                      tatccgcGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACGGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATc                                                                                                               <  2:1178314/144‑1 (MQ=255)                                          ctcGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCCATCCCAGTTTCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGt                                                                                                           <  1:1572986/150‑1 (MQ=255)                                            cGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTc                                                                                                          <  2:222027/149‑1 (MQ=255)                                            cGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTc                                                                                                          >  1:222027/1‑149 (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:625433/150‑1 (MQ=255)                                                                       ccAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGc                                                                              >  1:1857572/1‑150 (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                                                                                                                <  2:475124/94‑1 (MQ=255)                                                                                             gAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGAt                                                                                                                >  1:475124/1‑94 (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                                          <  2:1742414/150‑1 (MQ=255)                                                                                                           ccGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACcttc                                          <  1:213393/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:1702154/150‑1 (MQ=255)                                                                                                                               tctT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGaa                      <  1:1197351/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:846288/150‑1 (MQ=255)                                                                                                                                                 tGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGaa      <  1:508666/150‑1 (MQ=255)                                                                                                                                                  gCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAAGACCCCCGGCGTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAg     <  2:1894639/150‑1 (MQ=255)                                                                                                                                                     ggAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGGACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAGTAt  >  2:1503957/1‑150 (MQ=255)                                                                                                                                                       |                                                                                                                                                     ATCAGCTTGTTGAGACGGCTGGAAACCAGGTACTGCTTATTCTCGCCCAGCAGGATGCCACAGGCTTTCTCCAGGAAGACCCGGAACTGTTCGAAATCCAAATTACCCGTAGACACTACTGCCGCCTCTT‑‑TTCAATCACTGGTGCCGGAGGGCATGCCC‑CCGGCTGCTTCAATGCGCTGCCTTGATCCGGTCGACCACCCGCTGGGCCAGGTCATCCGGCTTGAACTTGGCCAGGAAGTCATCGGCACCGACCTTCTTGACCATTGCCTGGTTGAACACCCCGGAAAGCGAAGTAT  >  NC_002947/4984632‑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: ‑