Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,663,118 A→T 8.3% intergenic (+22/+16) 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,1180AT 8.3% 427.8 / 18.1 140intergenic (+22/+16)purT/yhjEphosphoribosylglycinamide formyltransferase/metabolite transport protein YhjE
Reads supporting (aligned to +/- strand):  ref base A (53/74);  new base T (5/7);  total (58/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 143 total aligned reads displayed.

CCTGTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGTG  >  NC_002947/1662970‑1663264
                                                                                                                                                    |                                                                                                                                                  
ccTGTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCaaa                                                                                                                                                   >  1:3849234/1‑150 (MQ=255)
  tGTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGc                                                                                                                                                 <  1:2961133/150‑1 (MQ=255)
   gTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCg                                                                                                                                                >  2:1004137/1‑150 (MQ=255)
    ttCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCgg                                                                                                                                               <  2:1039158/150‑1 (MQ=255)
    ttCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCgg                                                                                                                                               <  1:159442/150‑1 (MQ=255)
    ttCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCgg                                                                                                                                               <  2:2771648/150‑1 (MQ=255)
     tCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCACGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGc                                                                                                                                              <  2:4671655/150‑1 (MQ=255)
      cGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGcc                                                                                                                                             <  1:3913100/150‑1 (MQ=255)
       ggCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGccc                                                                                                                                            <  2:4531140/150‑1 (MQ=255)
       ggCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGccc                                                                                                                                            >  1:62368/1‑150 (MQ=255)
        gCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGcccc                                                                                                                                           <  2:4006150/150‑1 (MQ=255)
        gCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGcccc                                                                                                                                           <  2:2685339/150‑1 (MQ=255)
           aGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGt                                                                                                                                        <  2:62368/150‑1 (MQ=255)
            gCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGTCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTg                                                                                                                                       <  1:2063996/150‑1 (MQ=255)
             ccGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTgg                                                                                                                                      <  1:4173579/150‑1 (MQ=255)
              cGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGc                                                                                                                                     >  2:5069948/1‑150 (MQ=255)
               ggAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCt                                                                                                                                    >  1:132568/1‑150 (MQ=255)
               ggAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCt                                                                                                                                    >  1:1527373/1‑150 (MQ=255)
               ggAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGAGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCt                                                                                                                                    >  2:2877051/1‑150 (MQ=255)
                 aaaTCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTACTGGGGCCGCTTTGCGGCCCCAgcactctc                                                                                                                                  <  2:3555214/150‑8 (MQ=255)
                   aTCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACa                                                                                                                                <  2:4879794/150‑1 (MQ=255)
                   aTCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACa                                                                                                                                <  1:4639183/150‑1 (MQ=255)
                   aTCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACa                                                                                                                                <  1:4639158/150‑1 (MQ=255)
                       aTGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGt                                                                                                                            <  2:3622496/150‑1 (MQ=255)
                        tGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTg                                                                                                                           <  1:4248912/150‑1 (MQ=255)
                        tGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTg                                                                                                                           <  1:2732668/150‑1 (MQ=255)
                        tGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTACTGGGGCCGCTTTGCGGCCCCAgc                                                                                                                                        <  1:3501742/137‑2 (MQ=255)
                            aCCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTAt                                                                                                                       >  1:1617606/1‑150 (MQ=255)
                            aCCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTAt                                                                                                                       >  2:437439/1‑150 (MQ=255)
                              ccgccgCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTg                                                                                                                     >  1:2344671/1‑150 (MQ=255)
                              ccgccgCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTg                                                                                                                     >  1:1239601/1‑150 (MQ=255)
                                gccgcATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTgcg                                                                                                                   >  1:778422/1‑150 (MQ=255)
                                   gcATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGt                                                                                                                <  1:4112738/150‑1 (MQ=255)
                                      tGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGtctc                                                                                                             >  2:1238460/1‑150 (MQ=255)
                                      tGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGtctc                                                                                                             <  1:1432015/150‑1 (MQ=255)
                                      tGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGtctc                                                                                                             <  1:2773254/150‑1 (MQ=255)
                                          tgtgtgCCTGGCGCGTGACGACTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTc                                                                                                         <  2:3343612/150‑1 (MQ=255)
                                           gtgtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCa                                                                                                        <  1:529010/150‑1 (MQ=255)
                                           gtgtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCa                                                                                                        <  1:986857/150‑1 (MQ=255)
                                           gtgtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCa                                                                                                        <  1:596646/150‑1 (MQ=255)
                                           gtgtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCa                                                                                                        <  2:2516355/150‑1 (MQ=255)
                                             gtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAgg                                                                                                      <  1:1324581/150‑1 (MQ=255)
                                             gtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAgg                                                                                                      >  2:3491944/1‑150 (MQ=255)
                                             gtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAgg                                                                                                      >  2:3492226/1‑150 (MQ=255)
                                             gtgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAgg                                                                                                      >  1:787628/1‑150 (MQ=255)
                                              tgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGc                                                                                                     <  1:3321030/150‑1 (MQ=255)
                                              tgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGc                                                                                                     >  2:184423/1‑150 (MQ=255)
                                              tgCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGc                                                                                                     <  1:3321425/150‑1 (MQ=255)
                                                ccTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGcac                                                                                                   <  2:2649089/150‑1 (MQ=255)
                                                 cTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTCTGGCGCGTCTCTTTCAGGcaca                                                                                                  >  1:1869306/1‑150 (MQ=255)
                                                  tGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAg                                                                                                 >  2:4243317/1‑150 (MQ=255)
                                                   ggCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGc                                                                                                <  1:493246/150‑1 (MQ=255)
                                                    gcgcGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCa                                                                                               >  1:3036795/1‑150 (MQ=255)
                                                     cgcgTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCAc                                                                                              <  1:3097834/150‑1 (MQ=255)
                                                     cgcgTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCAc                                                                                              <  2:3802569/150‑1 (MQ=255)
                                                           aCGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAAt                                                                                        >  2:578576/1‑150 (MQ=255)
                                                            cGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATc                                                                                       <  2:1418003/150‑1 (MQ=255)
                                                             gAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGcccc                                                                                                                                           <  2:312153/97‑1 (MQ=255)
                                                             gAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGcccc                                                                                                                                           >  1:312153/1‑97 (MQ=255)
                                                             gAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCTTTGCGGcccc                                                                                                                                           >  1:311935/1‑97 (MQ=255)
                                                             gAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCa                                                                                      <  2:3974657/150‑1 (MQ=255)
                                                               aTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAgg                                                                                    >  2:1092191/1‑150 (MQ=255)
                                                               aTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAgg                                                                                    <  1:4998094/150‑1 (MQ=255)
                                                               aTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAgg                                                                                    <  1:3578536/150‑1 (MQ=255)
                                                               aTCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAgg                                                                                    <  1:2144026/150‑1 (MQ=255)
                                                                tCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGc                                                                                   <  1:2346608/150‑1 (MQ=255)
                                                                 cGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCt                                                                                  <  2:1763119/150‑1 (MQ=255)
                                                                 cGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCt                                                                                  <  1:3507070/150‑1 (MQ=255)
                                                                   gTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGa                                                                                >  1:537318/1‑150 (MQ=255)
                                                                    tCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGAt                                                                               <  1:4352194/150‑1 (MQ=255)
                                                                     cGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATc                                                                              <  2:758343/150‑1 (MQ=255)
                                                                       aaCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAc                                                                            >  1:1035057/1‑150 (MQ=255)
                                                                       aaCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAc                                                                            >  1:1034964/1‑150 (MQ=255)
                                                                       aaCTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAc                                                                            >  1:2569521/1‑150 (MQ=255)
                                                                          tGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAccgc                                                                         >  2:947129/1‑150 (MQ=255)
                                                                          tGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCAccgc                                                                         >  2:2464890/1‑150 (MQ=255)
                                                                          tGGCGCGGCCCAAGGCGACCCGCGCTTCGCAGGCGGTCAACGTCGAGTTCACATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACTGAAATCAGCCTCATCAcagc                                                                         >  2:2464642/1‑150 (MQ=255)
                                                                             cgcgGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACcgccgc                                                                      >  2:4574048/1‑150 (MQ=255)
                                                                              gcgGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCgccgcc                                                                     >  1:2751412/1‑150 (MQ=255)
                                                                                     aaGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGaca                                                              >  2:1267463/1‑150 (MQ=255)
                                                                                        gCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACAt                                                           <  1:3672948/150‑1 (MQ=255)
                                                                                        gCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACAt                                                           <  1:262325/150‑1 (MQ=255)
                                                                                         cGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCACAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATa                                                          >  2:1401548/1‑150 (MQ=255)
                                                                                          gACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTCCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATaa                                                         <  2:4846458/150‑1 (MQ=255)
                                                                                            cccGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTCATCACCGCCGCCACCGACACAtcggg                                                       >  1:323087/1‑146 (MQ=255)
                                                                                              cgcgCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAAcccc                                                     <  2:2942283/150‑1 (MQ=255)
                                                                                                                  ggTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAACcccccc                                 <  2:4013503/150‑1 (MQ=255)
                                                                                                                  ggTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAAcccccc                                 <  1:4829005/150‑1 (MQ=255)
                                                                                                                    tCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGGGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCCCAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCCCCCCGCGCGAAccccccgg                               >  2:2228586/1‑148 (MQ=255)
                                                                                                                        gTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTc                           <  1:3065446/150‑1 (MQ=255)
                                                                                                                            ggctCCGAGAGTGCTGGGGCCGCACCGCGGCCCCAGCGGCTTCCAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCAGGCTCAccc                       <  2:4250016/147‑2 (MQ=255)
                                                                                                                                   agagTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTg                <  1:3617040/150‑1 (MQ=255)
                                                                                                                                   agagTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTg                <  1:4347780/150‑1 (MQ=255)
                                                                                                                                    gagTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTgg               <  1:947129/150‑1 (MQ=255)
                                                                                                                                    cagcgccgggcCCGCAACCCCGCCCCAGGGGCTTACACCGTGGTATTGCGCGTCTCGTTCCGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTgg               <  2:4490008/143‑1 (MQ=255)
                                                                                                                                          tGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCATGGCAATCAGGCTGATCACCGCCGCCACCGACACATGACCCCCCCCCCAGCGCAAGCCCCCCGGGCTCACCAGCTGCGGGGCAAAg         >  2:3881675/1‑150 (MQ=255)
                                                                                                                                          tGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCAGCCAGCGCAAACCCCCCCTGCTCACCAGCGTCTGGGcgagg         >  2:4732137/1‑146 (MQ=255)
                                                                                                                                                cGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCGGCTCAAGCCCCCCCTGCGCACCACCTGCGGGGCCAAGAAGGGt   >  2:3011437/1‑150 (MQ=255)
                                                                                                                                                cGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCGCACACCCCCCCGGCTCCCCCGCGTCGGGGCCCAGGACGGt   >  2:4843260/1‑150 (MQ=255)
                                                                                                                                                 gCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATACCCCCCCACCCAGCTCACTCCCCCCATGCTCACCCGCTGCTGGGCGAAGGACGGtg  >  2:1204941/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CCTGTTCGGCAAGCCGGAAATCAATGGCACCCGCCGCATGGGTGTGTGCCTGGCGCGTGACGAATCGGTCGAACTGGCGCGGGCCAAGGCGACCCGCGCTTCGCAGGCGGTCAAGGTCGAGTTCTGATCCGAGAGTGCTGGGGCCGCAAAGCGGCCCCAGTGGCTTACAGCGTGGTATTGCGCGTCTCTTTCAGGCACAGCACGGCAATCAGGCTGATCACCGCCGCCACCGACACATAACCCCCCACCCAGCTCAAACCCCCCATGCTCACCAGCTTCTGGGCAAAGAACGGTG  >  NC_002947/1662970‑1663264

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