Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,859,120	C→G	7.6%	intergenic (+40/+61)	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,120	0	C	G	7.6%	203.6 / 5.2	65	intergenic (+40/+61)	PP_1661/PP_1662	dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base G (3/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

GTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAG  >  NC_002947/1858974‑1859267                                                                                                                                                     |                                                                                                                                                    gTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCg‑cag                                                                                                                                                   >  7:257228/1‑149 (MQ=255)   gAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCg‑cagca                                                                                                                                                 >  3:233230/1‑149 (MQ=255)      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)                                         ttACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCg‑c                                                                                                                                                     <  2:309783/107‑1 (MQ=255)                                         ttACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCg‑c                                                                                                                                                     >  1:309783/1‑107 (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                                                                                                                       <  7:158367/104‑1 (MQ=255)                                                                          ccTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAgg                                                                                                                       >  8:158367/1‑104 (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                             >  3:379647/1‑144 (MQ=255)                                                                                                                            cccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTc                             <  4:379647/144‑1 (MQ=255)                                                                                                                            cccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGtt                        <  3:373513/149‑1 (MQ=255)                                                                                                                            cccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGtt                        >  4:382154/1‑149 (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)                                                                                                                                                     |                                                                                                                                                    GTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAACATGCCCACCAG  >  NC_002947/1858974‑1859267

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: ‑