Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,663,117 A→T 5.0% intergenic (+21/+17) purT → / ← yhjE phosphoribosylglycinamide formyltransferase/metabolite transport protein YhjE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,663,1170AT 5.0% 239.3 / 3.5 80intergenic (+21/+17)purT/yhjEphosphoribosylglycinamide formyltransferase/metabolite transport protein YhjE
Reads supporting (aligned to +/- strand):  ref base A (32/44);  new base T (2/2);  total (34/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.60e-01

GCCTGTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGTGCT  >  NC_002947/1662969‑1663266
                                                                                                                                                    |                                                                                                                                                     
gCCTGTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCaa                                                                                                                                                      >  1:2402234/1‑150 (MQ=255)
gCCTGTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCaa                                                                                                                                                      >  1:1404622/1‑150 (MQ=255)
      tCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGc                                                                                                                                                <  1:846731/150‑1 (MQ=255)
          cAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGCCCCa                                                                                                                                            <  1:265294/150‑1 (MQ=255)
            aGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGt                                                                                                                                          <  2:1360264/150‑1 (MQ=255)
            aGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGt                                                                                                                                          <  1:926010/150‑1 (MQ=255)
             gCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTg                                                                                                                                         <  2:847913/150‑1 (MQ=255)
             gCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCCCAGCGCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTg                                                                                                                                         <  1:1733248/150‑1 (MQ=255)
                ggAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGGTGCTTGTGTTCTGGGGCCGCAAAGCGGCCCCAGTGGCt                                                                                                                                      >  2:2017538/1‑150 (MQ=255)
                ggAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCt                                                                                                                                      <  1:942531/150‑1 (MQ=255)
                  aaaTCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGCATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTa                                                                                                                                    <  1:1455255/150‑1 (MQ=255)
                        aTGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGt                                                                                                                              <  2:1963884/150‑1 (MQ=255)
                        aTGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGt                                                                                                                              >  1:2216125/1‑150 (MQ=255)
                        aTGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGt                                                                                                                              >  1:305094/1‑150 (MQ=255)
                         tGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGTCATCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTg                                                                                                                             >  1:22997/1‑150 (MQ=255)
                         tGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTg                                                                                                                             <  2:2413756/150‑1 (MQ=255)
                         tGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTg                                                                                                                             <  1:1196230/150‑1 (MQ=255)
                            cACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTa                                                                                                                          <  2:1404622/150‑1 (MQ=255)
                              cccGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGTTCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTAtt                                                                                                                        <  2:393500/150‑1 (MQ=255)
                               ccgccgCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGCCCCa                                                                                                                                            <  1:1238329/129‑1 (MQ=255)
                               ccgccgCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTg                                                                                                                       >  2:2372365/1‑150 (MQ=255)
                               ccgccgCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGACGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGCCCCa                                                                                                                                            >  2:1238329/1‑129 (MQ=255)
                                  ccgcATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTgcgc                                                                                                                    <  1:2044636/150‑1 (MQ=255)
                                  ccgcATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTAGCTTACAGCGTGGTATTgcgc                                                                                                                    <  1:2423900/150‑1 (MQ=255)
                                        gggTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTctct                                                                                                              >  1:2102356/1‑150 (MQ=255)
                                        gggTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTctct                                                                                                              <  1:1572359/150‑1 (MQ=255)
                                           tgtgtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTc                                                                                                           >  1:1495244/1‑150 (MQ=255)
                                            gtgtgCCTGGGGCGTGGCGAAGCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCgttta                                                                                                          <  2:748142/150‑3 (MQ=255)
                                                gCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGca                                                                                                      <  1:2372365/150‑1 (MQ=255)
                                                 ccTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGcac                                                                                                     <  1:2074584/150‑1 (MQ=255)
                                                        cgTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACgg                                                                                              <  1:2182005/150‑1 (MQ=255)
                                                               aaTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAg                                                                                       >  1:2436516/1‑150 (MQ=255)
                                                                aTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAgg                                                                                      <  1:2348929/150‑1 (MQ=255)
                                                                 tCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTg                                                                                                                             <  1:2397155/110‑1 (MQ=255)
                                                                 tCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTg                                                                                                                             >  2:2397155/1‑110 (MQ=255)
                                                                    gTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTGTGATCCGAGAGGGCGGGGGCCGCAAAGCGGCCCCAGTGGCTTCCAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGTAATCAGGCTGa                                                                                  >  2:1257102/1‑150 (MQ=255)
                                                                    gTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAAGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTAGAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGa                                                                                  >  2:1230012/1‑150 (MQ=255)
                                                                        aaCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGACCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAgg                                                                                      <  2:1341406/142‑1 (MQ=255)
                                                                        aaCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGACCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAgg                                                                                      >  1:1341406/1‑142 (MQ=255)
                                                                         aCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAcc                                                                             <  1:217527/150‑1 (MQ=255)
                                                                         aCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGGCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAcc                                                                             <  1:265075/150‑1 (MQ=255)
                                                                           tGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAccgc                                                                           >  1:209735/1‑150 (MQ=255)
                                                                            ggCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAccgcc                                                                          >  1:1007127/1‑150 (MQ=255)
                                                                               gcgGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGcaca                                                                                                    >  1:901404/1‑121 (MQ=255)
                                                                               gcgGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGcaca                                                                                                    <  2:901404/121‑1 (MQ=255)
                                                                                    ccAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCTAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGAAAAGCACGGCAATCAGGCTGATCACCGCCGCCACCGa                                                                  <  1:2017538/150‑1 (MQ=255)
                                                                                    ccAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGc                                                                                                                                       <  2:666296/81‑1 (MQ=255)
                                                                                    ccAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGc                                                                                                                                       >  1:666296/1‑81 (MQ=255)
                                                                                     cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAGAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGac                                                                 >  2:142261/1‑150 (MQ=255)
                                                                                     cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCa                                                                                                                                            <  2:586046/75‑1 (MQ=255)
                                                                                     cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCa                                                                                                                                            >  1:586046/1‑75 (MQ=255)
                                                                                     cAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGac                                                                 >  1:1567957/1‑150 (MQ=255)
                                                                                           gACCCGCGCTTCGCCGGCGGCCCCGGTCGCATTCTCATCCTGGGGTGCGGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATaa                                                           <  2:22997/150‑1 (MQ=255)
                                                                                                gcgcTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAgg                                                                                                        <  1:1181858/100‑1 (MQ=255)
                                                                                                gcgcTTCGCAGGCGGTCAAGGTCAAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAgg                                                                                                        >  2:1181858/1‑100 (MQ=255)
                                                                                                        cAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATACCCCcccacacggg                                              >  1:1604848/1‑145 (MQ=255)
                                                                                                           gCGGTCAAGGTCGAGTTCTGATCCGCGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGCGGTATTGCGCGACACTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCa                                           >  2:918483/1‑150 (MQ=255)
                                                                                                             ggTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTTCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGaca                                                                >  2:1433546/1‑127 (MQ=255)
                                                                                                             ggTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTTCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGaca                                                                <  1:1433546/127‑1 (MQ=255)
                                                                                                              gTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACAt                                                             <  2:1967655/129‑1 (MQ=255)
                                                                                                              gTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACAt                                                             >  1:1967655/1‑129 (MQ=255)
                                                                                                                 aaGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCTCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAAcccc                                     <  1:551395/150‑1 (MQ=255)
                                                                                                                 aaGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAAcccc                                     <  2:1567957/150‑1 (MQ=255)
                                                                                                                    gTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCa                                  >  1:955342/1‑150 (MQ=255)
                                                                                                                    gTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCa                                  <  2:209735/150‑1 (MQ=255)
                                                                                                                    gTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCa                                  <  2:2402234/150‑1 (MQ=255)
                                                                                                                    gTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCGCAAACCCCCCa                                  >  1:1821992/1‑150 (MQ=255)
                                                                                                                     tCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACAAACCACCCATCTCAAACACACCAt                                 >  1:71947/1‑150 (MQ=255)
                                                                                                                       gAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGc                                                                                     >  1:2452033/1‑96 (MQ=255)
                                                                                                                       gAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGc                                                                                     <  2:2452033/96‑1 (MQ=255)
                                                                                                                          ttCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCa                            <  2:2216125/150‑1 (MQ=255)
                                                                                                                          ttCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGATGATGGCCGCCTCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCa                            <  2:809997/150‑1 (MQ=255)
                                                                                                                              gATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCTGCCACCGACACATAACCGCCCACCCGGCCTAAACCCCCCATGCccaccat                        >  1:1317193/1‑149 (MQ=255)
                                                                                                                               aTCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGc                       <  2:1821992/150‑1 (MQ=255)
                                                                                                                                   gagaGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCGACCCAGCACACACCCCCCATACCcagcagactcc                   >  2:297263/1‑145 (MQ=255)
                                                                                                                                   gagaGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCa                         <  1:72470/144‑1 (MQ=255)
                                                                                                                                   gagaGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCACGCTCACCa                         >  2:72470/1‑144 (MQ=255)
                                                                                                                                    agagTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGAAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCGCAAACCCCCCATGCTCACCAGCTTCTg                  <  2:2102356/150‑1 (MQ=255)
                                                                                                                                     gagTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACCGCGCGGCATTCCGCGTACCCTTCAGGCTCACCCCGTGATCCAGTCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTgg                 <  2:1462870/150‑1 (MQ=255)
                                                                                                                                     gagTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTgg                 <  2:305094/150‑1 (MQ=255)
                                                                                                                                        cacTGGGGCCGCTTTGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCCCCCCGCTCAAACCCCCCATCCCCACCACCTTCTGGCCa              >  2:88857/3‑150 (MQ=255)
                                                                                                                                             gggCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCGCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGaa         >  1:987686/1‑150 (MQ=255)
                                                                                                                                                    aaaGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGtgct  >  1:2240374/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCTGTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGTGCT  >  NC_002947/1662969‑1663266

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