Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,955,263 A→G 8.1% intergenic (+50/‑143) PP_1752 → / → PP_1753 hypothetical protein/DMT superfamily permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,955,2630AG 8.1% 318.6 / 12.2 111intergenic (+50/‑143)PP_1752/PP_1753hypothetical protein/DMT superfamily permease
Reads supporting (aligned to +/- strand):  ref base A (58/44);  new base G (6/3);  total (64/47)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

Only 100 of 111 total aligned reads displayed.

CTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTATGC  >  NC_002947/1955118‑1955409
                                                                                                                                                 |                                                                                                                                                  
cTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGAACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGt                                                                                                                                                <  1:483121/150‑1 (MQ=255)
      gAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGcac                                                                                                                                          >  1:1825879/1‑150 (MQ=255)
        ggCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAg                                                                                                                                        <  1:415706/150‑1 (MQ=255)
        ggCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAg                                                                                                                                        <  1:848513/150‑1 (MQ=255)
        ggCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAg                                                                                                                                        <  1:191381/150‑1 (MQ=255)
          ccAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGc                                                                                                                                      <  1:629354/150‑1 (MQ=255)
          ccAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCAGCGAAAGGGTCAGCACAGGc                                                                                                                                      <  2:619257/150‑1 (MQ=255)
           cAATGACAAGCGCTGGGCGCGGGTGCAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCa                                                                                                                                     >  1:415131/1‑150 (MQ=255)
            aaTGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCCGCTTTCTCCGCCAAATGGACATCACAGGCaa                                                                                                                                    >  2:504024/1‑150 (MQ=255)
                  aaGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTc                                                                                                                              <  2:1825879/150‑1 (MQ=255)
                            cgcgGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCa                                                                                                                    >  1:573895/1‑150 (MQ=255)
                            cgcgGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCa                                                                                                                    >  1:76160/1‑150 (MQ=255)
                            cgcgGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCa                                                                                                                    <  2:76160/150‑1 (MQ=255)
                               gggTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCa                                                                                                                            >  1:1257685/1‑139 (MQ=255)
                                 gTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGccc                                                                                                               >  2:1105860/1‑150 (MQ=255)
                                   gAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcccc                                                                                                              <  1:1617331/149‑1 (MQ=255)
                                   gAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcccc                                                                                                              >  2:1617331/1‑149 (MQ=255)
                                      ggtggtGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAgg                                                                                                                                       >  2:1235626/1‑121 (MQ=255)
                                      ggtggtGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAgg                                                                                                                                       <  1:1235626/121‑1 (MQ=255)
                                         ggtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGATTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGaa                                                                                                       <  2:1906741/150‑1 (MQ=255)
                                          gtgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAAt                                                                                                      <  2:1234347/150‑1 (MQ=255)
                                          gtacgcACAATCAAGCAGGCACGGGAGGCGGCGTTTGCCAAGCACAGCGAATAGCCGAAGGTGAGGTTCCCTGTTCTGGCTGCTTCGCGGTTTTACCCGCAAAAGCGTCATAACAGGCAACACCTCCACCAAGCCAGGCCCCAACAGAAt                                                                                                      <  2:387959/147‑1 (MQ=255)
                                           tgCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTGTGGCCGCTTCGCGGGTTTACCCGAAAAAGGGGAAAACTACACAACACCTTCACCCAGCCATGCCCCAAACGAATg                                                                                                     >  1:1103163/1‑150 (MQ=255)
                                            gcgcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGa                                                                                                    <  1:1105860/150‑1 (MQ=255)
                                             cgcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                                          >  2:835394/1‑143 (MQ=255)
                                             cgcACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                                          <  1:835394/143‑1 (MQ=255)
                                               cacaATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCt                                                                                                 <  2:1249806/150‑1 (MQ=255)
                                                  aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                              >  1:1458346/1‑150 (MQ=255)
                                                  aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                              >  2:1571612/1‑150 (MQ=255)
                                                  aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                              >  1:329935/1‑150 (MQ=255)
                                                  aaTCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                              >  2:204821/1‑150 (MQ=255)
                                                    tCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgatga                                                                                            >  2:1727171/1‑150 (MQ=255)
                                                    tCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgatga                                                                                            >  2:2070182/1‑150 (MQ=255)
                                                         cAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCt                                                                                       >  1:587495/1‑150 (MQ=255)
                                                              gcTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACCGACCACACACCCCAGCAGCCATACCCCAAAAGAATGACCTGATGAATTCTTTtct                                                                                  >  2:22440/1‑150 (MQ=255)
                                                                    ggcgtcGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCAc                                                                            <  2:1880492/150‑1 (MQ=255)
                                                                     gcggcgTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCAGCTTCGCGGGTTGACCGGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCGct                                                                           >  2:840061/1‑150 (MQ=255)
                                                                       ggcgTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCActct                                                                         >  1:423094/1‑150 (MQ=255)
                                                                          gTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGGCAGCACAGGCAACACCTCCACCCAGCCATCAAGATAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCa                                                                      >  1:1193961/1‑150 (MQ=255)
                                                                           tttGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCAc                                                                     <  1:2070182/150‑1 (MQ=255)
                                                                               ccAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTtctc                                                                                 <  2:1834983/134‑1 (MQ=255)
                                                                               ccAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTtctc                                                                                 >  1:1834983/1‑134 (MQ=255)
                                                                                cAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtct                                                                >  1:139553/1‑150 (MQ=255)
                                                                                 aaGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGcc                                                                                                                <  2:762878/101‑1 (MQ=255)
                                                                                 aaGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                                          >  2:2241254/1‑107 (MQ=255)
                                                                                 aaGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCaaaa                                                                                                          <  1:2241254/107‑1 (MQ=255)
                                                                                 aaGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGtctc                                                               >  1:393149/1‑150 (MQ=255)
                                                                                      caAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                                  >  1:504394/1‑62 (MQ=255)
                                                                                      caAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAgg                                                                                                                                                  <  2:504394/62‑1 (MQ=255)
                                                                                      caAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcc                                                          <  2:393149/150‑1 (MQ=255)
                                                                                       aaGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcct                                                         >  2:1531029/1‑150 (MQ=255)
                                                                                         gAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGtcctcc                                                       >  1:1670764/1‑150 (MQ=255)
                                                                                          aaaTAGCCGAAGGTGAGGTTGCCTGTTCTGGACGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCAAAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCa                                                      >  1:71794/1‑150 (MQ=255)
                                                                                            aTAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAAc                                                    >  2:1290756/1‑150 (MQ=255)
                                                                                                ccGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCg                                                >  2:599951/1‑150 (MQ=255)
                                                                                                ccGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCAGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCg                                                >  2:118532/1‑150 (MQ=255)
                                                                                                  gAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGtg                                              <  1:1727171/150‑2 (MQ=255)
                                                                                                       tGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAGCGGCCAGAACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAAtt                                         <  1:88447/150‑1 (MQ=255)
                                                                                                       tGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACAACTCCACCCAGCCATGCCCCAAAAGAATGACCTGGTGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAAtt                                         >  1:1219155/1‑150 (MQ=255)
                                                                                                       aggtgTTGCCTGTGCTGACCCTTTCGCGGGTAAACCCGCGAAGCGGCCAGAACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAAtt                                         >  1:368452/5‑150 (MQ=21)
                                                                                                        gAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTc                                        <  2:1458346/150‑1 (MQ=255)
                                                                                                         aGGTTGCCTGTTCTGGCCGCTTGGCGGAGTTAACGGACAGAGGGTTAGCATAGGCCACAACTTGAGACCGCCATGCCGCAAAATAATGAAATAATGAATTCTTGTCTCGGCATTTCCGAACGTCCCGGACCTGTAACCCCATAGAATTca                                       >  2:353275/1‑150 (MQ=255)
                                                                                                           gTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTa                                              <  2:2130781/141‑1 (MQ=255)
                                                                                                           gTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTa                                              >  1:2130781/1‑141 (MQ=255)
                                                                                                              gCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCaacaac                                  >  1:489504/1‑150 (MQ=255)
                                                                                                                cTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAccc                                >  1:1607097/1‑150 (MQ=255)
                                                                                                                   ttCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGc                                                                                                                      >  2:1380397/1‑61 (MQ=255)
                                                                                                                   ttCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGc                                                                                                                      <  1:1380397/61‑1 (MQ=255)
                                                                                                                   ttCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTtc                                                                                   <  2:2206826/96‑1 (MQ=255)
                                                                                                                   ttCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTtc                                                                                   >  1:2206826/1‑96 (MQ=255)
                                                                                                                        gCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGAACAGGCAACAGCTCCACTCTGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACg                        >  2:328984/1‑150 (MQ=255)
                                                                                                                         ccGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGGAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCTGCACTCTCCACCGTCTCTGTACTCCAACCCCGTAGAATTCAACAACCCCCACCACgg                       >  2:1352267/1‑150 (MQ=255)
                                                                                                                          cctttCGCGGGTAAACCCGCGAAGGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTctct                                                              <  2:1837406/107‑1 (MQ=38)
                                                                                                                          cctttCGCGGGTAAACCCGCGAAGGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTctct                                                              >  1:1837406/4‑110 (MQ=38)
                                                                                                                          cctttCGCGGGTAAACCCGCGAAGGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACggg                      >  2:1948228/4‑150 (MQ=255)
                                                                                                                           gCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAcccc                               >  2:145853/1‑140 (MQ=255)
                                                                                                                            tttCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGc                                                                                                                      <  2:2252258/51‑1 (MQ=255)
                                                                                                                            tttCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGc                                                                                                                      >  1:2252258/2‑52 (MQ=255)
                                                                                                                            tttCGCGGGTAAACCCGCGAAGGGGTCAGCACAGGCAAGACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCt                    >  2:1345729/2‑150 (MQ=255)
                                                                                                                             ttCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCtt                   >  1:538584/1‑150 (MQ=255)
                                                                                                                              tCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAcc                                 >  1:1115890/1‑135 (MQ=255)
                                                                                                                              tCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAAcc                                 <  2:1115890/135‑1 (MQ=255)
                                                                                                                               cgcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGa                 <  2:242366/150‑1 (MQ=255)
                                                                                                                                gcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                              >  1:1612662/1‑72 (MQ=255)
                                                                                                                                gcgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCtgat                                                                                              <  2:1612662/72‑1 (MQ=255)
                                                                                                                                 cggttaaaCCCGCGAAGCGGCCAGAACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATc               <  2:368452/143‑1 (MQ=255)
                                                                                                                                 cgGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATc               <  1:204821/150‑1 (MQ=255)
                                                                                                                                   ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa             <  2:415131/150‑1 (MQ=255)
                                                                                                                                   ggTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCACACGAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGa             <  1:22440/150‑1 (MQ=255)
                                                                                                                                          ccGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCg                                                <  2:2092359/108‑1 (MQ=255)
                                                                                                                                          ccGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCg                                                >  1:2092359/1‑108 (MQ=255)
                                                                                                                                          ccGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTcaa                                      >  1:1553964/1‑118 (MQ=255)
                                                                                                                                          ccGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTcaa                                      <  2:1553964/118‑1 (MQ=255)
                                                                                                                                          ccGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGAtttt      >  1:53664/1‑150 (MQ=255)
                                                                                                                                          ccGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGAtttt      >  2:733578/1‑150 (MQ=255)
                                                                                                                                          ccGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGAtttt      >  2:2229658/1‑150 (MQ=255)
                                                                                                                                            gcgAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGaa                                           <  2:328732/111‑1 (MQ=255)
                                                                                                                                            gcgAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGaa                                           >  1:328732/1‑111 (MQ=255)
                                                                                                                                             cgAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTATg   <  2:423094/150‑1 (MQ=255)
                                                                                                                                              gAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTATGc  >  1:535712/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
CTGGTGGAGGCCAATGACAAGCGCTGGGCGCGGGTGAAGGTGGTGCGCACAATCAACCAGGCGCTTGAGGCGGCGTTTGCCAAGCACAAGAAATAGCCGAAGGTGAGGTTGCCTGTTCTGGCCGCTTCGCGGGTTTACCCGCGAAAGGGTCAGCACAGGCAACACCTCCACCCAGCCATGCCCCAAAAGAATGACCTGATGAATTCTTTTCTCGGCACTCTCCACCGTCTCTGTCCTCCAACCCCGTAGAATTCAACAACCCCCACCACGGGCTTGATCGAGGATTTTATGC  >  NC_002947/1955118‑1955409

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: