Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,859,108	C→G	7.3%	intergenic (+28/+73)	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,108	0	C	G	7.3%	216.7 / 5.3	68	intergenic (+28/+73)	PP_1661/PP_1662	dehydrogenase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (35/28);  new base G (3/2);  total (38/30)
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.79e-01

GCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAA  >  NC_002947/1858960‑1859255                                                                                                                                                     |                                                                                                                                                      gcgTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTc                                                                                                                                                       <  2:95137/149‑1 (MQ=255) gcgTGGCCAGCAAGGGGAGTTTCAGCTGTGCTTCAATGACGACCTGGACGGGCGTTACGGGGCGGGCCTGATCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCTGGGCTGTGTAGGCCCTGTCATGGCAACCCCACGATGCTc                                                                                                                                                       <  8:355029/149‑1 (MQ=255)   gTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCg                                                                                                                                                     >  5:263892/1‑149 (MQ=255)       ccAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCa                                                                                                                                                  >  6:61863/1‑148 (MQ=255)        cAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAAt                                                                                                                                                >  6:279680/1‑149 (MQ=255)         aGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATg                                                                                                                                               >  4:307494/1‑149 (MQ=255)               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:114389/1‑149 (MQ=255)                        aGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTggg                                                                                                                              <  5:24068/149‑1 (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                                                                                                                                         >  1:309783/1‑107 (MQ=255)                                                       ttACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCg‑c                                                                                                                                         <  2:309783/107‑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                                                                                                           <  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                                                                                 <  8:314244/117‑1 (MQ=255)                                                                                                     gTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGc                                                                                 >  7:314244/1‑117 (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            <  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           <  6:114389/149‑1 (MQ=255)                                                                                                                                           ccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc           <  6:112904/149‑1 (MQ=255)                                                                                                                                           ccACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTc           <  8:329802/149‑1 (MQ=255)                                                                                                                                            cACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCg          >  4:344172/1‑149 (MQ=255)                                                                                                                                            cACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCg          >  5:176776/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)                                                                                                                                                     |                                                                                                                                                      GCGTGGCCAGCAAGGTGAGTTGCAGCTGTGCATCAATGACGACCTGGACGGGCGTTACGGGGCTGGCCTGAGCGATAACCTGGGCAGCCTGTCGGTAGAGGTGCGTTTCGGGGCGGTGTAGGCCCTGTCATGGCAACCCCACGATGCTCCGGCAATG‑CCG‑CAGCATTTTGGGGTTGCTCACTGAGCCCAGGCGAGTCAGGGGTTCGCCGCAATGGGCAGGGTCATTGTGTGGCGCGCTGGTTTTCCAGGTCTTCAAGCTCAACGCGGTATTGGTGTGCTTCGTGTTCGTTATGGAA  >  NC_002947/1858960‑1859255

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