Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,691,828 T→G 8.9% intergenic (+53/‑118) PP_1487 → / → PP_1488 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,691,8280TG 8.9% 241.2 / 9.3 87intergenic (+53/‑118)PP_1487/PP_1488hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (43/36);  new base G (4/4);  total (47/40)
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.91e-01

GAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATC  >  NC_002947/1691680‑1691975
                                                                                                                                                    |                                                                                                                                                   
gAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTg                                                                                                                                                    >  1:1476659/1‑150 (MQ=255)
 aGCGTGCCTGCGTGGACGGGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTgc                                                                                                                                                   <  2:484826/150‑1 (MQ=255)
       ccTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCt                                                                                                                                             >  2:1593303/1‑150 (MQ=255)
         tACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTcc                                                                                                                                           <  1:2245771/150‑1 (MQ=255)
          aCCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCg                                                                                                                                          <  2:740797/150‑1 (MQ=255)
             tGGACGAGCTGGCCCGCAAGCTGGTGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGtct                                                                                                                                       <  2:1476659/150‑1 (MQ=255)
                  gAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAg                                                                                                                                  <  1:2449897/150‑1 (MQ=255)
                  gAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAg                                                                                                                                  <  2:1893780/150‑1 (MQ=255)
                     cTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTcc                                                                                                                                           >  2:256416/1‑138 (MQ=255)
                     cTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTcc                                                                                                                                           <  1:256416/138‑1 (MQ=255)
                     cTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTcgat                                                                                                                                         <  1:1403400/140‑4 (MQ=38)
                     cTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTcgat                                                                                                                                         >  2:1403400/1‑137 (MQ=39)
                       ggCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTTGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGgc                                                                                                                                <  1:2213986/147‑1 (MQ=255)
                       ggCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGgc                                                                                                                                >  2:2213986/1‑147 (MQ=255)
                          ccGCCAGCTGCGGCTGGACCCGTGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTTGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCaaaa                                                                                                                          >  2:562900/1‑150 (MQ=255)
                          ccGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCaaaa                                                                                                                          >  2:1659960/1‑150 (MQ=255)
                          ccGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGAGCCTGTGAGATCGAGCGCCGCCCGTGCGACGCTCCGTCTCACAGCCGCaaaa                                                                                                                          >  2:86192/1‑150 (MQ=255)
                           cGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAAt                                                                                                                         >  1:2043908/1‑150 (MQ=255)
                              cAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTcc                                                                                                                      <  2:1737493/150‑1 (MQ=255)
                               aGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCt                                                                                                                     <  1:1826044/150‑1 (MQ=255)
                                    cGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTgacga                                                                                                                <  1:2314521/150‑1 (MQ=255)
                                       cTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGCGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACaa                                                                                                             <  1:2588685/150‑1 (MQ=255)
                                       cTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACaa                                                                                                             >  1:1427463/1‑150 (MQ=255)
                                           aCCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATcac                                                                                                         >  1:49631/1‑150 (MQ=255)
                                             ccGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACAt                                                                                                       <  1:2630781/150‑1 (MQ=255)
                                              cGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACTGCCGATTTAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATc                                                                                                      <  1:1593303/150‑1 (MQ=255)
                                              cGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATc                                                                                                      <  2:49810/150‑1 (MQ=255)
                                                        ccAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTaa                                                                                            >  2:2796827/1‑150 (MQ=255)
                                                        ccAGCAACTGGAAAGCCAGGTGCGTTTGGCTGCAGCGTAGTTGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGACATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTaa                                                                                            <  1:1484388/150‑1 (MQ=255)
                                                            gAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGAcc                                                                                        >  2:2180564/1‑150 (MQ=255)
                                                             aaCTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACct                                                                                       >  1:3039879/1‑150 (MQ=255)
                                                               cTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACGTCAATa                                                                                                  >  2:2021708/1‑137 (MQ=255)
                                                               cTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATa                                                                                                  <  1:2021708/137‑1 (MQ=255)
                                                               cTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACctct                                                                                     >  1:2690266/1‑150 (MQ=255)
                                                                tGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTg                                                                                    >  2:2878495/1‑150 (MQ=255)
                                                                   aaaGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAAt                                                                                                            <  2:409136/123‑1 (MQ=255)
                                                                   aaaGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAAt                                                                                                            >  1:409136/1‑123 (MQ=255)
                                                                       ccAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGtct                                                                                                                                       <  2:2935235/92‑1 (MQ=255)
                                                                       ccAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGtct                                                                                                                                       >  1:2935235/1‑92 (MQ=255)
                                                                              cGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGc                                                                      <  2:1936360/150‑1 (MQ=255)
                                                                                gTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGtct                                                                                                                                       >  2:28484/1‑83 (MQ=255)
                                                                                gTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGtct                                                                                                                                       <  1:28484/83‑1 (MQ=255)
                                                                                  tGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGAtt                                                                  <  2:3039879/150‑1 (MQ=255)
                                                                                      tGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTc                                                              >  1:608077/1‑150 (MQ=255)
                                                                                         aGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCg                                                           <  1:2180564/150‑1 (MQ=255)
                                                                                         aGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTcgat                                                                                                                                         >  1:72120/1‑69 (MQ=17)
                                                                                         aGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTcgat                                                                                                                                         <  2:72120/72‑4 (MQ=17)
                                                                                          gCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATc                                                                                                      <  2:2077994/106‑1 (MQ=255)
                                                                                          gCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATc                                                                                                      >  1:2077994/1‑106 (MQ=255)
                                                                                               gTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGtt                                                     <  1:2796827/150‑1 (MQ=255)
                                                                                                    ggTGCTTGCCACAGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGaa                                                >  1:2585261/1‑150 (MQ=255)
                                                                                                     gTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCg                                                                        <  2:1231906/125‑1 (MQ=255)
                                                                                                     gTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCg                                                                        >  1:1231906/1‑125 (MQ=255)
                                                                                                     gTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAAt                                               >  2:609808/1‑150 (MQ=255)
                                                                                                      tGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTaa                                                                                            <  2:1369173/104‑1 (MQ=255)
                                                                                                      tGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTaa                                                                                            >  1:1369173/1‑104 (MQ=255)
                                                                                                       gCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTc                                             >  1:1796771/1‑150 (MQ=255)
                                                                                                        cTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTcc                                            >  1:1524377/1‑150 (MQ=255)
                                                                                                         ttGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCt                                           >  2:1912745/1‑150 (MQ=255)
                                                                                                           gCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGa                                         <  2:1796771/150‑1 (MQ=255)
                                                                                                             cACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAgg                                       >  1:2944135/1‑150 (MQ=255)
                                                                                                                   cTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGAttttt                                                               >  2:2687290/1‑120 (MQ=255)
                                                                                                                   cTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGaa                                 >  2:2928078/1‑150 (MQ=255)
                                                                                                                   cTATTGAGCGCCTGTGAGACCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGAttttt                                                               <  1:2687290/120‑1 (MQ=255)
                                                                                                                    tATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAt                                <  2:2585261/150‑1 (MQ=255)
                                                                                                                    tATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAt                                >  2:992669/1‑150 (MQ=255)
                                                                                                                              cTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACt                      >  2:2071705/1‑150 (MQ=255)
                                                                                                                              cTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACt                      <  2:49631/150‑1 (MQ=255)
                                                                                                                              cTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCTAAATTCCTGACGACAATCACATCAATACCGTaa                                                                                            >  1:1677245/1‑80 (MQ=14)
                                                                                                                              cTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCTAAATTCCTGACGACAATCACATCAATACCGTaa                                                                                            <  2:1677245/80‑1 (MQ=14)
                                                                                                                              cTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCTAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAgg                                       <  2:2129387/133‑1 (MQ=37)
                                                                                                                              cTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAGGCGCTAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAgg                                       >  1:2129387/1‑133 (MQ=25)
                                                                                                                                 tgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGTGCAAAATTCCTGACGACAATCACATCAATACCGTAACACCTCTCGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTCATTTCCTGAGGGCTGAATGTGAAGAACTGGa                   <  1:562900/150‑1 (MQ=255)
                                                                                                                                 tgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAt                                >  2:1686419/1‑137 (MQ=255)
                                                                                                                                 tgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAt                                <  1:1686419/137‑1 (MQ=255)
                                                                                                                                 tgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGa                   >  2:1906808/1‑150 (MQ=255)
                                                                                                                                   agaTCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGAcc                 >  2:2888496/1‑150 (MQ=255)
                                                                                                                                     aTCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTtgt                                                   >  1:52270/1‑114 (MQ=255)
                                                                                                                                     aTCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTtgt                                                   <  2:52270/114‑1 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGTGTGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACGGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTgcg            >  2:2205052/1‑150 (MQ=255)
                                                                                                                                        gAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACATTgag            >  1:256845/1‑148 (MQ=255)
                                                                                                                                         aGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTgcgc           >  2:1693871/1‑150 (MQ=255)
                                                                                                                                            gccgccCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACGTTGCGCCaa        >  1:1125611/1‑150 (MQ=255)
                                                                                                                                             ccgccCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAAc       <  1:2888496/150‑1 (MQ=255)
                                                                                                                                                 ccGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGAt   <  1:1510079/150‑1 (MQ=255)
                                                                                                                                                  cGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGCGGTGTGAATGTGAAGAACTGGACCTTGCGCCAACGGATc  >  1:940651/1‑150 (MQ=255)
                                                                                                                                                  cGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATc  >  1:459275/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGGCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATC  >  NC_002947/1691680‑1691975

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: