Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,297,474 C→A 7.4% R151L (CGG→CTG)  PP_4671 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,297,4740CA 7.4% 253.9 / 8.6 81R151L (CGG→CTG) PP_4671hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/37);  new base A (2/4);  total (40/41)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.90e-01

GTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGCATCGCAGGCCTGGGCGACCT  >  NC_002947/5297328‑5297621
                                                                                                                                                  |                                                                                                                                                   
gTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCg                                                                                                                                                    >  5:252005/1‑148 (MQ=255)
   gggCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcagg                                                                                                                                                >  2:223808/1‑149 (MQ=255)
   gggCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                                                <  1:71545/149‑1 (MQ=255)
     gcgcACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggca                                                                                                                                              <  5:124752/149‑1 (MQ=255)
         aCGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAAc                                                                                                                                          <  7:90839/149‑1 (MQ=255)
         aCGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcaggtgacct                                                                                                                                          <  8:131444/149‑7 (MQ=255)
               cGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTg                                                                                                                                    <  3:98248/149‑1 (MQ=255)
                gCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGc                                                                                                                                   >  6:52628/1‑149 (MQ=255)
                  ttGCGGGTC‑GCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGa                                                                                                                                >  5:40306/1‑149 (MQ=255)
                     cGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGAcc                                                                                                                              <  6:369271/149‑1 (MQ=255)
                       ggTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGt                                                                                                                            <  2:340828/149‑1 (MQ=255)
                         tCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTg                                                                                                                           >  6:129849/1‑148 (MQ=255)
                             cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAgg                                                                                                                                       >  5:43950/1‑132 (MQ=255)
                             cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAgg                                                                                                                                       <  6:43950/132‑1 (MQ=255)
                             cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                                      <  6:40306/149‑1 (MQ=255)
                             cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                                      <  1:223808/149‑1 (MQ=255)
                             cTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATgg                                                                                                                      <  6:252005/149‑1 (MQ=255)
                                 aTGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTACCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCa                                                                                                                  >  5:62210/1‑149 (MQ=255)
                                       gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGc                                                                                                             >  2:189237/1‑148 (MQ=255)
                                       gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                                                >  5:367316/1‑113 (MQ=255)
                                       gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                                                >  5:362078/1‑113 (MQ=255)
                                       gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                                                <  6:367316/113‑1 (MQ=255)
                                       gctCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCGGGTCCGCGACCAgcagg                                                                                                                                                <  6:362078/113‑1 (MQ=255)
                                         tCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCAGcgag                                                                                                          >  5:76276/1‑146 (MQ=255)
                                           cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                                        >  7:51874/1‑149 (MQ=255)
                                           cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                                        <  4:134900/149‑1 (MQ=255)
                                           cGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGc                                                                                                        <  2:267547/149‑1 (MQ=255)
                                                ggCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATg                                                                                                   <  5:52628/149‑1 (MQ=255)
                                                   gccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCt                                                                                                 >  3:276124/1‑148 (MQ=255)
                                                   gccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                                                <  6:171720/149‑1 (MQ=255)
                                                    ccgcGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCa                                                                                               <  4:353808/149‑1 (MQ=255)
                                                        gACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCAc                                                                                           >  5:74253/1‑149 (MQ=255)
                                                             ttGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGt                                                                                      <  8:13366/149‑1 (MQ=255)
                                                                 agcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgcc                                                                                  <  6:320672/149‑1 (MQ=255)
                                                                  gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                                                >  3:60449/1‑134 (MQ=255)
                                                                  gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTc                                                                                                <  4:60449/134‑1 (MQ=255)
                                                                  gcaCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgccg                                                                                 >  4:151736/1‑149 (MQ=255)
                                                                   caCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCACAGCACCGGGTGgccgc                                                                                >  5:305333/1‑149 (MQ=255)
                                                                    aCGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGgccgcc                                                                               <  8:51874/149‑1 (MQ=255)
                                                                      gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                                             >  8:98927/1‑149 (MQ=255)
                                                                      gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                                             <  6:208112/149‑1 (MQ=255)
                                                                      gCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCaa                                                                             <  6:213236/149‑1 (MQ=255)
                                                                        tGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATa                                                                           <  6:74253/149‑1 (MQ=255)
                                                                          gggTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGATGATGCTCAGCACCGGGTGGCCGCCAATACg                                                                         <  3:52133/149‑1 (MQ=255)
                                                                            gTGAATGCGCACGCCAGAGCGTGCCCCCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGCCCGGGCATGGGTCATCCGCCTAGCTGATGCTCAGCACCGGGTGGCCTCCACTACGCt                                                                       >  8:79097/1‑149 (MQ=255)
                                                                            gTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCt                                                                       >  4:234366/1‑149 (MQ=255)
                                                                             tGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTg                                                                      <  3:151736/149‑1 (MQ=255)
                                                                                    gcACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGc                                                               >  7:124427/1‑149 (MQ=255)
                                                                                     cACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGcc                                                              <  3:208747/149‑1 (MQ=255)
                                                                                     cACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGcc                                                              <  2:313043/149‑1 (MQ=255)
                                                                                       cGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGa                                                            <  5:15306/149‑1 (MQ=255)
                                                                                       cGCCAGAGCGTACCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGa                                                            <  6:62210/149‑1 (MQ=255)
                                                                                          cAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcg                                                         >  6:333990/1‑149 (MQ=255)
                                                                                           agagCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGCCCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcgc                                                        >  4:156291/1‑149 (MQ=255)
                                                                                            gagCGTGCCCCCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATcgcg                                                       >  1:163117/1‑149 (MQ=255)
                                                                                               cGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCAGCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTg                                                    <  7:79097/149‑1 (MQ=255)
                                                                                                 tGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGc                                                  >  4:131786/1‑149 (MQ=255)
                                                                                                   ccACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCAACCGTGCATGGGCCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGccc                                                >  4:90506/1‑149 (MQ=255)
                                                                                                    cACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCt                                               >  7:285108/1‑149 (MQ=255)
                                                                                                               gaTTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCTGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTg                                    <  7:314789/149‑1 (MQ=255)
                                                                                                                   gACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggc                                                                                                                                               >  2:236687/1‑38 (MQ=255)
                                                                                                                   gACTTCAAAGGTCACCTGCTGGTCGCGGACCCgcaggc                                                                                                                                               <  1:236687/38‑1 (MQ=255)
                                                                                                                   gACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCaa                                >  8:82612/1‑149 (MQ=255)
                                                                                                                     cTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATg                              >  1:120344/1‑149 (MQ=255)
                                                                                                                      ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCTTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGGCAATg                              >  7:171400/1‑148 (MQ=255)
                                                                                                                      ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATg                              >  5:33470/1‑148 (MQ=255)
                                                                                                                      ttCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGAAGAACAGGCTGCAGACCGTGCATGGGTCATACGCCTCGCTGATGATCAGCACCGTGTGGCCGCCAATACGCTGGAACAGCCGATAGCGCTGGCCCTCGTCAAGCTTGCCAAtc                              >  4:238422/1‑147 (MQ=255)
                                                                                                                       tCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGGCAATg                              >  3:242810/1‑147 (MQ=255)
                                                                                                                        cAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGt                             <  3:238422/147‑1 (MQ=255)
                                                                                                                         aaaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTaga                          <  2:120344/149‑1 (MQ=255)
                                                                                                                         aaaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTa                            >  5:275699/1‑147 (MQ=255)
                                                                                                                          aaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTagag                         <  6:33470/149‑1 (MQ=255)
                                                                                                                          aaGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTagag                         <  3:353608/149‑1 (MQ=255)
                                                                                                                             gTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAggg                       >  7:248779/1‑148 (MQ=255)
                                                                                                                              tCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCTCTGGTCAAGCTTGCCAATGTAGAGGGCa                     >  1:359187/1‑149 (MQ=255)
                                                                                                                              tCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGc                      >  4:81272/1‑148 (MQ=255)
                                                                                                                                aCCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCt                                                                                                       >  7:150774/1‑65 (MQ=255)
                                                                                                                                aCCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCt                                                                                                       <  8:150774/65‑1 (MQ=255)
                                                                                                                                   tgctgGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTgg                                                   <  6:120290/114‑1 (MQ=255)
                                                                                                                                   tgctgGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTgg                                                   >  5:120290/1‑114 (MQ=255)
                                                                                                                                     ctgGTCGCGGACCCGCAGGCAGAACAGGCTTCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGTCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGcc                                  <  3:367777/129‑1 (MQ=255)
                                                                                                                                                 ccGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGCATCGCAGGCCTGGGCGACCt  <  1:189237/149‑1 (MQ=255)
                                                                                                                                                 ccGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGCATCGCAGGCCTGGGCGACCt  >  1:156054/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
GTCGGGCGCACGGTACGCTTGCGGGTCGGCTTCATGGCTGCTCCGCCCGGCGCCGCGACAATTGCAGCACGCTGGGGTGAATGCGCACGCCAGAGCGTGCCACCGAGTCGAGATTGACTTCAAAGGTCACCTGCTGGTCGCGGACCCGCAGGCAGAACAGGCTGCCGACCGTGCATGGGTCATCCGCCTCGCTGATGCTCAGCACCGGGTGGCCGCCAATACGCTGGAACAGCCGATCGCGCTGGCCCTGGTCAAGCTTGCCAATGTAGAGGGCATCGCAGGCCTGGGCGACCT  >  NC_002947/5297328‑5297621

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