Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,637,742	A→C	6.5%	I108L (ATC→CTC)	PP_3204 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,637,742	0	A	C	6.5%	190.5 / 2.9	61	I108L (ATC→CTC)	PP_3204	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (25/32);  new base C (2/2);  total (27/34)
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.98e-01

CCCTGGACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACGAGGCGACAT  >  NC_002947/3637596‑3637890                                                                                                                                                   |                                                                                                                                                     cccTGGACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATc                                                                                                                                                    <  7:32278/149‑1 (MQ=255)   cTGGACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCgg                                                                                                                                                  <  5:156539/149‑1 (MQ=255)     ggACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggc                                                                                                                                                 >  3:423301/1‑148 (MQ=255)      gACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcga                                                                                                                                               <  1:409553/149‑1 (MQ=255)         aaCGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcgaggc                                                                                                                                            <  5:435871/149‑1 (MQ=255)         aaCGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcgaggc                                                                                                                                            <  6:397983/149‑1 (MQ=255)         aaCGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcgaggc                                                                                                                                            >  6:435871/1‑149 (MQ=255)           cGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCggcgaggc                                                                                                                                            >  3:174770/1‑147 (MQ=255)            ggCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGATCTCggcgaagacag                                                                                                                                         <  1:82682/149‑7 (MQ=255)             gCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGt                                                                                                                                        >  5:228202/1‑149 (MQ=255)             gCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCCTCGCCGATCTCggcgaagacagc                                                                                                                                        <  8:366267/149‑8 (MQ=255)                     ccGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCAc                                                                                                                                >  8:248828/1‑149 (MQ=255)                       gACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACgg                                                                                                                              <  5:263607/149‑1 (MQ=255)                        aCAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACAGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCTGCACggg                                                                                                                             >  4:267203/1‑149 (MQ=255)                          aGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGc                                                                                                                            >  2:203597/1‑148 (MQ=255)                            cgcAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCaaa                                                                                                                         <  4:423301/149‑1 (MQ=255)                              cAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGc                                                                                                                       <  2:28307/149‑1 (MQ=255)                                gCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGcc                                                                                                                      >  6:190335/1‑148 (MQ=255)                                    cTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTc                                                                                                                 <  7:248828/149‑1 (MQ=255)                                       cAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGGTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGaa                                                                                                              <  8:185065/149‑1 (MQ=255)                                       cAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGaa                                                                                                              <  7:252211/149‑1 (MQ=255)                                               cccTGGCCCGCTGTGCATCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGc                                                                                                      <  6:340445/149‑1 (MQ=255)                                                 cTGGCCCGCGGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTg                                                                                                    <  5:190335/149‑1 (MQ=255)                                                    gCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGAt                                                                                                  <  4:174770/148‑1 (MQ=255)                                                       cGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATgcac                                                                                              >  8:354971/1‑149 (MQ=255)                                                       cGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATgcac                                                                                              >  6:224900/1‑149 (MQ=255)                                                          tGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATgcaccgc                                                                                           <  1:370928/149‑1 (MQ=255)                                                                  gTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAAc                                                                                   >  7:106145/1‑149 (MQ=255)                                                                    gATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAAccc                                                                                 >  1:32292/1‑149 (MQ=255)                                                                     aTTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAAccc                                                                                 <  8:106145/148‑1 (MQ=255)                                                                      ttCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGTACCGCACCGAAACCCTg                                                                               <  3:322067/149‑1 (MQ=255)                                                                        cgcgTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTgg                                                                              <  7:354971/148‑1 (MQ=255)                                                                         gcgTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGAc                                                                            >  8:215423/1‑149 (MQ=255)                                                                               tggACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGc                                                                      <  1:203597/149‑1 (MQ=255)                                                                                  aCATACCGCCCGACTGCGTGCAGATCCAGGTCAGTCCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATc                                                                   <  2:196427/149‑1 (MQ=255)                                                                                      aCCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACACGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTa                                                               >  3:63186/1‑149 (MQ=255)                                                                                        cGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTAcc                                                             >  4:209270/1‑149 (MQ=255)                                                                                        cGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATACAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTAcc                                                             >  1:409017/1‑149 (MQ=255)                                                                                           ccGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGATATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAg                                                          <  7:215423/149‑1 (MQ=255)                                                                                                   gTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGt                                                  >  4:312550/1‑149 (MQ=255)                                                                                                     gCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTct                                                >  6:167872/1‑149 (MQ=255)                                                                                                     gCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTct                                                >  3:120827/1‑149 (MQ=255)                                                                                                      cAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTctg                                               >  1:175847/1‑149 (MQ=255)                                                                                                      cAGAACCAGGTCAGTGCCGAAGATACAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTctg                                               <  2:409017/149‑1 (MQ=255)                                                                                                          aCCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCg                                                      <  8:375131/138‑1 (MQ=255)                                                                                                          aCCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCg                                                      >  7:375131/1‑138 (MQ=255)                                                                                                           ccAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTctggctgg                                          >  1:352287/1‑149 (MQ=255)                                                                                                            cAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGt                                         <  3:209270/149‑1 (MQ=255)                                                                                                              ggTCAGTGCCGAAGATGCATCCTCTGTCTTCGCCGAGATCGGCGATGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGc                                       <  6:349975/149‑1 (MQ=255)                                                                                                                tCATTGCCGAACAGGCAGCCGCTGTCTTCGCCGAGATCGCCGAGGCGAGTGCCGGCACTGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACCTCATCGTTACCGAGGGCGAGGTCTGGCTTGTGCTg                                     <  6:228202/149‑1 (MQ=255)                                                                                                                tCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTg                                     <  2:32292/149‑1 (MQ=255)                                                                                                                    tGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGAc                                  >  7:304413/1‑148 (MQ=255)                                                                                                                    tGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGGTACCGAGGGCGAGGTCTGGCTGGTGCTGGACg                                 <  3:267203/149‑1 (MQ=255)                                                                                                                         aaGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGc                            <  3:12492/149‑1 (MQ=255)                                                                                                                              ctgcacGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAg                        <  4:294673/143‑1 (MQ=255)                                                                                                                                  ccGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCAc                   <  3:312550/149‑1 (MQ=255)                                                                                                                                      tcgcctCGCCGATCTCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTTCCGAGGGCGAGGTCTGGCTGGTGCTGGCCGATGGCGAAGTCCACCTc                >  5:213777/6‑148 (MQ=255)                                                                                                                                        tCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCaa              <  2:352287/148‑1 (MQ=255)                                                                                                                                          ctCGCCGATCTCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAAc            >  2:425832/2‑148 (MQ=255)                                                                                                                                           tCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAAc            >  7:230099/1‑147 (MQ=255)                                                                                                                                              ccGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACGAGGc       >  3:392999/1‑149 (MQ=255)                                                                                                                                                   aTCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACGAGGCGACAt  >  7:257377/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CCCTGGACAACGGCGATGATCCGACAGCGCAGCCGCTGCAATTGAGCCCTGGCCCGCTGGGCAGCGTGATTCGCGTGGTGGACATACCGCCCGACAGCGTGCAGAACCAGGTCAGTGCCGAAGATGCAGCCGCTGTCTTCGCCGAGATCGGCGAGGCGAGTGCCGGCACGGGCAAAGCCGATTCGAAGCACAAGCTGATGCACCGCACCGAAACCCTGGACTACGGCATCGTTACCGAGGGCGAGGTCTGGCTGGTGCTGGACGATGGCGAAGTCCACCTCAAACGAGGCGACAT  >  NC_002947/3637596‑3637890

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