Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,633,500	T→G	16.4%	intergenic (+100/‑82)	hupB → / → PP_2304	DNA binding regulator subunit beta B/peptidyl‑prolyl cis‑trans isomerase D

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,633,500	0	T	G	16.4%	156.1 / 26.9	73	intergenic (+100/‑82)	hupB/PP_2304	DNA binding regulator subunit beta B/peptidyl‑prolyl cis‑trans isomerase D
Reads supporting (aligned to +/- strand):  ref base T (33/28);  new base G (5/7);  total (38/35)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.17e-01

GTTCCAGGTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATCGTCGT  >  NC_002947/2633353‑2633649                                                                                                                                                     |                                                                                                                                                       gTTCCAGGTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCg                                                                                                                                                      <  2:29146/150‑1 (MQ=255)      aGGTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATgc                                                                                                                                                 >  2:1615196/1‑150 (MQ=255)        ttttCAAGGCCGGCAAAGGCCTCAAAGATGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCGAGGATgcgc                                                                                                                                               <  1:318622/149‑1 (MQ=255)        gTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCGAGGATgcgc                                                                                                                                               <  1:164547/150‑1 (MQ=255)              aGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑ctttct                                                                                                                                        >  2:2066470/1‑150 (MQ=255)              aGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCGAGGATGCGC‑ct                                                                                                                                            <  1:106143/146‑1 (MQ=255)              aGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCGAGGATGCGC‑ct                                                                                                                                            >  2:106143/1‑146 (MQ=255)                 ccGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCtttt                                                                                                                                     >  2:1624815/1‑150 (MQ=255)                        aGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGa                                                                                                                              <  1:623374/150‑1 (MQ=255)                             tGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTa                                                                                                                         >  2:332221/1‑150 (MQ=255)                             tGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTa                                                                                                                         >  2:774061/1‑150 (MQ=255)                                aaGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTAccc                                                                                                                      >  1:2001835/1‑150 (MQ=255)                                 aGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCa                                                                                                                     <  1:774061/150‑1 (MQ=255)                                     gCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATgc                                                                                                                                                 <  2:1536299/119‑1 (MQ=255)                                     gCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATgc                                                                                                                                                 >  1:1536299/1‑119 (MQ=255)                                     gCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCt                                                                                                                 >  1:659378/1‑150 (MQ=255)                                       cGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGGCACGTTTGCTAAGCAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTcc                                                                                                               >  2:1302362/1‑150 (MQ=255)                                          cAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAACACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCcgcg                                                                                                            <  1:1619637/150‑1 (MQ=255)                                             cTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGTATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGtt                                                                                                         >  1:738333/1‑150 (MQ=255)                                             cTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGtt                                                                                                         <  1:1893660/150‑1 (MQ=255)                                             cTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGtt                                                                                                         >  2:572015/1‑150 (MQ=255)                                                  tcgtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGa                                                                                                    <  1:332221/150‑1 (MQ=255)                                                    gtcTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACg                                                                                                  >  1:1370354/1‑150 (MQ=255)                                                         tCAGCCGGGCCCGGACTGGCCTGGTCTGAGCACGGTGCTGCCGGGGCGCAAACGTTCCCGCATGACACGTTTGCTAAGAACAGGAGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTg                                                                                             <  1:493029/150‑1 (MQ=255)                                                           aGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGt                                                                                           <  1:1132094/150‑1 (MQ=255)                                                               ggACCAGGCCAGGCCGGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAg                                                                                       >  1:566920/1‑150 (MQ=255)                                                                    cGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGa                                                                                                                              <  2:2163781/106‑1 (MQ=255)                                                                    cGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGa                                                                                                                              >  1:2163781/1‑106 (MQ=255)                                                                          gatctggTCCGAGCACTGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGATT‑ATAACCCTTTCCTCAAACCGCTGATGCGC‑CTTGCTAATATAAGGATTCTTCCCAGGATCAGCGGGGGTCGTCGCAGTGGTACAGCCGTTTCTggg                                                                            <  1:2388924/147‑1 (MQ=255)                                                                           gcctggTCCGAGCACGGTGATGGCTGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTgggg                                                                           >  2:1621214/1‑150 (MQ=255)                                                                              tggTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGAc                                                                        <  2:2181219/150‑1 (MQ=255)                                                                              tggTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGAc                                                                        <  2:311079/150‑1 (MQ=255)                                                                                gTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCATGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGc                                                                      >  1:2107614/1‑150 (MQ=255)                                                                                gTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTggggagata                                                                      >  2:2388924/1‑145 (MQ=255)                                                                                   cGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCAtg                                                                   <  1:2066470/150‑1 (MQ=255)                                                                                    gAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCAtgc                                                                  >  1:1085822/1‑150 (MQ=255)                                                                                       cACGGTGATGGCGGGGCGCAAACGTTCCCGCGTGACACGTTTGCTAAGAAAAGGCGC‑CTACTCGAAAG‑GCGAAGGCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCAtgctgc                                                               <  1:1615196/150‑1 (MQ=255)                                                                                         cGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGCCGCATGCTGCaa                                                             >  2:1310331/1‑150 (MQ=255)                                                                                              aTGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAatat                                                        >  2:1837940/1‑150 (MQ=255)                                                                                                  cGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAgg                                                    <  2:2001835/150‑1 (MQ=255)                                                                                                  cGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAgg                                                    <  2:1167405/150‑1 (MQ=255)                                                                                                  cGGGGCCCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAgg                                                    <  1:152499/150‑1 (MQ=255)                                                                                                         cAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCGAGGATGCGC‑cttttct                                                                                                                                       <  1:1078/60‑4 (MQ=21)                                                                                                         cAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCGAGGATGCGC‑cttttct                                                                                                                                       >  2:1078/1‑57 (MQ=21)                                                                                                           aaCGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑AGCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTca                                           <  2:738333/150‑1 (MQ=255)                                                                                                               ttCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAg                                       <  1:579999/150‑1 (MQ=255)                                                                                                                tCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAgg                                      <  2:659378/150‑1 (MQ=255)                                                                                                                   cGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTg                                   >  1:1907618/1‑150 (MQ=255)                                                                                                                      cTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGAt                                >  1:1550209/1‑150 (MQ=255)                                                                                                                          cacGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGcc                            <  2:566920/150‑1 (MQ=255)                                                                                                                             gTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAg                         >  1:1333177/1‑150 (MQ=255)                                                                                                                             gTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAg                         >  2:1926152/1‑150 (MQ=255)                                                                                                                                  cTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATACCAGGGACAATTCACAAGGTTGGa                                 >  2:684951/1‑137 (MQ=255)                                                                                                                                  cTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATACCAGGGACAATTCACAAGGTTGGa                                 <  1:684951/137‑1 (MQ=255)                                                                                                                                   tAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACcatc                   <  1:1811646/150‑1 (MQ=255)                                                                                                                                   tAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACcatc                   >  2:2314533/1‑150 (MQ=255)                                                                                                                                    aaGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACcatca                  >  2:954165/1‑150 (MQ=255)                                                                                                                                     aagaaaGGCGC‑ATCCGAGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACcatcat                 >  2:2305963/4‑150 (MQ=255)                                                                                                                                       gaaaGGCGC‑ATCCGAGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGtt                                    >  1:581284/2‑129 (MQ=255)                                                                                                                                       gaaaGGCGC‑ATCCGAGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGtt                                    <  2:581284/128‑1 (MQ=255)                                                                                                                                       aaaaGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCg               >  2:2245531/1‑150 (MQ=255)                                                                                                                                       aaaaGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCg               >  1:1160420/1‑150 (MQ=255)                                                                                                                                        aaaGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCgg              >  2:1274352/1‑150 (MQ=255)                                                                                                                                         aaGTCGC‑ATCCTCGGATGCGC‑CTTTCTTCTATCAGGATTCTACCAACGCTCCGCGGTTGTCGACGCAGTGGTACAGCAGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGc             <  2:1085822/150‑1 (MQ=255)                                                                                                                                         aaGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCATGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGc             <  2:2107614/150‑1 (MQ=255)                                                                                                                                         aaGGCGC‑ATCCGAGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGc             >  2:2393912/1‑150 (MQ=255)                                                                                                                                         aaGGCGC‑ATCCGAGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGc             <  1:2305963/150‑1 (MQ=255)                                                                                                                                          aGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGcc            <  2:1333177/150‑1 (MQ=255)                                                                                                                                               c‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATACCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTAtc       <  1:1624815/150‑1 (MQ=255)                                                                                                                                                 aTCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGAtt                               >  2:1634500/1‑125 (MQ=255)                                                                                                                                                 aTCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGAtt                               <  1:1634500/125‑1 (MQ=255)                                                                                                                                                 aTCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTAtcg      >  2:169179/1‑150 (MQ=255)                                                                                                                                                    cTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGTCGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTAtcgtcg   >  2:998781/1‑150 (MQ=255)                                                                                                                                                     tCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATcgtcgt  >  2:66512/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                       GTTCCAGGTTTCAAGGCCGGCAAAGGCCTGAAAGACGCCGTCAACTAAGTCGTCTTTCAGCCGGACCCGGCCTGGCCTGGTCCGAGCACGGTGATGGCGGGGCGCAAACGTTCCCGCCTGACACGTTTGCTAAGAAAAGGCGC‑ATCCTCGGATGCGC‑CTTTCTTTTATCAGGATTCTACCCACGCTCCGCGGTTGTCGACGCAGTGGTACAGCCGTTTCTGGGGGACGCATGCTGCAAAATATCAGGGACAATTCACAAGGTTGGATTGCCAAGACCATCATCGGCCTTATCGTCGT  >  NC_002947/2633353‑2633649

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