Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,647,606 A→G 5.7% Q185R (CAG→CGG) ‡ PP_4107 → LysR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,647,6060AG 5.7% 1323.3 / 40.8 401Q185R (CAG→CGG) ‡PP_4107LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (180/198);  new base G (14/9);  total (194/207)
Fisher's exact test for biased strand distribution p-value = 2.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 410 total aligned reads displayed.

ATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTCCG  >  NC_002947/4647457‑4647755
                                                                                                                                                     |                                                                                                                                                     
atacGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGACGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                       >  2:1758328/4‑150 (MQ=255)
aTTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                       >  1:6055431/1‑150 (MQ=255)
 ttCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAg                                                                                                                                                      >  2:4179867/1‑150 (MQ=255)
 ttCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAg                                                                                                                                                      >  1:4128478/1‑150 (MQ=255)
     aCCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGgcgc                                                                                                                                                  >  1:2282748/1‑150 (MQ=255)
        tGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgt                                                                                                                                               <  2:5985343/150‑2 (MQ=255)
        tGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgt                                                                                                                                               <  2:6302990/150‑2 (MQ=255)
        tGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctga                                                                                                                                               >  2:1865223/1‑150 (MQ=255)
         ggTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCTctgtc                                                                                                                                              <  2:6355784/150‑3 (MQ=255)
         ggTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgtc                                                                                                                                              <  2:665077/150‑3 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGTCAGCGCCTGGGCGctgtcg                                                                                                                                             >  1:5371411/1‑147 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgact                                                                                                                                             <  1:3987739/150‑1 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgact                                                                                                                                             <  2:7221969/150‑1 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgact                                                                                                                                             >  1:2486504/1‑150 (MQ=255)
          gTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgact                                                                                                                                             <  2:5891428/150‑1 (MQ=255)
            aCAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgactga                                                                                                                                           <  1:5739061/150‑1 (MQ=255)
            aCAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgactga                                                                                                                                           <  1:2582951/150‑1 (MQ=255)
            aCAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctgactga                                                                                                                                           <  2:7292706/150‑1 (MQ=255)
             cAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGTCAGCGCCTGGGCGctgtcgaaa                                                                                                                                          <  2:5371411/150‑7 (MQ=255)
              aGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGc                                                                                                                                         >  1:7266282/1‑150 (MQ=255)
                gCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCAc                                                                                                                                       >  1:7201198/1‑150 (MQ=255)
                 ccAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCAcc                                                                                                                                      >  1:3380454/1‑150 (MQ=255)
                      gCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTg                                                                                                                                 <  2:911118/150‑1 (MQ=255)
                        ggACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                               <  2:803708/150‑1 (MQ=255)
                               gCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCTGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGc                                                                                                                        >  1:1464590/1‑150 (MQ=255)
                                 cTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgt                                                                                                                      <  1:1723137/150‑1 (MQ=255)
                                  tttttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgtg                                                                                                                     >  2:1124583/1‑150 (MQ=255)
                                  tttttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgtg                                                                                                                     >  1:356262/1‑150 (MQ=255)
                                   ttttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgtgt                                                                                                                    <  1:7128152/150‑1 (MQ=255)
                                      tcttcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTcgc                                                                                                                 <  2:1305056/150‑1 (MQ=255)
                                       cccgggCCAGGTGCCCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTcgcg                                                                                                                >  2:5147348/5‑150 (MQ=255)
                                        ttcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAg                                                                                                                                          <  2:4101395/123‑1 (MQ=255)
                                        ttcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCgg                                                                                                               <  1:2989135/150‑1 (MQ=255)
                                        ttcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCgg                                                                                                               <  2:3130282/150‑1 (MQ=255)
                                         tcGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGc                                                                                                              <  1:7418999/150‑1 (MQ=255)
                                            gCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCCCCGCCTGTTGCAATGCGTGTCGCGGCCgg                                                                                                           >  2:467768/1‑150 (MQ=255)
                                            gCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCgg                                                                                                           >  2:5831092/1‑150 (MQ=255)
                                             ccAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGa                                                                                                          <  2:5031521/150‑1 (MQ=255)
                                                ggtgcccTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGc                                                                                                       >  1:6753095/6‑150 (MQ=255)
                                                gggCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGc                                                                                                       <  2:6496206/150‑1 (MQ=255)
                                                 ggCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGcc                                                                                                      >  1:3718317/1‑150 (MQ=255)
                                                   cACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTg                                                                                                    <  2:2758929/150‑1 (MQ=255)
                                                    aCCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTgg                                                                                                   <  2:2651863/150‑1 (MQ=255)
                                                     ccTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGc                                                                                                  <  1:6284191/150‑1 (MQ=255)
                                                     ccTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGc                                                                                                  <  2:3731709/150‑1 (MQ=255)
                                                      cTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCa                                                                                                 >  1:5665096/1‑150 (MQ=255)
                                                         gCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGa                                                                                              <  1:1800776/150‑1 (MQ=255)
                                                         gCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGa                                                                                              <  2:2993871/150‑1 (MQ=255)
                                                          cccGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGaa                                                                                             <  2:2259245/150‑1 (MQ=255)
                                                          cccGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGaa                                                                                             <  2:3867440/150‑1 (MQ=255)
                                                                     cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                  >  1:4676538/1‑150 (MQ=255)
                                                                     cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                  >  2:1443055/1‑150 (MQ=255)
                                                                           cGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTgg                                                                            <  1:5048237/150‑1 (MQ=255)
                                                                                 gttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCa                                                                      <  2:4016504/150‑1 (MQ=255)
                                                                                   tCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                    <  1:7053798/150‑1 (MQ=255)
                                                                                    cAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAg                                                                   <  1:4534233/150‑1 (MQ=255)
                                                                                     aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                  >  1:4013512/1‑150 (MQ=255)
                                                                                     aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                  >  2:2424367/1‑150 (MQ=255)
                                                                                        gaagaaGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTg                                                                                                    <  1:795944/113‑1 (MQ=255)
                                                                                        gaagaaGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAg                                                               <  1:99825/150‑1 (MQ=255)
                                                                                          agaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCACCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGaa                                                             <  2:6832090/150‑1 (MQ=255)
                                                                                           gaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGtt                                                                         <  1:6718617/137‑1 (MQ=255)
                                                                                            aagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACa                                                           <  1:1155109/150‑1 (MQ=255)
                                                                                            aagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACa                                                           <  1:1370199/150‑1 (MQ=255)
                                                                                            aagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACa                                                           <  1:2342762/150‑1 (MQ=255)
                                                                                                  tgCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTg                                                                                    >  2:3306430/1‑119 (MQ=255)
                                                                                                  tgCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTg                                                                                    <  1:3306430/119‑1 (MQ=255)
                                                                                                   gCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCa                                                    >  1:4466814/1‑150 (MQ=255)
                                                                                                        tggctggCCAGCAGCTGCGGGCTGCACACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGc                                               >  1:621953/7‑150 (MQ=255)
                                                                                                            cAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCcgcg                                           <  2:1545719/150‑1 (MQ=255)
                                                                                                            cAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCcgcg                                           >  1:7341832/1‑150 (MQ=255)
                                                                                                             aGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGa                                                                                              <  1:3210475/98‑1 (MQ=255)
                                                                                                             aGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCgcgc                                          >  2:4428633/1‑150 (MQ=255)
                                                                                                              gCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCt                                         >  1:4877438/1‑150 (MQ=255)
                                                                                                              gCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCt                                         >  1:3440775/1‑150 (MQ=255)
                                                                                                                ccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTc                                       <  2:7248829/150‑1 (MQ=255)
                                                                                                                ccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTc                                       <  1:3536873/150‑1 (MQ=255)
                                                                                                                  gCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGa                                     >  1:2157505/1‑150 (MQ=255)
                                                                                                                     gctgctGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGacac                                  <  1:5876381/150‑1 (MQ=255)
                                                                                                                       tgctgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGt                                <  2:551630/150‑1 (MQ=255)
                                                                                                                        gctgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGtt                               <  2:2032666/150‑1 (MQ=255)
                                                                                                                         ctgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTc                              <  1:481370/150‑1 (MQ=255)
                                                                                                                         ctgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTc                              >  1:4979298/1‑150 (MQ=255)
                                                                                                                            gccagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTAc                           <  2:5568885/150‑1 (MQ=255)
                                                                                                                            gccagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTAc                           >  1:1022891/1‑150 (MQ=255)
                                                                                                                            gccagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTAc                           <  2:7126625/150‑1 (MQ=255)
                                                                                                                              cagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCt                         >  2:5241431/1‑150 (MQ=255)
                                                                                                                               agccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtg                        >  2:3518935/1‑150 (MQ=255)
                                                                                                                               agccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtg                        <  2:4575784/150‑1 (MQ=255)
                                                                                                                                gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGGCGCGCTTCGACACGTTCTACCtgc                       <  1:2895278/150‑2 (MQ=255)
                                                                                                                                gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                       <  1:4188558/150‑1 (MQ=255)
                                                                                                                                gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                       <  1:1965025/150‑1 (MQ=255)
                                                                                                                                  caCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTgtgt                     >  2:6012950/1‑150 (MQ=255)
                                                                                                                                      gTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATcc                 >  1:6596301/1‑150 (MQ=255)
                                                                                                                                          cGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcgcgc             >  2:7362561/1‑150 (MQ=255)
                                                                                                                                          agtcAGCGCCTGGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCTTAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcgcgc             >  2:6955952/4‑150 (MQ=255)
                                                                                                                                           gtcAGCGCCTGGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcgcgcg            >  2:4805658/3‑150 (MQ=255)
                                                                                                                                            aCAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCgg           >  2:1613109/1‑150 (MQ=255)
                                                                                                                                               gcgcCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCa        >  1:5233629/1‑150 (MQ=255)
                                                                                                                                                   ccAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGc                                                                       >  2:358234/1‑83 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTCCG  >  NC_002947/4647457‑4647755

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: