Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,408,883	A→G	5.6%	Y12H (TAC→CAC)	PP_3892 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,408,883	0	A	G	5.6%	208.1 / 2.1	71	Y12H (TAC→CAC)	PP_3892	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (31/36);  new base G (2/2);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

TTCTCCATGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGG  >  NC_002947/4408737‑4409028                                                                                                                                                       |                                                                                                                                                  ttctCCATGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGc                                                                                                                                                 >  8:123002/1‑149 (MQ=255) ttctCCATGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGc                                                                                                                                                 >  2:135736/1‑149 (MQ=255)   ctcCATGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCgg                                                                                                                                               >  2:369624/1‑149 (MQ=255)   ctcCATGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCg                                                                                                                                                <  6:241469/148‑1 (MQ=255)      cATGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTAt                                                                                                                                            >  1:211011/1‑149 (MQ=255)      cATGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTAt                                                                                                                                            >  7:47416/1‑149 (MQ=255)       aTGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTAtt                                                                                                                                           <  1:139482/149‑1 (MQ=255)       aTGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTAtt                                                                                                                                           >  2:212964/1‑149 (MQ=255)         gAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGCAGAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTc                                                                                                                                          >  1:337131/1‑149 (MQ=255)                aaTAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACtt                                                                                                                                  >  8:4748/1‑149 (MQ=255)                  tAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCt                                                                                                                                >  7:297947/1‑149 (MQ=255)                    aCGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTgg                                                                                                                              <  7:59619/149‑1 (MQ=255)                       cAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTggcgg                                                                                                                           <  5:80161/149‑1 (MQ=255)                           aGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGgcagc                                                                                                                       <  2:211011/149‑1 (MQ=255)                              tGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGgcagcagc                                                                                                                    <  7:123002/149‑1 (MQ=255)                               gCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGgcagcagca                                                                                                                   >  7:71069/1‑149 (MQ=255)                               gCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGgcagcagca                                                                                                                   >  6:403429/1‑149 (MQ=255)                                cTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCAt                                                                                                                  <  3:383678/149‑1 (MQ=255)                                 tCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGAGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCAt                                                                                                                  >  4:53622/1‑148 (MQ=255)                                     ccACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCTGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGaa                                                                                                             <  8:98001/149‑1 (MQ=255)                                      cacaGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGTCTTCTGGCGGCAGCAGCATCAGAAc                                                                                                            <  5:403429/149‑1 (MQ=255)                                      cacaGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGaa                                                                                                             >  7:332963/1‑148 (MQ=255)                                           gcgcTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCT‑GGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTcc                                                                                                      >  2:53777/1‑149 (MQ=255)                                               tGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTGGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCct                                                                                                   >  4:248859/1‑149 (MQ=255)                                               tGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTGGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCct                                                                                                   >  4:254084/1‑149 (MQ=255)                                                 gCCGTGTGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAGCCCCTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTccc                                                                                                     <  1:135736/148‑1 (MQ=255)                                                   cGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCctggct                                                                                               >  2:414810/1‑149 (MQ=255)                                                     tGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCctggctgg                                                                                             <  7:283667/149‑1 (MQ=255)                                                       aGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGCCTCTGCGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCCGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCctgcctgg                                                                                             >  3:343679/1‑147 (MQ=255)                                                                   aaTCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCa                                                                                <  8:385250/149‑1 (MQ=255)                                                                   aaTCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCa                                                                                >  3:37847/1‑149 (MQ=255)                                                                     tCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAgttgtt                                                                                     <  3:319441/142‑1 (MQ=255)                                                                     tCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAgttgtt                                                                                     >  4:319441/1‑142 (MQ=255)                                                                              cGGTAGAGCGCAATTGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTGGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGggag                                                                     <  3:248859/149‑1 (MQ=255)                                                                              cGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAgttg                                                                                       >  7:234629/1‑131 (MQ=255)                                                                              cGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAgttg                                                                                       <  8:240915/131‑1 (MQ=255)                                                                              cGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAgttg                                                                                       <  8:234629/131‑1 (MQ=255)                                                                              cGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAgttg                                                                                       >  7:240915/1‑131 (MQ=255)                                                                              cGGTAGAGCGCAATGGGCGGCTCTGGGTTGGGGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTGGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGggag                                                                     <  3:254084/149‑1 (MQ=255)                                                                                gTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGCCAGGTCCAGCAG‑‑‑GTCGATGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGggagtc                                                                   >  3:230950/1‑147 (MQ=255)                                                                                    agCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTg                                                               <  5:42900/149‑1 (MQ=255)                                                                                     gcgcAATGGGCGGCTGTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTg                                                               <  4:230950/148‑1 (MQ=255)                                                                                     gcgcAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGc                                                              <  8:288256/149‑1 (MQ=255)                                                                                     gcgcAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGc                                                              <  8:71069/149‑1 (MQ=255)                                                                                       tcAATGGGCGGCTGTGGGTTGGTGGACAGGTCCTGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGTGGACCTGCGa                                                            <  4:343679/148‑1 (MQ=255)                                                                                             ggcggcTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGGAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGt                                                      <  8:297947/149‑1 (MQ=255)                                                                                               cggcTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGtt                                                     >  2:124280/1‑148 (MQ=255)                                                                                               cggcTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTg                                                    <  4:188799/149‑1 (MQ=255)                                                                                                ggcTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGtt                                                     >  6:166393/1‑147 (MQ=255)                                                                                                 gcTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAg                                                  <  1:369624/149‑1 (MQ=255)                                                                                                        gTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGctgctgc                                           >  4:21568/1‑149 (MQ=255)                                                                                                            gtggACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGctgctgcctgc                                       >  6:414438/1‑149 (MQ=255)                                                                                                                aCAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCctgctgct                                   >  7:175204/1‑149 (MQ=255)                                                                                                                  aGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCtgctgctgc                                 >  2:364929/1‑149 (MQ=255)                                                                                                                    gTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGAGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTc                               >  6:188256/1‑149 (MQ=255)                                                                                                                    gTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCGGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTc                               <  7:4748/149‑1 (MQ=255)                                                                                                                       cagcag‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGAGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTcgc                             <  5:188256/148‑1 (MQ=255)                                                                                                                        agcag‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACGTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCAGCAGCAGGCAGCAGCTcaacc                           <  3:106492/149‑5 (MQ=255)                                                                                                                                   gCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGc                   <  4:37847/149‑1 (MQ=255)                                                                                                                                     tGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTcc                      >  7:181991/1‑144 (MQ=255)                                                                                                                                     tGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTcc                      <  8:181991/144‑1 (MQ=255)                                                                                                                                       cGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTa                <  8:47416/148‑1 (MQ=255)                                                                                                                                       cGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAg               >  5:224247/1‑149 (MQ=255)                                                                                                                                         gCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGtt             <  5:166393/149‑1 (MQ=255)                                                                                                                                          cAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTTCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTg            <  2:129372/149‑1 (MQ=255)                                                                                                                                          cAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTg            <  6:224247/149‑1 (MQ=255)                                                                                                                                           aGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTgg           <  1:212964/149‑1 (MQ=255)                                                                                                                                           aGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTgg           >  5:81006/1‑149 (MQ=255)                                                                                                                                            gTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGc          <  3:21568/149‑1 (MQ=255)                                                                                                                                                 gcgcGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAAcc     <  6:196288/149‑1 (MQ=255)                                                                                                                                                    cgTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCggg  <  8:332963/149‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  TTCTCCATGAAGGTCAATAACGCAGAAGGTGCTCGGCCACAGGCGCTGGCCGTGAGCGATGGC‑GAAATCTTTGTTGCGGTAGAGCGCAATGGGCGGCTCTGGGTTGGTGGACAGGTCCAGCAG‑‑‑GTCGCTGCGGCAGTACAGCGCGTAGCGGTATTCGGTGACTTCTGGCGGCAGCAGCATCAGAACCTCTCCCCTGGCTGGAAGTTGTTCAGCATCGAGCGGGAGGACCTGCGAGCCGGTTGAGCTGCTGCCTGCTGCTGCTCGCGAGCTCCGGCGTAGTTGGCAAACCGGG  >  NC_002947/4408737‑4409028

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