Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,338,538 A→T 12.4% intergenic (+114/+28) PP_3812 → / ← PP_5632 lipase/allo‑isoleucine metabolism protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,338,5380AT12.4% 233.5 / 22.0 88intergenic (+114/+28)PP_3812/PP_5632lipase/allo‑isoleucine metabolism protein
Reads supporting (aligned to +/- strand):  ref base A (40/37);  new base T (7/4);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.49e-01

GGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAA  >  NC_002947/4338391‑4338685
                                                                                                                                                   |                                                                                                                                                   
ggCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTTGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAc                                                                                                                                                   <  2:1052808/150‑1 (MQ=255)
  cGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAccc                                                                                                                                                 <  2:1471397/150‑1 (MQ=255)
    ggCAACGAACTTAAAGTGCACCTCCAGTGAGGCGACGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTAGTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCcgc                                                                                                                                               <  1:1361698/150‑1 (MQ=255)
         cGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAg                                                                                                                                          >  2:110780/1‑150 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                         <  2:2398435/150‑1 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                                                                                                         <  1:158024/150‑1 (MQ=255)
          gAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAgg                                                                                                                                         >  2:1136966/1‑150 (MQ=255)
                  aGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCaa                                                                                                                                 >  1:1831859/1‑150 (MQ=255)
                    tGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGCAAGGGCCGGCAAAt                                                                                                                               >  2:966259/1‑150 (MQ=255)
                          tGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCt                                                                                                                         >  1:1403590/1‑150 (MQ=255)
                            cAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGcc                                                                                                                                      <  2:191496/135‑1 (MQ=255)
                            cAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTTTACCCGCGAAAGGGcc                                                                                                                                      >  1:191496/1‑135 (MQ=255)
                            cAGTGAGGCGATGAGCAATAGACACGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAAAAGGCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTACGGAAAGTTCGCGGGTAAGCCTGTGAAAGGGCGGGGAAATCGTGCtaa                                                                                                                       >  2:1734721/1‑148 (MQ=255)
                             aGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAg                                                                                                                      >  2:929681/1‑150 (MQ=255)
                              gTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAgg                                                                                                                     <  2:1403590/150‑1 (MQ=255)
                                 aGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTGTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGACAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCATGGCt                                                                                                                  >  1:466010/1‑150 (MQ=255)
                                 aGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAccc                                                                                                                                                 <  1:1840390/119‑1 (MQ=255)
                                 aGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAAccc                                                                                                                                                 >  2:1840390/1‑119 (MQ=255)
                                        aaGCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAAcac                                                                                                           >  2:2304938/2‑148 (MQ=255)
                                          gCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCtt                                                                                                         >  1:1421373/1‑150 (MQ=255)
                                          gCAAAAGACCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCGtt                                                                                                         >  1:1610483/1‑150 (MQ=255)
                                           caaacgccccGCCCAGGCGGGTTTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTc                                                                                                        <  1:1474827/143‑1 (MQ=255)
                                            aaaaGACCCGCCCAGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                       >  1:2275750/4‑150 (MQ=255)
                                            aaTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                       >  1:154193/1‑150 (MQ=255)
                                             aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca                                                                                                      <  2:1831859/150‑1 (MQ=255)
                                             aTAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCca                                                                                                      >  1:2046980/1‑150 (MQ=255)
                                              tAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCg                                                                                                                                     >  1:1036227/1‑118 (MQ=255)
                                              tAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCg                                                                                                                                     <  2:1036227/118‑1 (MQ=255)
                                                gACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACGCGCGTAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACaa                                                                                                   >  1:2070339/1‑150 (MQ=255)
                                                 aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCCGCTCAAGCCCTTTCACGGGTAAACCCGCGAAAGGGCCGGCAAATCGTGCTCCGGGCTGAAACCCTTCCACAAt                                                                                                  >  1:42470/1‑150 (MQ=255)
                                                 aCCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                  <  2:773130/150‑1 (MQ=255)
                                                 aCCCGCCCAGGCGGGTCTATTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCGCAGGGCTCAAACCCTTCCACAAt                                                                                                  >  1:613203/1‑150 (MQ=255)
                                                  cccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATg                                                                                                 <  1:959886/150‑1 (MQ=255)
                                                  cccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATg                                                                                                 <  1:1598548/150‑1 (MQ=255)
                                                   ccGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGGAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGa                                                                                                >  2:1316840/1‑150 (MQ=255)
                                                       ctcggCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGGCAGCGGTGCGGCAGGTCCGGCTCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCAc                                                                                            <  1:1677353/147‑1 (MQ=255)
                                                             gggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTg                                                                                                                           <  1:867261/113‑1 (MQ=255)
                                                             gggTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTg                                                                                                                           >  2:867261/1‑113 (MQ=255)
                                                                       cTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAcc                                                                            >  1:1304547/1‑150 (MQ=255)
                                                                        ttCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAccc                                                                           >  1:3994/1‑150 (MQ=255)
                                                                         tCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCAcccc                                                                          >  1:724855/1‑150 (MQ=255)
                                                                             tttGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTgg                                                                      <  2:1421373/150‑1 (MQ=255)
                                                                             tttGCGGAAACAGACCGCACGTTTCAATGAATCTGCCAGCGGTGCGGGAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTgg                                                                      <  2:1158843/150‑1 (MQ=255)
                                                                                 cGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGAt                                                                  <  2:2180070/150‑1 (MQ=255)
                                                                                    aaaCAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGc                                                               >  2:980371/1‑150 (MQ=255)
                                                                                             gCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCct                                                      <  2:434383/150‑1 (MQ=255)
                                                                                                 gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCGGGCTGCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                  <  2:154193/150‑1 (MQ=255)
                                                                                                 gCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTg                                                  <  1:989649/150‑1 (MQ=255)
                                                                                                  cTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTgg                                                 <  2:2046980/150‑1 (MQ=255)
                                                                                                   ttCAATGAATCTGCCAGCGGTTCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGt                                                >  1:2000836/1‑150 (MQ=255)
                                                                                                         gAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAgcg                                          >  2:922453/1‑150 (MQ=255)
                                                                                                           atcggtcAGCGGTGCGGCAGGTCCGGCCCTTTCGTGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGcgcg                                        <  2:613203/144‑1 (MQ=255)
                                                                                                           aTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATgg                                                                >  2:2150486/1‑126 (MQ=255)
                                                                                                           aTCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATgg                                                                <  1:2150486/126‑1 (MQ=255)
                                                                                                            tCTGCAAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTACTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCgg                                       >  2:1730824/1‑150 (MQ=255)
                                                                                                              tGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATCATCACCTCAGCCTTCGCCACCCCtt                                                                        <  2:1837306/115‑1 (MQ=255)
                                                                                                              tGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATACTCACCTCAGCCTTCGCCACCCCtt                                                                        >  1:1837306/1‑115 (MQ=255)
                                                                                                                ccAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCg                                   <  1:476754/150‑1 (MQ=255)
                                                                                                                 cAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCgg                                  <  2:2275750/150‑1 (MQ=255)
                                                                                                                  aGCGGGTCGTCAGGTCCGGCCCTTGTGCTGGTAAACCCGCGAAAGGGCCGGCACATCCTGCTCAGGGCTCAAACCCTTCCACACTGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGt                                 <  1:2300158/150‑1 (MQ=255)
                                                                                                                   gCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGCTa                                <  2:1844128/150‑1 (MQ=255)
                                                                                                                   gCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCAGAGCGGTa                                >  1:1319208/1‑150 (MQ=255)
                                                                                                                       tGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAg                            <  1:1316840/150‑1 (MQ=255)
                                                                                                                        gCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAgg                           <  2:1320987/150‑1 (MQ=255)
                                                                                                                           gccggTCCGGCCCTTTCGCGGGTAAACCTGCGAAAGGGCCGGCAAATCCTCCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGa                        <  2:42470/147‑1 (MQ=255)
                                                                                                                           gCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGa                        >  1:2304300/1‑150 (MQ=255)
                                                                                                                             aGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACg                      <  1:922453/150‑1 (MQ=255)
                                                                                                                               gTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTaa                               <  1:2052332/139‑1 (MQ=255)
                                                                                                                               gTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTaa                               >  2:2052332/1‑139 (MQ=255)
                                                                                                                                tCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGCAACAGGCGACGGcc                   >  1:2388849/1‑150 (MQ=255)
                                                                                                                                      ccTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                  >  1:1428346/1‑65 (MQ=255)
                                                                                                                                      ccTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAAt                                                                                                  <  2:1428346/65‑1 (MQ=25)
                                                                                                                                       cgttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                       <  2:1667903/57‑1 (MQ=255)
                                                                                                                                       cTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAg                                     >  1:900057/1‑125 (MQ=255)
                                                                                                                                       cTTTCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAg                                     <  2:900057/125‑1 (MQ=255)
                                                                                                                                       cTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTcc                                                                                                       >  1:1667903/1‑59 (MQ=255)
                                                                                                                                       cTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGt            >  2:1151093/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  1:1908205/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  2:1156108/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATCCTACTCAGGGCTCAAACCCTTCCACAATGATCACCACAGCCTTCGCCCCCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  1:1976689/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTTTACCCGCGAAAGGGCCGGCAAATACTGCTCAGGGCTCAAACCCTTCCACAATGATCACTTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  1:1976523/1‑150 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGttag                              <  1:1892639/131‑4 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAAc                              >  2:1892639/1‑131 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           <  2:2304300/150‑1 (MQ=255)
                                                                                                                                        tttCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTa           >  2:814634/1‑150 (MQ=255)
                                                                                                                                            gcgGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGaa       >  1:121028/1‑150 (MQ=255)
                                                                                                                                               ggTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGGATCg    >  1:200561/1‑150 (MQ=255)
                                                                                                                                                 tAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGaa  <  1:929681/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGCGGGCAACGAACTTAAAGTGCATCTGCAGTGAGGCGATGAGCAATAGACCCGCCTGGGCGGGTCTTTTGCTTCTATTTGCGGAAACAGCCCGCAGGCTTCAATGAATCTGCCAGCGGTGCGGCAGGTCCGGCCCTTTCGCGGGTAAACCCGCGAAAGGGCCGGCAAATCCTGCTCAGGGCTCAAACCCTTCCACAATGATCACCTCAGCCTTCGCCACCCCTTGGCGATGGCTGCATGCCTCTTGGTACAGCGCGGAGCGGTAACAGGCGACGGCCTGCGCGTACGAATCGAA  >  NC_002947/4338391‑4338685

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