Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,859,127	T→G	8.8%	intergenic (+47/+54)	PP_1661 → / ← PP_1662	dehydrogenase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,859,127	0	T	G	8.8%	196.5 / 5.9	66	intergenic (+47/+54)	PP_1661/PP_1662	dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base G (3/3);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAGTT  >  NC_002947/1858979‑1859269                                                                                                                                                       |                                                                                                                                               ttGCAGCTGTGCATCAATGACGACCTGGACGGGCTTTAGGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCAttt                                                                                                                                                <  1:146169/149‑1 (MQ=255) ttGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCAttt                                                                                                                                                >  5:259416/1‑149 (MQ=255)     aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTggg                                                                                                                                            >  5:112904/1‑149 (MQ=255)     aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTggg                                                                                                                                            <  5:24068/149‑1 (MQ=255)     aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTggg                                                                                                                                            >  5:114389/1‑149 (MQ=255)       cTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGt                                                                                                                                          <  5:217783/149‑1 (MQ=255)             aTCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCa                                                                                                                                    <  6:263892/149‑1 (MQ=255)                aaTGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACt                                                                                                                                  >  3:307673/1‑148 (MQ=255)                    acgacCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGc                                                                                                                              <  3:307494/148‑1 (MQ=255)                     cgacCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGcc                                                                                                                             >  4:126249/1‑148 (MQ=255)                               gggCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGt                                                                                                                    >  8:144057/1‑147 (MQ=255)                                ggCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAg                                                                                                                 <  7:144057/149‑1 (MQ=255)                                 gCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGGGTCAgg                                                                                                                <  6:121709/149‑1 (MQ=255)                                 gCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCACGCGAGTCAgg                                                                                                                <  6:26474/149‑1 (MQ=255)                                        ggggCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTcgc                                                                                                         <  3:126249/149‑1 (MQ=255)                                              ggCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAAt                                                                                                   <  6:259416/149‑1 (MQ=255)                                                   gAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGc                                                                                               >  5:318228/1‑148 (MQ=255)                                                          aaCCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAt                                                                                       <  5:13893/149‑1 (MQ=255)                                                           aCCTGGGCAGCCTGTCGGTAGATGTGCGTTTCGGGGGGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCATTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAtt                                                                                      <  1:81231/149‑1 (MQ=255)                                                           aCCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CAG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCAt                                                                                       <  7:163128/148‑1 (MQ=255)                                                                ggCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgtg                                                                                 <  2:61414/149‑1 (MQ=255)                                                                ggCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgt                                                                                  <  8:257228/148‑1 (MQ=255)                                                                ggCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATtgtgt                                                                                  <  4:307673/148‑1 (MQ=255)                                                                 gCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTgg                                                                                >  4:373513/1‑149 (MQ=255)                                                                     ccTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAgg                                                                                                                         >  8:158367/1‑104 (MQ=255)                                                                     ccTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAgg                                                                                                                         <  7:158367/104‑1 (MQ=255)                                                                     ccTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGgcgcg                                                                            >  2:319443/1‑149 (MQ=255)                                                                       tGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGTGCAGGGTCATTGTGTGGCGCGCt                                                                          >  5:393376/1‑149 (MQ=255)                                                                         tCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTg                                                                         <  5:61863/148‑1 (MQ=255)                                                                           ggTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAgg                                                                                                                <  5:99164/107‑1 (MQ=255)                                                                           ggTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAgg                                                                                                                >  6:99164/1‑107 (MQ=255)                                                                            gTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGTGTCAGGGGTTCGCGGCAATGGGCAGGGTGATTGTGTGGCGCGCTGGttt                                                                     >  5:393021/1‑149 (MQ=255)                                                                            gTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGttt                                                                     >  6:152520/1‑149 (MQ=255)                                                                                 ggTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATGCCCGACAGAATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGAAGGGTCATTGGTTGGCGCGCTGGTTTTc                                                                   >  1:111608/1‑148 (MQ=255)                                                                                  gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGc                                                                                               >  7:314244/1‑117 (MQ=255)                                                                                  gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGc                                                                                               <  8:314244/117‑1 (MQ=255)                                                                                           ggggCTGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGTGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAg                                                       <  5:110611/148‑1 (MQ=255)                                                                                              gCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTc                                                    >  1:267259/1‑148 (MQ=255)                                                                                                   gtAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGCGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAAcgcg                                              <  4:233230/149‑1 (MQ=255)                                                                                                       gCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGTGCAGGGTCATTGTGTGGCGCGCTGGTTTTACAGGTCTTCAAGCTCAACGCGGTAt                                          >  1:219649/1‑149 (MQ=255)                                                                                                         ccTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTg                                        >  3:312966/1‑149 (MQ=255)                                                                                                               aTGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGc                                  >  7:214959/1‑149 (MQ=255)                                                                                                                gagCAACCCCAAAATGCTGCGGCATTG‑CCG‑GAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCt                                 >  2:306925/3‑149 (MQ=37)                                                                                                                  gCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc                               >  7:329802/1‑149 (MQ=255)                                                                                                                     aCCCCAAAATGCTGCGGCATTG‑CCG‑GAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCgtg                            >  3:30367/1‑149 (MQ=38)                                                                                                                      ccccAAAATGCTGCGGCATTG‑CCG‑GAGCATCGTGGGGTTGC‑CACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGtt                          >  3:174396/1‑149 (MQ=38)                                                                                                                       cccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc                               <  4:379647/144‑1 (MQ=255)                                                                                                                       cccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc                               >  3:379647/1‑144 (MQ=255)                                                                                                                       cccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGtt                          >  4:382154/1‑149 (MQ=255)                                                                                                                       cccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGtt                          <  3:373513/149‑1 (MQ=255)                                                                                                                       cccACGATGCTCAGGCAATG‑CAG‑CAGCATTTTGGGGTTGCTCACTGAGACCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGtt                          >  2:220603/1‑149 (MQ=255)                                                                                                                        ccaaaaTGCTGCGGCATTG‑CCG‑GAGCATCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                         <  1:306925/144‑1 (MQ=38)                                                                                                                        ccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                         <  8:329802/149‑1 (MQ=255)                                                                                                                        ccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                         <  6:114389/149‑1 (MQ=255)                                                                                                                        ccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc                         <  6:112904/149‑1 (MQ=255)                                                                                                                         cACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCg                        >  5:176776/1‑149 (MQ=255)                                                                                                                         cACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCg                        >  4:344172/1‑149 (MQ=255)                                                                                                                            aaTGCTGCGGCATTG‑CCG‑GAGCATCGTGGGGTTGC‑CACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTAt                    <  4:174396/148‑1 (MQ=38)                                                                                                                               gCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATg                   >  2:308486/1‑148 (MQ=255)                                                                                                                                 tCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGaa                >  1:188469/1‑149 (MQ=255)                                                                                                                                    ggCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACAt             <  6:176776/149‑1 (MQ=255)                                                                                                                                        aTG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGccc         <  3:344172/149‑1 (MQ=255)                                                                                                                                          g‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCccac       >  5:324073/1‑149 (MQ=255)                                                                                                                                          g‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCccac       >  2:332062/1‑149 (MQ=255)                                                                                                                                            ccG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGTTTCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCccacc      <  1:319443/149‑1 (MQ=255)                                                                                                                                              g‑cagcaTTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAg    >  6:6589/1‑149 (MQ=255)                                                                                                                                                 agcaTTTTGGGGTTGCTCACTGAGGCGAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAGtt  <  2:111608/149‑1 (MQ=255)                                                                                                                                                 agcaTCGTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAGtt  <  4:30367/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                               TTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAGTT  >  NC_002947/1858979‑1859269

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