Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,848,334	C→G	13.5%	A33A (GCC→GCG)	PP_3397 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,848,334	0	C	G	13.5%	204.0 / 20.4	74	A33A (GCC→GCG)	PP_3397	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (36/28);  new base G (4/6);  total (40/34)
Fisher's exact test for biased strand distribution p-value = 4.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.33e-01

GTCACCAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAC  >  NC_002947/3848186‑3848476                                                                                                                                                     |                                                                                                                                               gTCACCAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGcgcc                                                                                                                                                >  4:162001/1‑149 (MQ=255)       aGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGATCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGATACTCTCTTACCGTTACAAGCgcggcgaac                                                                                                                                          >  2:336000/1‑146 (MQ=255)            cAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTcgccgc                                                                                                                                     >  6:41686/1‑149 (MQ=255)                gtttgtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAACGTCTCAGCCTTCAAGGATGCAGCATGAACACTTGGTCGCGCCGTGCGCTGGTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                 <  3:252847/145‑1 (MQ=255)                ggtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                 <  8:5063/149‑1 (MQ=255)                ggtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                 <  1:336000/149‑1 (MQ=255)                ggtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACATGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                 <  3:247923/149‑1 (MQ=255)                 gtggtAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGt                                                                                                                                >  1:234193/1‑149 (MQ=255)                         gACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTc                                                                                                                        >  2:192185/1‑149 (MQ=255)                          aCTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTc                                                                                                                        >  1:319815/1‑148 (MQ=255)                           cTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGa                                                                                                                      <  2:234193/149‑1 (MQ=255)                              cGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATc                                                                                                                    <  7:279957/148‑1 (MQ=255)                              cGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTGTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATcc                                                                                                                   >  4:285654/1‑149 (MQ=255)                                  gACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCgggg                                                                                                               <  3:162001/149‑1 (MQ=255)                                      tggtgctggtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCt                                                                                                           <  2:319815/149‑1 (MQ=255)                                             ggtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                 <  4:90329/120‑1 (MQ=255)                                             ggtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                 >  3:90329/1‑120 (MQ=255)                                              gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCtt                                                                                                    <  2:211278/148‑1 (MQ=255)                                              gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                   >  1:62036/1‑149 (MQ=255)                                              gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                   <  5:90262/149‑1 (MQ=255)                                              gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACATGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                   >  1:270527/1‑149 (MQ=255)                                                   aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTc                                                                                                                                          >  5:207611/1‑105 (MQ=255)                                                   aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTc                                                                                                                                          <  6:207611/105‑1 (MQ=255)                                                   aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCca                                                                                               >  4:152420/1‑148 (MQ=255)                                                    tGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTac                                                                                                         <  6:174255/137‑1 (MQ=255)                                                    tGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTac                                                                                                         >  5:174255/1‑137 (MQ=255)                                                     gAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAg                                                                                            >  6:102705/1‑149 (MQ=255)                                                     gAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAg                                                                                            >  1:248659/1‑149 (MQ=255)                                                      aaCAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAgg                                                                                           >  1:318968/1‑149 (MQ=255)                                                         aGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGCTCCGGGGCCCTACACCTTCGCCACAGGTgg                                                                                        >  8:199004/1‑149 (MQ=255)                                                         aGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTgg                                                                                        <  1:192185/149‑1 (MQ=255)                                                           aCAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGc                                                                                       >  5:67359/1‑148 (MQ=255)                                                                     aGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGAtgtg                                                                            <  4:299439/149‑1 (MQ=255)                                                                     aGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGAtgtg                                                                            <  5:142809/149‑1 (MQ=255)                                                                      gCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTgg                                                                           >  2:362600/1‑149 (MQ=255)                                                                      gCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTgg                                                                           <  2:318968/149‑1 (MQ=255)                                                                            aaGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAg                                                                     <  2:248659/149‑1 (MQ=255)                                                                              ggaggCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGa                                                                   <  5:41686/149‑1 (MQ=255)                                                                                    agcaTGAACGCTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCt                                                             >  4:336091/1‑149 (MQ=255)                                                                                    agcaTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCt                                                             >  1:171250/1‑149 (MQ=255)                                                                                    agcaTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCt                                                             >  5:19853/1‑149 (MQ=255)                                                                                     gcaTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGCTGTGGTACAAGGACAGCCTc                                                            <  6:67359/149‑1 (MQ=255)                                                                                      caTGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCg                                                           <  2:62036/149‑1 (MQ=255)                                                                                             acTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATc                                                    >  8:46321/1‑149 (MQ=255)                                                                                                 ggTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATccc                                                  <  8:254738/147‑1 (MQ=255)                                                                                                  gTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTgg                                               <  3:152420/149‑1 (MQ=255)                                                                                                   tCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGa                                                          >  3:336743/1‑137 (MQ=255)                                                                                                   tCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGa                                                          <  4:336743/137‑1 (MQ=255)                                                                                                   tCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGa                                              <  1:137679/149‑1 (MQ=255)                                                                                                     gcgcCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGc                                            >  3:169435/1‑149 (MQ=255)                                                                                                      cgcCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCt                                           >  3:351136/1‑149 (MQ=255)                                                                                                             cgcTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAg                                    >  8:341517/1‑149 (MQ=255)                                                                                                                    ggCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTg                             >  4:257915/1‑149 (MQ=255)                                                                                                                      ccgccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGc                            >  8:223650/1‑148 (MQ=255)                                                                                                                        gccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                     <  3:68786/137‑1 (MQ=255)                                                                                                                        gccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccac                         >  4:4849/1‑149 (MQ=255)                                                                                                                        gccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccac                         >  3:308331/1‑149 (MQ=255)                                                                                                                        gccgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACATTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                     >  4:68786/1‑137 (MQ=255)                                                                                                                          cgGTTTCTCTCTTACCGTTACCAGCGCCGGGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccacca                       <  2:270527/149‑1 (MQ=255)                                                                                                                          cgGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccacca                       <  1:362600/149‑1 (MQ=255)                                                                                                                                 tctctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                     >  3:36738/1‑128 (MQ=255)                                                                                                                                 tctctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                     <  4:36738/128‑1 (MQ=255)                                                                                                                                   tctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTg              >  7:345519/1‑149 (MQ=255)                                                                                                                                    ctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTgc             >  8:286993/1‑149 (MQ=255)                                                                                                                                     ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGt                              <  8:207261/131‑1 (MQ=255)                                                                                                                                     ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGt                              >  7:207261/1‑131 (MQ=255)                                                                                                                                        ccGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTgcgc           >  7:29122/1‑147 (MQ=255)                                                                                                                                         atcagacAGCGCGGCGAACGCGGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAg        >  1:137333/7‑149 (MQ=255)                                                                                                                                          gTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAgg       <  8:345519/149‑1 (MQ=255)                                                                                                                                              acAGCGCGGCGAACGCGGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccac                         <  6:22684/126‑1 (MQ=255)                                                                                                                                              acAGCGCGGCGAACGCGGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGccac                         >  5:22684/2‑127 (MQ=255)                                                                                                                                              acAGCGCGGCGAACGCGGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCa   >  2:120691/2‑149 (MQ=255)                                                                                                                                               cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCaa                                                      <  4:77112/97‑1 (MQ=255)                                                                                                                                               cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCaa                                                      >  3:77112/1‑97 (MQ=255)                                                                                                                                               cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGTCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAc  >  7:276039/1‑149 (MQ=255)                                                                                                                                               cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGCGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAc  <  3:285654/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               GTCACCAGCGCCAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAC  >  NC_002947/3848186‑3848476

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