Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,349,030	A→T	6.5%	T502S (ACC→TCC)	PP_4704 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,349,030	0	A	T	6.5%	171.7 / 2.7	62	T502S (ACC→TCC)	PP_4704	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/28);  new base T (2/2);  total (30/32)
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.25e-01

CAGGCCTACATCAGCCCGTTCATCCTGCTGGAGCTGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCTTCATCATCAACGAAGACCGCACCACCAAT  >  NC_002947/5348886‑5349176                                                                                                                                                  |                                                                                                                                                     cAGGCCTACATCAGCCCGTTCATCCTGCTGGAGCTGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCAt                                                                                                                                                  <  7:248248/149‑1 (MQ=255) cAGGCCTACATCAGCCCGTTCATCCTGCTGGAGCTGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCAt                                                                                                                                                  <  3:200638/149‑1 (MQ=255)            cAGCCCGTTCATCCTGCTGGAGCTGCGCAGCAGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGAAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCcgc                                                                                                                                        >  8:235213/1‑148 (MQ=255)              gCCCGTTCATCCTGCTGGAGCTGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGAc                                                                                                                                     <  2:270453/149‑1 (MQ=255)                ccGTTCATCCTGCTGGAGCTGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACtt                                                                                                                                   >  2:26166/1‑149 (MQ=255)                 cGTTCATCCTTCTGGAGCTGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAATAACACCGCGCC‑GCTGGCCTTCAATTTGACCGGCAAGGCGCAGGTCAGCCATTTGCATACCCTTTACACCATCAAGAGCCGCGACTTc                                                                                                                                  <  7:235213/149‑1 (MQ=255)                          tgctgGAGCTGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCA‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGACACTTGCAAACCCTCGACACCATAAAGAGCCGCGACTTCGTCAAATg                                                                                                                          >  4:290948/1‑148 (MQ=255)                                   tGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGt                                                                                                                 >  8:145370/1‑148 (MQ=255)                                            gCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGAc                                                                                                       <  2:239325/149‑1 (MQ=255)                                              atcctctCGAGCGACCTCAAGGTCGACCTGAAGAACACAGCGCCGGCTGCCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACACTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGAGAACCTCGACg                                                                                                      >  1:265327/4‑149 (MQ=255)                                                gCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCt                                                                                                   >  4:74689/1‑149 (MQ=255)                                                 ctctCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCtt                                                                                                  <  6:85044/149‑1 (MQ=255)                                                 ctctCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCt                                                                                                   <  3:74689/148‑1 (MQ=255)                                                 ctctCGAGCGACCTCAAGCTCGACCTGAACAACACCGCGCC‑GCTGGCCTTCAATCTGACCGGCAAGGCGCAGGTCAGCCAGTTCCATACCCTCGACCCCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCt                                                                                                   <  6:310792/148‑1 (MQ=255)                                                    tcGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTc                                                                                               >  1:61757/1‑149 (MQ=255)                                                       aGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTa                                                                                            >  3:2212/1‑149 (MQ=255)                                                        gCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTAc                                                                                           <  6:355778/149‑1 (MQ=255)                                                           aCCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTg                                                                                        >  6:49165/1‑149 (MQ=255)                                                             cTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCa                                                                                      >  3:16764/1‑149 (MQ=255)                                                              tCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCAc                                                                                     >  8:194397/1‑149 (MQ=255)                                                                   ctCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGa                                                                                <  6:91362/148‑1 (MQ=255)                                                                          tGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCTAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTg                                                                       >  1:120353/1‑149 (MQ=255)                                                                          tGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTg                                                                       >  5:45312/1‑149 (MQ=255)                                                                          tGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTg                                                                       >  6:176939/1‑149 (MQ=255)                                                                            aagaaCACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGt                                                                      <  6:45312/148‑1 (MQ=255)                                                                            aagaaCACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGt                                                                      <  5:176939/148‑1 (MQ=255)                                                                             agaaCACCGCGAC‑GCTGGCCTTATATGTGACCGGCAATGCGCAGGTCAGCCAGTTGCCTAGCCTCCACACCATCCCGAGCCGGGTCTTGGTCAAATGGCTGGAGGTGACCGTCGACGGCTTGTCCTACGTGCAGGGCGATGCT‑‑CTGtc                                                                     <  5:104458/148‑1 (MQ=255)                                                                              gaaCACCGCGCC‑GATGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTCCATACCCACGACAGCATCAAGAGCCGCGCCTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGACCTACGAGCACGGCGATGCT‑‑CTGtcga                                                                   >  4:25405/1‑149 (MQ=255)                                                                                acacCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGtcgatc                                                                 <  8:4486/149‑1 (MQ=255)                                                                                  acCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGtcgatcga                                                               >  2:159251/1‑149 (MQ=255)                                                                                  acCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGtcgatcga                                                               >  4:119722/1‑149 (MQ=255)                                                                                    cgcgCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCg                                                                                 <  3:180029/131‑1 (MQ=255)                                                                                    cgcgCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCg                                                                                 >  4:180029/1‑131 (MQ=255)                                                                                    cgcgCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGAc                                                              >  1:357963/1‑148 (MQ=255)                                                                                       gcc‑gcTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAg                                                          <  7:132436/149‑1 (MQ=255)                                                                                           gcTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGa                                                       <  2:61757/149‑1 (MQ=255)                                                                                             tGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGAcc                                                     >  2:106976/1‑149 (MQ=255)                                                                                               gCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGAcc                                                     <  1:26166/147‑1 (MQ=255)                                                                                                 cTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCctgc                                                 >  8:319920/1‑149 (MQ=255)                                                                                                  ttCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACTCCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGtcgatc                                                                 <  7:69536/132‑1 (MQ=255)                                                                                                  ttCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACTCCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGtcgatc                                                                 >  8:69536/1‑132 (MQ=255)                                                                                                  ttCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACTCCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGtcgatc                                                                 >  4:101264/1‑132 (MQ=255)                                                                                                  ttCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACTCCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGtcgatc                                                                 <  3:101264/132‑1 (MQ=255)                                                                                                  ttCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCctgct                                                >  5:255631/1‑149 (MQ=255)                                                                                                  ttCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCctgct                                                >  2:184230/1‑149 (MQ=255)                                                                                                  ttCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCctgct                                                >  8:180719/1‑149 (MQ=255)                                                                                                    cAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCctgct                                                <  4:16764/147‑1 (MQ=255)                                                                                                     aaTGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCa                                             <  4:2212/149‑1 (MQ=255)                                                                                                           cccGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCt                                        <  2:16829/147‑1 (MQ=255)                                                                                                           aCCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCt                                        <  2:22587/148‑1 (MQ=255)                                                                                                             cGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATg                                     <  2:357963/149‑1 (MQ=255)                                                                                                               gCAAGGCGCAGGTAAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCa                                   <  3:290948/149‑1 (MQ=255)                                                                                                                cAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGc                                    >  6:101082/1‑147 (MQ=255)                                                                                                                  aGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGc                                <  7:145370/149‑1 (MQ=255)                                                                                                                    gcgcAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCtt                              <  7:319920/149‑1 (MQ=255)                                                                                                                     cgcAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCtt                              >  3:280304/1‑148 (MQ=255)                                                                                                                      gcAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCTtca                            >  7:268258/1‑149 (MQ=255)                                                                                                                      gcAGGGCAGCCCGTTGCATACCCTCGACCCCATCAAGAGCCGCGCCTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCTtca                            <  4:268695/149‑1 (MQ=255)                                                                                                                        aGGTCAGCCAGTTGCATACCCGCGACACCATCAAGAGCCGCGACTTCTTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCATCCCTATGCACGCTtcatg                          <  2:265327/149‑2 (MQ=255)                                                                                                                              gCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCTTCATCATCAAc                    <  8:268258/149‑1 (MQ=255)                                                                                                                                       aTACCCTCGACACCATCAACAGCCGCGACTTCGTCAAATTGCAGCACGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATTCTCCCTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCTTCATCATCAACGAAGAc              <  1:106976/148‑1 (MQ=255)                                                                                                                                             tCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCTTCATCATCAACGAAGACCGcaccac      >  7:93389/1‑148 (MQ=255)                                                                                                                                                acacCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCTTCATCATCAACGAAGACCGCACCACCAAt  <  6:255631/149‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                     CAGGCCTACATCAGCCCGTTCATCCTGCTGGAGCTGCGCAGCGGCATGCTCTCGAGCGACCTCAAGGTCGACCTGAAGAACACCGCGCC‑GCTGGCCTTCAATGTGACCGGCAAGGCGCAGGTCAGCCAGTTGCATACCCTCGACACCATCAAGAGCCGCGACTTCGTCAAATGGCAGCAGGTGAACGTCGACGGCTTGTCCTACGTGCACGGCGATGCT‑‑CTGTCGATCGACAAAGTGACCCTGCTGCAGCCCTATGCACGCTTCATCATCAACGAAGACCGCACCACCAAT  >  NC_002947/5348886‑5349176

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