Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,642,657 A→G 7.4% 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,642,6570AG 7.4% 198.4 / 4.1 67Q185R (CAG→CGG) ‡PP_4107LysR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base G (2/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

ATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTCCG  >  NC_002947/4642508‑4642806
                                                                                                                                                      |                                                                                                                                                     
aTTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                       >  1:1178574/1‑150 (MQ=255)
      ccTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCGCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGct                                                                                                                                                 >  1:549237/1‑150 (MQ=255)
         ggTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgta                                                                                                                                              <  2:135465/150‑3 (MQ=255)
             cAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgtcgaaa                                                                                                                                          <  1:1014100/150‑7 (MQ=255)
             cAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAg                                                                                                                                          <  2:858861/150‑1 (MQ=255)
                   aaGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGc                                                                                                                                    >  1:18972/1‑150 (MQ=255)
                        ggACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                               <  2:738145/150‑1 (MQ=255)
                            gTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCaa                                                                                                                           >  1:1034499/1‑150 (MQ=255)
                            gTCGCCTTTTTCTACCGCAAGCGAACCTGGCCCCGGTCGACCTGCATCGAGTTGTTCAGAGAAGTAGGACGATGTGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCaa                                                                                                                           <  1:542639/150‑1 (MQ=255)
                               gCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctga                                                                                                                                               <  2:954355/127‑1 (MQ=255)
                               gCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGctga                                                                                                                                               >  1:954355/1‑127 (MQ=255)
                                   ttttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCgtgt                                                                                                                    <  2:142986/150‑1 (MQ=255)
                                    tttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCaa                                                                                                                           >  2:162236/1‑142 (MQ=255)
                                    tttCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCaa                                                                                                                           <  1:162236/142‑1 (MQ=255)
                                               aGGGCACCTGGCCCGGGGCGGCCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAgg                                                                                                        <  1:1261854/150‑1 (MQ=255)
                                                gggCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGc                                                                                                       >  1:404421/1‑150 (MQ=255)
                                                    aCCTGGCCCGGGGCGATCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTgg                                                                                                   >  1:70242/1‑150 (MQ=255)
                                                          cccGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGaa                                                                                             >  1:592512/1‑150 (MQ=255)
                                                                    ccTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTgg                                                                                   <  2:18972/150‑1 (MQ=255)
                                                                     cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCt                                                                                     <  2:1054841/147‑1 (MQ=255)
                                                                     cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCt                                                                                     >  1:1054841/1‑147 (MQ=255)
                                                                     cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                  >  1:217279/1‑150 (MQ=255)
                                                                     cTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                  >  2:793355/1‑150 (MQ=255)
                                                                     cTGCATCGAGTGGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTggg                                                                                  >  1:142918/1‑150 (MQ=255)
                                                                            gAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                       >  1:1097603/1‑74 (MQ=255)
                                                                            gAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCa                                                                                                                                                       <  2:1097603/74‑1 (MQ=255)
                                                                              gttgttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGtt                                                                         <  1:793355/150‑1 (MQ=255)
                                                                                 gttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCa                                                                      <  2:1034499/150‑1 (MQ=255)
                                                                                  ttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgtc                                                                                                                                              >  2:1239576/1‑75 (MQ=255)
                                                                                  ttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCTGGGCGctgtc                                                                                                                                              <  1:1239576/77‑3 (MQ=255)
                                                                                  ttCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCAt                                                                     >  2:341894/1‑150 (MQ=255)
                                                                                   tCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                    <  1:1280478/150‑1 (MQ=255)
                                                                                     aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGc                                                                                                              >  2:438675/1‑106 (MQ=255)
                                                                                     aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGc                                                                                                              <  1:438675/106‑1 (MQ=255)
                                                                                     aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                  <  2:70242/150‑1 (MQ=255)
                                                                                     aGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                  >  2:747554/1‑150 (MQ=255)
                                                                                         aagaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCCCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGa                                                                                              <  1:288019/118‑1 (MQ=255)
                                                                                         aagaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGa                                                                                              >  2:288019/1‑118 (MQ=255)
                                                                                          agaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGaa                                                             <  2:404421/150‑1 (MQ=255)
                                                                                           gaagTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGGCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAAc                                                            >  1:1103831/1‑150 (MQ=255)
                                                                                              gtggtgCCGGTGTGCAGCCCGCAGCTGCCGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGc                                                         >  1:482484/1‑150 (MQ=255)
                                                                                                  tgCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                    <  2:604386/135‑1 (MQ=255)
                                                                                                  tgCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATa                                                                    >  1:604386/1‑135 (MQ=255)
                                                                                                     cGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACg                                                  >  2:1193406/1‑150 (MQ=255)
                                                                                                      ggTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACgg                                                 <  2:1178574/150‑1 (MQ=255)
                                                                                                      ggTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACgg                                                 >  2:1203304/1‑150 (MQ=255)
                                                                                                      ggTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACgg                                                 >  1:673530/1‑150 (MQ=255)
                                                                                                        gtgtGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGAGAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACggg                                                >  2:738496/1‑150 (MQ=255)
                                                                                                              aGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCgcgc                                          <  2:549237/150‑1 (MQ=255)
                                                                                                                cccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCCGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCtt                                        >  2:1150769/1‑150 (MQ=255)
                                                                                                                 ccGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCACAACAGCTACCACGGGCCGCGCTTc                                       <  1:747554/150‑1 (MQ=255)
                                                                                                                  cGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCtt                                        <  1:62594/148‑1 (MQ=255)
                                                                                                                  cGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCtt                                        >  2:62594/1‑148 (MQ=255)
                                                                                                                           tgGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCt                             <  2:243550/150‑1 (MQ=255)
                                                                                                                               cagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                  >  1:229716/1‑109 (MQ=255)
                                                                                                                               cagccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGc                                                                  <  2:229716/109‑1 (MQ=255)
                                                                                                                               cagccaCCGTTGCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGTTTGCATAGCCAGAACAGCTACCACGGGCCGCGCCTCCACACGTTCTACCt                         <  2:142918/150‑1 (MQ=255)
                                                                                                                                 gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                       <  1:1003055/150‑1 (MQ=255)
                                                                                                                                 gccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                       <  2:482484/150‑1 (MQ=255)
                                                                                                                                  ccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTggg                                                                           >  2:114665/1‑97 (MQ=255)
                                                                                                                                  ccaCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTggg                                                                           <  1:114665/97‑1 (MQ=255)
                                                                                                                                   caCCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATGGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTgtgt                     <  1:1201918/150‑1 (MQ=255)
                                                                                                                                      cGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATc                  <  1:1193406/150‑1 (MQ=255)
                                                                                                                                        tttCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCcg                >  1:995720/1‑150 (MQ=255)
                                                                                                                                            gACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                               <  2:1031161/35‑1 (MQ=255)
                                                                                                                                            gACAGCGCCCAGGCGCTGACTGAGCACCGCCTGtt                                                                                                                               >  1:1031161/1‑35 (MQ=255)
                                                                                                                                             tcAGCGCCTGGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCAGGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAAGAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCgg           >  1:683595/2‑150 (MQ=255)
                                                                                                                                             aCAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                       <  2:57114/138‑1 (MQ=255)
                                                                                                                                             aCAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCtgt                       >  1:57114/1‑138 (MQ=255)
                                                                                                                                                    ccAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCAtctc    <  1:1150769/150‑1 (MQ=255)
                                                                                                                                                      aGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTCCg  >  2:947511/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
ATTCGACCTGGTACAGGCCAAGGCGGACGTCGCCTTTTTCTTCGGCCAGGGCACCTGGCCCGGGGCGACCTGCATCGAGTTGTTCAGCGAAGAAGTGGTGCCGGTGTGCAGCCCGCAGCTGCTGGCCAGCCACCGTTTCGACAGCGCCCAGGCGCTGACTGAGCACCGCCTGTTGCAATGCGTGTCGCGGCCGGAGGCCTGGCATGAATGGTTCCTGGGGCTTGGGTTGCATAGCCAGAACAGCTACCACGGGCCGCGCTTCGACACGTTCTACCTGTGTATCCGCGCGGCCATCTCCG  >  NC_002947/4642508‑4642806

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