Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,848,339	T→A	13.3%	F35Y (TTC→TAC) ‡	PP_3397 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,848,339	0	T	A	13.3%	188.2 / 17.2	75	F35Y (TTC→TAC) ‡	PP_3397	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (36/29);  new base A (4/6);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 5.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.83e-01

CAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGGCGC  >  NC_002947/3848197‑3848485                                                                                                                                               |                                                                                                                                                   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                                                                                                                                          >  3:90329/1‑120 (MQ=255)                                  ggtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCGGCGAACGCGGCGCTg                                                                                                                                          <  4:90329/120‑1 (MQ=255)                                   gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCtt                                                                                                             <  2:211278/148‑1 (MQ=255)                                   gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                            <  5:90262/149‑1 (MQ=255)                                   gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                            >  1:62036/1‑149 (MQ=255)                                   gtgcCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACATGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTc                                                                                                            >  1:270527/1‑149 (MQ=255)                                        aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTc                                                                                                                                                   <  6:207611/105‑1 (MQ=255)                                        aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTc                                                                                                                                                   >  5:207611/1‑105 (MQ=255)                                        aTGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCca                                                                                                        >  4:152420/1‑148 (MQ=255)                                         tGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTac                                                                                                                  >  5:174255/1‑137 (MQ=255)                                         tGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTac                                                                                                                  <  6:174255/137‑1 (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                                                                                     <  5:142809/149‑1 (MQ=255)                                                          aGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGAtgtg                                                                                     <  4:299439/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                                                                   <  4:336743/137‑1 (MQ=255)                                                                                        tCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGa                                                                   >  3:336743/1‑137 (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                                              <  4:36738/128‑1 (MQ=255)                                                                                                                      tctctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCa                                              >  3:36738/1‑128 (MQ=255)                                                                                                                        tctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTg                       >  7:345519/1‑149 (MQ=255)                                                                                                                         ctTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTgc                      >  8:286993/1‑149 (MQ=255)                                                                                                                          ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGt                                       >  7:207261/1‑131 (MQ=255)                                                                                                                          ttACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGt                                       <  8:207261/131‑1 (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                                                               >  3:77112/1‑97 (MQ=255)                                                                                                                                    cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCaa                                                               <  4:77112/97‑1 (MQ=255)                                                                                                                                    cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGTCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAc           >  7:276039/1‑149 (MQ=255)                                                                                                                                    cAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGCGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCAc           <  3:285654/149‑1 (MQ=255)                                                                                                                                             gTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGgcgc  <  7:286993/149‑1 (MQ=255)                                                                                                                                               |                                                                                                                                                   CAACGGTGGTAGCGACTCCGAAGACGTGGTGCTGGTGCCATGAACAGACAAAGGCTCAGCCTTCAAGGAGGCAGCATGAACACTTGGTCGCGCCGTGCGCTGTTGGCCGCCGGTTTCTCTCTTACCGTTACCAGCGCCGCGTTCGCCGCGCTGTCTGATGTCGATCCGGGGCCCTACACCTTCGCCACAGGTGGCTATCCGATGTGGTACAAGGACAGCCTCGACCAATCCCTGGAGCTGTGCCAGTCGCGTGCCACCAGCACCCGTGCGCCAGGCGCACCCGGGGCGC  >  NC_002947/3848197‑3848485

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