Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,984,218	A→T	6.8%	L140Q (CTG→CAG)	PP_3512 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,984,218	0	A	T	6.8%	200.1 / 3.2	69	L140Q (CTG→CAG)	PP_3512	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (31/33);  new base T (2/3);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

CACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGCACCACACCACCAGG  >  NC_002947/3984072‑3984364                                                                                                                                                   |                                                                                                                                                   cacGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACcagc                                                                                                                                                  <  8:106892/149‑1 (MQ=255) cacGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACcag                                                                                                                                                   >  3:94951/1‑148 (MQ=255)   cGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACcagcag                                                                                                                                                <  6:61256/149‑1 (MQ=255)      cccccTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTg                                                                                                                                             >  5:225423/1‑149 (MQ=255)         ccTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCg                                                                                                                                          >  3:309832/1‑149 (MQ=255)            cGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGc                                                                                                                                       >  8:170432/1‑149 (MQ=255)             gTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCTGGGCt                                                                                                                                      <  8:153301/149‑1 (MQ=255)             gTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCt                                                                                                                                      <  4:295499/149‑1 (MQ=255)                 gAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTg                                                                                                                                   <  4:159222/148‑1 (MQ=255)                 gAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGt                                                                                                                                  <  3:162113/149‑1 (MQ=255)                    cGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCg                                                                                                                                >  1:298250/1‑148 (MQ=255)                    cGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGa                                                                                                                               <  4:94951/149‑1 (MQ=255)                       gTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACg                                                                                                                            <  4:178256/149‑1 (MQ=255)                             ccGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGt                                                                                                                      >  3:285450/1‑149 (MQ=255)                              cgcgACTAGGCCACACACGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCTCATTAGCGGGCAAGCCCGCGCCGGGTGACAGATTCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAg                                                                                                                       <  8:154383/147‑1 (MQ=255)                              cgcgAATAGGCCACACCAGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGAGGGCAAGCACGCGCCGGGTGACAGAATCATGCCACACGCCGCTGAACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGta                                                                                                                     >  7:154383/1‑148 (MQ=255)                                    tAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCTGCTGGTGcagcg                                                                                                                                        >  6:78806/1‑120 (MQ=255)                                    tAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCTGCTGGTGcagcg                                                                                                                                        <  7:92343/124‑5 (MQ=255)                                    tAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCTGCTGGTGcagcg                                                                                                                                        >  8:92343/1‑120 (MQ=255)                                    tAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCTGCTGGTGcagcg                                                                                                                                        <  5:78806/124‑5 (MQ=255)                                    tAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGAc                                                                                                               >  2:17941/1‑149 (MQ=255)                                       gCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAgg                                                                                                                                               >  1:306463/1‑114 (MQ=255)                                       gCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAgg                                                                                                                                               <  2:306463/114‑1 (MQ=255)                                         cacaCCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCAGCGCCGGGTGACAGAATCATGCAACACGCCGCTGCACCAGCAGGTGCCGTGCTTTGTCGAACGCCCAGTTGTAGACCACgg                                                                                                          >  8:288190/1‑149 (MQ=255)                                              ccGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAg                                                                                                     >  5:320965/1‑149 (MQ=255)                                                 gcatccaATAATTCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGc                                                                                                  <  6:225423/145‑1 (MQ=255)                                                 ccagccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAggg                                                                                                   >  8:18377/1‑148 (MQ=255)                                                  ctgccaATAATGCCTGGGATATCACCGTTTTCAGTGCAGCCTCATTAGCGGGAAAGCCCGCGCCGGGTGACAGAATCATGCCACATGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACTCCCTGTTGTAGACCACGGTGTAGGGCa                                                                                                 <  7:202274/147‑1 (MQ=255)                                                  cagccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCa                                                                                                 <  3:314461/149‑1 (MQ=255)                                                     ccaATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGga                                                                                              <  5:290665/149‑1 (MQ=255)                                                        ataatGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGgatga                                                                                           >  6:297329/1‑149 (MQ=255)                                                           atGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATc                                                                                         <  2:100065/148‑1 (MQ=255)                                                               cTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACgaa                                                                                    <  1:17941/149‑1 (MQ=255)                                                                               gTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACgaagaa                                                                                 <  7:227802/136‑1 (MQ=255)                                                                               gTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACgaagaa                                                                                 >  8:227802/1‑136 (MQ=255)                                                                                  aGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCAt                                                                  >  1:118356/1‑148 (MQ=255)                                                                                      cAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCTGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAAt                                                              <  7:170432/148‑1 (MQ=255)                                                                                       aGCCGCATTAGCGGGCAAGCCAGCGACGGGTGACAGAATCATGCCACACGAAGCTGCACCAGCAGGTGCCGGGGTTTGTCGAACGCACAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGATGAAGCCCAGCTCACCCATGAATGc                                                            >  8:49009/1‑149 (MQ=255)                                                                                         ccGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCt                                                          >  6:267034/1‑149 (MQ=255)                                                                                         ccGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCt                                                          >  7:78126/1‑149 (MQ=255)                                                                                          cGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTg                                                         <  4:309832/149‑1 (MQ=255)                                                                                           gCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGa                                                                <  8:247963/141‑1 (MQ=255)                                                                                           gCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGa                                                                >  7:247963/1‑141 (MQ=255)                                                                                            cATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTgca                                                       >  2:321645/1‑149 (MQ=255)                                                                                            cATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTgca                                                       <  4:285450/149‑1 (MQ=255)                                                                                             aTTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTgcag                                                      <  6:320965/149‑1 (MQ=255)                                                                                             aTTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGGGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTgcag                                                      <  7:18377/149‑1 (MQ=255)                                                                                                 gCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCagcagc                                                  <  5:267034/149‑1 (MQ=255)                                                                                                   gggCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCg                                                 >  1:85992/1‑148 (MQ=255)                                                                                                         gCCAGCGCCGGGTGACAGAATCATGCCACACGCCTCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGTATGCCTTCAGCAGCTAAAtctc                                          <  7:49009/149‑1 (MQ=255)                                                                                                          cccGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAAtctc                                          >  3:255206/1‑148 (MQ=255)                                                                                                           ccGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACg                                                                                      >  8:6895/1‑103 (MQ=255)                                                                                                           ccGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACg                                                                                      <  7:6895/103‑1 (MQ=255)                                                                                                           ccGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAAtctc                                          >  7:30109/1‑147 (MQ=255)                                                                                                               gcCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCAt                                    >  8:214512/1‑149 (MQ=255)                                                                                                                   ggTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCaa                                <  2:298250/149‑1 (MQ=255)                                                                                                                    gTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGTGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAg                               <  7:288190/149‑1 (MQ=255)                                                                                                                     tgccAGAATCATGCCACACGCCGCTGCACCAGCACGTGCCGGGCTATGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGCATGATCACGAAGAAGCCCAGCTCCACCATCAATGCCTGCagcagc                                                  >  2:319809/4‑129 (MQ=255)                                                                                                                     tGACATAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGGAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCagcagc                                                  <  1:319809/129‑1 (MQ=255)                                                                                                                       aCAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAgcc                             <  7:214512/148‑1 (MQ=255)                                                                                                                        cAGAATCATGCCACGCGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAgccgc                           >  7:274543/1‑149 (MQ=255)                                                                                                                            aTCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCAcc                       >  6:278846/1‑149 (MQ=255)                                                                                                                               aTGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCggc                    <  1:321645/149‑1 (MQ=255)                                                                                                                                  ccACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACTGTGTAGGGc                                                                                                  <  5:31852/68‑1 (MQ=255)                                                                                                                                  ccACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGc                                                                                                  >  6:31852/1‑68 (MQ=255)                                                                                                                                   cacaCGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCggcaggc                <  8:30109/149‑1 (MQ=255)                                                                                                                                     cacGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGcac              <  2:85992/149‑1 (MQ=255)                                                                                                                                        gccgcTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGcacca            >  4:279225/1‑148 (MQ=255)                                                                                                                                        gccgcTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGcacca            >  2:189898/1‑148 (MQ=255)                                                                                                                                          cgcTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGcaccacac         <  5:297329/149‑1 (MQ=255)                                                                                                                                              gCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGCaccacaccac      >  1:286088/1‑148 (MQ=255)                                                                                                                                                 ccAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGCACCACACCACCAgg  <  1:189898/149‑1 (MQ=255)                                                                                                                                                  cagcagGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCggcaggc                <  5:64595/134‑1 (MQ=255)                                                                                                                                                  cagcagGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCggcaggc                >  6:64595/1‑134 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACAGAATCATGCCACACGCCGCTGCACCAGCAGGTGCCGGGCTTTGTCGAACGCCCAGTTGTAGACCACGGTGTAGGGCAGGATGATCACGAAGAAGCCCAGCTCCACCATGAATGCCTGCAGCAGCGAAATCTCCAGCATCCAAGCCGCCACCGGCAGGCACCACACCACCAGG  >  NC_002947/3984072‑3984364

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