Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,287,681	T→A	8.8%	Y270N (TAC→AAC)	PP_2016 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,287,681	0	T	A	8.8%	190.2 / 6.6	68	Y270N (TAC→AAC)	PP_2016	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base A (4/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

TGAACTACCCGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTT  >  NC_002947/2287535‑2287828                                                                                                                                                      |                                                                                                                                                    tGAACTACCCGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTAc                                                                                                                                                   <  2:65854/149‑1 (MQ=255)    aCTACCCGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAg                                                                                                                                                <  4:352080/149‑1 (MQ=255)        cccGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCAACGACACCGGTGTGAGCAATATTCCCAAG‑‑T‑ACACGGGCTACCAGtggc                                                                                                                                            >  8:170182/1‑146 (MQ=255)                  gCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCa                                                                                                                                  <  6:24951/149‑1 (MQ=255)                  gCGTACCCGAAGAAGATTCGCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCa                                                                                                                                  <  6:265682/149‑1 (MQ=255)                    gTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCAACTGGTAGGCCGTgtacttgc                                                                                                                                <  1:187721/149‑8 (MQ=255)                      aCCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTa                                                                                                                               >  1:55255/1‑148 (MQ=255)                        ccGAAGTTGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                            <  6:83439/149‑1 (MQ=255)                          gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                          <  3:78292/149‑1 (MQ=255)                          gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                          >  8:135475/1‑149 (MQ=255)                          gaagaaGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATATTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAg                                                                                                                          >  3:228450/1‑149 (MQ=255)                                gaTTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATTTTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTg                                                                                                                    <  1:93231/149‑1 (MQ=255)                                  ttGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGca                                                                                                                  <  4:59878/149‑1 (MQ=255)                                        gCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCg                                                                                                            >  6:52924/1‑149 (MQ=255)                                           tCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTAc                                                                                                         >  2:304084/1‑149 (MQ=255)                                              aCCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCg                                                                                                       >  1:39/1‑148 (MQ=255)                                                cTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaa                                                                                                     <  7:21398/148‑1 (MQ=255)                                                cTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaa                                                                                                     >  8:21398/1‑148 (MQ=255)                                                    aCAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGGCCTGCACATGCGTACCGAACAAcc                                                                                                <  6:70983/149‑1 (MQ=255)                                                      aaGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                            <  7:63673/119‑1 (MQ=255)                                                      aaGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTaaaa                                                                                                                            >  8:63673/1‑119 (MQ=255)                                                        ggACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCaa                                                                                            <  7:135475/149‑1 (MQ=255)                                                                       gggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAg                                                                             <  6:178627/149‑1 (MQ=255)                                                                       gggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAg                                                                             <  5:179086/149‑1 (MQ=255)                                                                        ggCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAgg                                                                            <  8:110926/149‑1 (MQ=255)                                                                         gCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGa                                                                           >  1:22474/1‑149 (MQ=255)                                                                          cACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGAt                                                                          >  5:28537/1‑149 (MQ=255)                                                                           aCGCTGCTGGGCATCAAGGGCCATTAACTTATCTTCGACACTGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCATTTGGCCGTGCTCCAGTAAAAATGACCTGCACATTTGTACCGAACAACCGCAAGTGATc                                                                                      <  7:56088/136‑1 (MQ=255)                                                                           aCGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACTGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATc                                                                                      >  8:56088/1‑136 (MQ=255)                                                                           aCGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGAt                                                                          >  7:175791/1‑148 (MQ=255)                                                                                    ggCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGgc                                                                  <  6:365267/147‑1 (MQ=255)                                                                                      cATCAAGGGCCAGTACCTGATCTTCGACTCCTGTGTGATCAATATTAGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACCTGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCg                                                              <  4:228450/149‑1 (MQ=255)                                                                                      cATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCTCAAGCCTCTCACGTAGGGCGAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGgcg                                                                 <  1:312112/149‑1 (MQ=255)                                                                                       aTCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCg                                                              >  4:368142/1‑148 (MQ=255)                                                                                            gggCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTAc                                                        <  3:226510/149‑1 (MQ=255)                                                                                                cAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGa                                                    >  5:368626/1‑149 (MQ=255)                                                                                                  gTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATc                                                  <  6:368626/149‑1 (MQ=255)                                                                                                   tACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATc                                                  >  4:82699/1‑148 (MQ=255)                                                                                                      cTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgac                                               >  2:363634/1‑148 (MQ=255)                                                                                                       tGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaacaa                                                                                                  >  5:276356/1‑96 (MQ=255)                                                                                                       tGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGaacaa                                                                                                  <  6:276356/96‑1 (MQ=255)                                                                                                         aTCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgacga                                             <  5:52924/147‑1 (MQ=255)                                                                                                         aTCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATcgacga                                             <  2:55255/147‑1 (MQ=255)                                                                                                               gACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCAcc                                     >  2:79817/1‑149 (MQ=255)                                                                                                                 cacCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTCCACATCCGTACCGAACAACCGCAAGTGATCTACCTAAAGGATTACCAGGCGCCGGAGAACCTGATCAACGACACGCACCt                                    >  7:245953/1‑148 (MQ=255)                                                                                                                  acCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                  <  1:362525/149‑1 (MQ=255)                                                                                                                   ccGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGAc                                 <  3:182952/149‑1 (MQ=255)                                                                                                                    cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                  <  4:325095/147‑1 (MQ=255)                                                                                                                    cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGa                                  >  3:325095/1‑147 (MQ=255)                                                                                                                    cGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGAcc                                >  5:274340/1‑149 (MQ=255)                                                                                                                         tgATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCg                            <  3:297461/148‑1 (MQ=255)                                                                                                                         tgATCAATATACGCAAG‑‑T‑AGACGGCCTACCTTTTGGCCGTGCTCCAGTAATAAGGACCTGCACATGCGTACCGAACAACCGCTAGTGATCTACCTCATGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGa                           <  7:170182/149‑1 (MQ=255)                                                                                                                           aTCATTATTCGGAAG‑‑T‑ACTCGGTCTAGCAGTTGGCCGTGCTCCATTAAAAAGGAGCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGTATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAg                          <  8:245953/148‑1 (MQ=255)                                                                                                                              aaTATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGt                      <  3:368142/149‑1 (MQ=255)                                                                                                                                tatTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTc                    <  6:28537/149‑1 (MQ=255)                                                                                                                                    cGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAg                 <  8:175791/148‑1 (MQ=255)                                                                                                                                    cGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAg                 <  2:22474/148‑1 (MQ=255)                                                                                                                                      cAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGAc              >  5:84803/1‑149 (MQ=255)                                                                                                                                        aG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACgaga                                           >  1:179588/1‑118 (MQ=255)                                                                                                                                        aG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACgaga                                           <  2:179588/118‑1 (MQ=255)                                                                                                                                         ctggagcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGAc              >  2:30662/7‑149 (MQ=255)                                                                                                                                          tggagcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCAAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGATGTTCGAGGAc              >  1:109682/6‑148 (MQ=255)                                                                                                                                              gcacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATTCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCacac         >  8:143455/2‑149 (MQ=255)                                                                                                                                               cacGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGCCGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACg        >  6:242463/1‑149 (MQ=255)                                                                                                                                                acGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGa                  >  5:240056/1‑138 (MQ=255)                                                                                                                                                acGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGa                  <  6:240056/138‑1 (MQ=255)                                                                                                                                                acGGCCAACTGGTAGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGc       <  1:30662/149‑1 (MQ=255)                                                                                                                                                     cTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGtt  <  1:304084/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    TGAACTACCCGGTCGAGGCGTACCCGAAGAAGATTGTCAGCTTCAACCTCGACAAGGACCCGGTGGCCGAGGGCACGCTGCTGGGCATCAAGGGCCAGTACCTGATCTTCGACACCGGTGTGATCAATATTCGCAAG‑‑T‑ACACGGCCTACCAGTTGGCCGTGCTCCAGTAAAAAGGACCTGCACATGCGTACCGAACAACCGCAAGTGATCTACCTCAAGGATTACCAGGCGCCGGAGTACCTGATCGACGAGACGCACCTGACCTTCGAGCTGTTCGAGGACCACACGCTGGTT  >  NC_002947/2287535‑2287828

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑