Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,832,307	T→A	5.8%	N295I (AAT→ATT)	pvdY ←	hydroxyproline acetylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,832,307	0	T	A	5.8%	213.2 / 2.4	69	N295I (AAT→ATT)	pvdY	hydroxyproline acetylase
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base A (2/2);  total (33/36)
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.56e-01

AACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGT  >  NC_002947/4832161‑4832454                                                                                                                                                    |                                                                                                                                                        aaCACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCAtt                                                                                                                                                        <  4:125239/148‑1 (MQ=255)    accacACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCg                                                                                                                                                    >  6:51488/1‑149 (MQ=255)    accacACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGGCGACCTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGtag                                                                                                                                                    >  7:73018/1‑147 (MQ=255)      cacaCTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGg‑c                                                                                                                                                 >  2:197257/1‑149 (MQ=255)        cacTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGg‑cgc                                                                                                                                               <  4:242614/149‑1 (MQ=255)         acTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGg‑cgcg                                                                                                                                              >  1:140062/1‑149 (MQ=255)           tAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑cgcgcg                                                                                                                                            >  8:289834/1‑149 (MQ=255)           tAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑cgcgcg                                                                                                                                            >  3:274310/1‑149 (MQ=255)             gCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCgg                                                                                                                                           <  3:17290/148‑1 (MQ=255)               cccTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTc                                                                                                                                        <  2:115170/149‑1 (MQ=255)                ccTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCATGTCTGGAAAGTCGAACACCTTGTCCCAAT‑GGAAACATTCGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑AGAGCGGCTCg                                                                                                                                       >  8:255737/1‑149 (MQ=255)                     gTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGccac                                                                                                                                  <  5:113023/149‑1 (MQ=255)                           ttGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCg                                                                                                                            >  8:208704/1‑149 (MQ=255)                             gCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCt                                                                                                                          >  4:311365/1‑149 (MQ=255)                               aaCGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCTCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTgg                                                                                                                        <  6:57331/149‑1 (MQ=255)                               aaCGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTgg                                                                                                                        <  2:62814/149‑1 (MQ=255)                                  gATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTg                                                                                                                     >  5:14290/1‑149 (MQ=255)                                    tCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTgcg                                                                                                                   >  1:148573/1‑149 (MQ=255)                                      gAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGcgcg                                                                                                                 >  2:224446/1‑149 (MQ=255)                                      gAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCAATGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGcgcg                                                                                                                 >  2:213366/1‑149 (MQ=255)                                       aaaaaCCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCgg                                                                                                                <  5:51488/149‑1 (MQ=255)                                        aaaaCCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCgg‑‑g                                                                                                             <  6:14290/149‑1 (MQ=255)                                         aaaCCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑Gt                                                                                                            <  2:148573/149‑1 (MQ=255)                                             cGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATc                                                                                                        >  3:84481/1‑149 (MQ=255)                                                  ccGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGGCCGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAg                                                                                                     <  1:224446/148‑1 (MQ=255)                                                  ccGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGa                                                                                                   <  7:289834/149‑1 (MQ=255)                                                    gCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGaaa                                                                                                 <  1:197257/149‑1 (MQ=255)                                                    gCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCAATGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCa                                                                                                      <  3:261467/144‑1 (MQ=255)                                                    gCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCAATGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGAGCGCGG‑‑GTCATCCa                                                                                                      >  4:261467/1‑144 (MQ=255)                                                      cccAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAg                                                                                               <  1:114187/149‑1 (MQ=255)                                                       ccAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAgg                                                                                              <  3:311365/149‑1 (MQ=255)                                                        cAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAgg                                                                                              >  2:323405/1‑148 (MQ=255)                                                           aaTCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTTGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGt                                                                                          >  4:206048/1‑149 (MQ=255)                                                                  aaCGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATcc                                                                                                       >  3:161786/1‑129 (MQ=255)                                                                  aaCGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATcc                                                                                                       <  4:161786/129‑1 (MQ=255)                                                                  aaCGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCaa                                                                                  >  3:125381/1‑149 (MQ=255)                                                                    cgccgcCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACg                                                                                >  6:107717/1‑149 (MQ=255)                                                                        gcccgcTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGa                                                                            <  7:208704/149‑1 (MQ=255)                                                                        gcccgcTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGa                                                                            <  1:323405/149‑1 (MQ=255)                                                                         cccgcTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGa                                                                            >  2:266481/1‑148 (MQ=255)                                                                          ccggTTGTGGGGAAAGTCGAACTCCTTGTCGCTAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACGAACGCCGATa                                                                          <  7:255737/149‑1 (MQ=255)                                                                          ccgcTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCTTCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATa                                                                          <  1:244109/149‑1 (MQ=255)                                                                          ccgcTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATa                                                                          <  3:196280/149‑1 (MQ=255)                                                                            gcTTGTGTGGAAAGTCGAACTCCTTGTTGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGc                                                                        <  3:206048/149‑1 (MQ=255)                                                                                gtgtgGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGc                                                                    >  3:149867/1‑149 (MQ=255)                                                                                 tgtgGAAAGTGGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATGTTGCCATTGTGGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGc                                                                    <  8:73018/148‑1 (MQ=255)                                                                                  gtgGAAAGTCGAACTCCTTGTCGCAATGGGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑Ac                                                                                     >  1:177492/1‑131 (MQ=255)                                                                                  gtgGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑Ac                                                                                     <  2:177492/130‑1 (MQ=255)                                                                                     gAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGc                                                                                        >  3:122266/1‑125 (MQ=255)                                                                                     gAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGc                                                                                        <  4:122266/125‑1 (MQ=255)                                                                                     gAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGc                                                                >  3:169146/1‑148 (MQ=255)                                                                                              actactTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGCTACCAACGCCGATAGCCAGCTGGCCACCTtgtg                                                       <  1:39530/149‑1 (MQ=255)                                                                                              aCTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGc                                                      <  1:266481/149‑1 (MQ=255)                                                                                                  cTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCgg                                                    <  4:149867/147‑1 (MQ=255)                                                                                                         cAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTa                                                                                            <  6:247669/101‑1 (MQ=255)                                                                                                         cAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTa                                                                                            >  5:247669/1‑101 (MQ=255)                                                                                                          aaT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTgg                                          <  5:107717/149‑1 (MQ=255)                                                                                                                  aCATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGctctcct                                   <  4:169146/149‑1 (MQ=255)                                                                                                                        gTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCcca                             >  1:132540/1‑149 (MQ=255)                                                                                                                         tCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCccac                            <  2:140062/149‑1 (MQ=255)                                                                                                                         tCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCccac                            <  4:84481/149‑1 (MQ=255)                                                                                                                                   aGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATg                   <  4:125381/148‑1 (MQ=255)                                                                                                                                    gCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGa                 >  1:119567/1‑149 (MQ=255)                                                                                                                                     ccAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGGGTGTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATg                   >  1:89049/1‑148 (MQ=255)                                                                                                                                     ccAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGa                 >  4:50012/1‑148 (MQ=255)                                                                                                                                     ccAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGAt                <  6:53697/149‑1 (MQ=255)                                                                                                                                      cAATCATCTTGCCATTGTCAG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGAt                >  4:268905/1‑148 (MQ=255)                                                                                                                                       aaTCATCTTGCCAATGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGATAc              <  1:213366/149‑1 (MQ=255)                                                                                                                                                gCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGATACCACGGTCg      >  5:119281/1‑148 (MQ=255)                                                                                                                                                   aTTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGt  >  1:159926/1‑149 (MQ=255)                                                                                                                                                    ttGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAATGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGctctcctc                                  <  6:125172/116‑1 (MQ=255)                                                                                                                                                    ttGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGctctcctc                                  >  5:125172/1‑116 (MQ=255)                                                                                                                                                    |                                                                                                                                                        AACACCACACTAGCCCCTACGTCAGCTTGCAACGATCGAAAAACCGCTCCCGCCCCAGAATCATCAACGCCGCCCGCTTGTGTGGAAAGTCGAACTCCTTGTCGCAAT‑GGAAACATTGGTCGTGCATGTAGCCAATCATCTTGCCATTGTCGG‑CGCGCGGCTCGGCCACCACCCGCTGGGTGCGCGG‑‑GTCATCCAGGAAAAGGTAGTGC‑ACCAACGCCGATAGCCAGCTGGCCACCTTGTGCGGCCCTCGGTGGCTCTCCTCTCCCACCAGCATATGGATACCACGGTCGTAGT  >  NC_002947/4832161‑4832454

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