Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,745,790 T→G 9.2% intergenic (+32/+720) gltA → / ← yeiW citrate synthase/putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,745,7900TG 9.2% 683.6 / 47.7 244intergenic (+32/+720)gltA/yeiWcitrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (92/125);  new base G (11/12);  total (104/140)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.12e-01

Only 100 of 245 total aligned reads displayed.

GGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTG  >  NC_002947/4745641‑4745938
                                                                                                                                                     |                                                                                                                                                    
gggCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTATAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCt                                                                                                                                                      <  2:7889985/150‑1 (MQ=255)
 ggCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTc                                                                                                                                                     >  1:1084636/1‑150 (MQ=255)
  gCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCg                                                                                                                                                    >  2:357842/1‑150 (MQ=255)
   cTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCgg                                                                                                                                                   <  2:803012/150‑1 (MQ=255)
    tGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCggg                                                                                                                                                  >  1:2679525/1‑150 (MQ=255)
    tGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCggg                                                                                                                                                  <  1:3327371/150‑1 (MQ=255)
      gATCTCGCACTGGAAAGATATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCa                                                                                                                                                >  2:4698062/1‑150 (MQ=255)
      gATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCa                                                                                                                                                >  1:8434043/1‑150 (MQ=255)
        tctcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCGAGGGCAGc                                                                                                                                              <  1:1671646/150‑1 (MQ=255)
         ctcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGcc                                                                                                                                             >  1:8077346/1‑150 (MQ=255)
          tcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                            >  1:2812177/1‑150 (MQ=255)
          tcGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCt                                                                                                                                            >  2:3789638/1‑150 (MQ=255)
           cGCACTGGAAAGATATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtt                                                                                                                                           <  1:4698062/150‑1 (MQ=255)
            gCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttt                                                                                                                                          <  1:3242525/150‑1 (MQ=255)
             cACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttt                                                                                                                                         <  2:4264425/150‑1 (MQ=255)
              ccgggAACGAACGGCTCGCCCGCCCGTACACGCTCGGCCGTCCGCGCCAGCTGTCCACCGGCGAGCAGAAGCGCGACATCGGTGCCCGGCAAGACCGCAAGTAAGCTGCAGGAGCCGAACGCAAAAAGGCTGCCCGCGGGCAGCcttctt                                                                                                                                        <  2:1684199/147‑1 (MQ=255)
              aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCttttt                                                                                                                                        >  2:7525780/1‑150 (MQ=255)
              aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGACAGCCttttt                                                                                                                                        >  2:8391196/1‑150 (MQ=255)
              aCTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCcttctt                                                                                                                                        <  2:1826829/150‑1 (MQ=255)
               cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCtttttt                                                                                                                                       <  2:2271605/150‑1 (MQ=255)
               cTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAAAAGGCTGCCCGAGGGCAGCcttcttg                                                                                                                                       <  1:3579341/150‑2 (MQ=255)
                 ggAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAt                                                                                                                                     <  1:5278943/150‑1 (MQ=255)
                  gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAtt                                                                                                                                    <  1:2879339/150‑1 (MQ=255)
                  gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTGTTAtt                                                                                                                                    >  1:5982648/1‑150 (MQ=255)
                  gaaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTGTTTAtt                                                                                                                                    >  1:3074538/1‑150 (MQ=255)
                   aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAttt                                                                                                                                   <  1:6373203/150‑1 (MQ=255)
                   aaagaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTAttt                                                                                                                                   <  1:7708706/150‑1 (MQ=255)
                      gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGACGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTc                                                                                                                                <  2:4669744/150‑1 (MQ=255)
                      gaaaTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTc                                                                                                                                <  2:6769919/150‑1 (MQ=255)
                         aTGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCtg                                                                                                                             >  1:8679048/1‑150 (MQ=255)
                          tGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCtgt                                                                                                                            <  2:4423165/150‑1 (MQ=255)
                              ctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCg                                                                                                                        <  2:3514760/150‑1 (MQ=255)
                              ctcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCg                                                                                                                        <  1:3642791/150‑1 (MQ=255)
                               tcCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCTAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTGTATTTGTCCTGTGCCgg                                                                                                                       >  2:2568286/1‑150 (MQ=255)
                                    cccGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCctct                                                                                                                  <  2:4894702/150‑1 (MQ=255)
                                         aCAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTcgcg                                                                                                             <  1:3769476/150‑1 (MQ=255)
                                          cAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCACTTCGCgg                                                                                                            >  1:5591844/1‑150 (MQ=255)
                                           aaGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGc                                                                                                           <  1:887867/150‑1 (MQ=255)
                                             gATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCt                                                                                                         <  1:1368146/150‑1 (MQ=255)
                                              aTCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCtt                                                                                                        <  2:5359065/150‑1 (MQ=255)
                                               tCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCgggcaag                                                                                                       <  2:7157102/150‑6 (MQ=255)
                                                  gCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGccc                                                                                                    <  1:7294707/150‑1 (MQ=255)
                                                    cGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCcgc                                                                                                  <  1:6780591/150‑1 (MQ=255)
                                                       ccGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                               <  1:7736659/150‑1 (MQ=255)
                                                        cgcgCCAGCTGTACACCGGCGAGTAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                              <  1:927346/150‑1 (MQ=255)
                                                        cgcgCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAag                                                                                              <  1:8237297/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  2:4067577/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  2:5135778/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  1:7120608/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAaga                                                                                             <  1:7013818/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                             <  1:5818092/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                             <  1:6115857/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                             <  1:4622080/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                             <  1:4124241/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAaga                                                                                             <  1:7902075/150‑1 (MQ=255)
                                                         gcgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCGTCGCGGCCGTGCCCGCGAaga                                                                                             >  2:3142893/1‑150 (MQ=255)
                                                          cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTTCTGTGCCGGCCTCTTCGCGGGCAAGGCCGCGAagag                                                                                            <  1:7016726/150‑1 (MQ=39)
                                                          cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                            <  1:2141336/150‑1 (MQ=255)
                                                          cgcCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAagag                                                                                            <  2:6228095/150‑1 (MQ=255)
                                                            ccAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTCGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGc                                                                                          <  1:1910579/150‑1 (MQ=255)
                                                              aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCg                                                                                        <  2:6010670/150‑1 (MQ=255)
                                                              aGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGACCCCGAAGAGGCCg                                                                                        >  1:4732291/1‑150 (MQ=255)
                                                               gCTGTACACCGGCGAGCAGAAGCGCGACAGCGTGGCCCTGAAAGACCGCAAGTAAGCTGCCGTAGCCGAACACAAGACGGCGGCCCTGGGGCAGCCTTTTTTATGTGTCGTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCgg                                                                                       <  2:6183521/150‑1 (MQ=255)
                                                                cTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGAACGCAAGTAAGCTGCAGTAGCAGAACGCAAGAAGGCTGCCCACCGGCAGACATTTTTATTTGTCCTGTGCCCGCCAAATACCGGCCCATCCCCCGAAAAGGCCGGc                                                                                      >  1:6313844/1‑150 (MQ=255)
                                                                 tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                               >  2:6259114/1‑140 (MQ=255)
                                                                 tGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGaa                                                                                               <  1:6259114/140‑1 (MQ=255)
                                                                        cGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTGATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGCAGAGGCCGGCACAGacac                                                                              >  2:8219762/1‑150 (MQ=255)
                                                                            gagaAGACGCGCGACATCGTGGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCGGGCCTCTTCGCGGCCTTGCCCGCGAAGTGGCCGGCACAGACACcaca                                                                          <  2:7772806/150‑1 (MQ=255)
                                                                                   aGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTCGCCGAACGAAAGATGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGGGGCCTTGCCCGCGAAGTGGCCGGCACAGACACCACAAAACCTc                                                                   <  1:6586344/150‑1 (MQ=255)
                                                                                    gcgcGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTAGCGGCCTTGCCCGCGAAGAGGCCGGCACCGACACCACAAAACCTCa                                                                  >  2:3653107/1‑150 (MQ=255)
                                                                                      gcgACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCGTAGCGGCCATGCCCCCGCAGAGGCCCGCACACACACCCCaaaaacacacg                                                                >  1:764690/1‑143 (MQ=255)
                                                                                          cATCGTTGCCCTGAATGACCGCAAGTAAGGTGCTGTTGCCGAACGCAAGATGGGTGCCCTCGGGCAGCCTTTTTTATTTGTCGTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGTGGCCGGCACAGACACCACAAAACCTCAGttttt                                                            <  2:5822207/150‑1 (MQ=255)
                                                                                           aTCGTTGCCCTGAAAGACCGAAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCATTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAATAGGCCGGCACAGACACCACAAAACCTCAGTTTTtc                                                           <  2:4172742/150‑1 (MQ=255)
                                                                                                tGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCACGAAGGCTGCCCACGGCCAGCATTTTTTATTTGTCCTGTGCCGGCATCTAGGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTccgc                                                      <  2:6058998/150‑1 (MQ=255)
                                                                                                  cccTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCACAGTTTTTATccgccc                                                    >  1:3482841/1‑150 (MQ=255)
                                                                                                     tGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCATCGCGGCCGTGCCCGCGCAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTATccgcccgcg                                                 >  2:1728971/1‑149 (MQ=255)
                                                                                                              gCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCATCGCGGCCGTGCCCGCGACGAGGCCGGCACAGACACCACAAAACCACAGTTTTTGTGCGCCCGCACACCCCGcc                                        >  1:1231508/1‑150 (MQ=255)
                                                                                                                   tAAGCTACAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCATGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGCCCGTTc                                   >  2:1320265/1‑150 (MQ=255)
                                                                                                                    aaGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                  <  2:8734196/150‑1 (MQ=255)
                                                                                                                    aaGCTGCAGTAGCCGAACGCAAGAAGGCGGCCCCCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTGGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCa                                  <  2:1980289/150‑1 (MQ=255)
                                                                                                                      gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTGCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAgg                                <  2:5591844/150‑1 (MQ=255)
                                                                                                                      gCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTAGCGGCCTTGCCCGCGCAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTATACCCCCGCGCACCCAGcccggtccgg                                >  2:2920008/1‑143 (MQ=255)
                                                                                                                          cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTAt                            >  2:2111930/1‑150 (MQ=255)
                                                                                                                          cAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTAAGTTTTTCTCCGCCCGCTCAAGCAGCCCATTCAGGGTAt                            >  2:6420937/1‑150 (MQ=255)
                                                                                                                             tAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCCCGGCCTTTCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCACCGCCCGCTCACGCAGCCCGGTCAGGGGAATGa                         >  2:5510816/1‑150 (MQ=255)
                                                                                                                                  gAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgc                    <  2:2812177/150‑1 (MQ=255)
                                                                                                                                    aataaaaaagGCTGCCCGAGGGCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGgcgc                  >  2:8680894/8‑150 (MQ=255)
                                                                                                                                    aataaaaaagGCTGCCCGAGGGCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCGTTCCGGGTATTGAACGgcgc                  >  2:8443618/8‑150 (MQ=255)
                                                                                                                                      gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACCAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCAt                >  2:5097845/1‑150 (MQ=255)
                                                                                                                                      gCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCATTCAGGGTATTGAACGGCGCAt                >  2:1444960/1‑150 (MQ=255)
                                                                                                                                       cAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGGCCTGTGCCGCCCTCTGCGCGGCCTTGCCCGCGAAGCGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACgggctgga               <  2:1990714/150‑7 (MQ=255)
                                                                                                                                       aaaaaagGCTGCCCGAGGGCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATc               >  1:3554162/5‑150 (MQ=255)
                                                                                                                                        aagaagGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCGCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCGTTCAGGGTATTGAACGGCGCATcc              >  2:8196859/1‑150 (MQ=255)
                                                                                                                                        aaaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGCAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCGGTCAGGGTATTGAACGGCGCCTcc              >  2:7075825/4‑150 (MQ=255)
                                                                                                                                        aaaaagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCGTCGCGGCCTTGCCCGCGAAGGGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCGCACGCAGCCCTTTCAGGGTATTGAACGGCGCAtac              >  1:5725339/4‑148 (MQ=255)
                                                                                                                                           aagGCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGGCTCTTCGCGGCCGTGCCCGCGCAGAGGCCGGCCCAGACACCACAAAAACGCAGTTTTTGTCCGCCCGCTCACGCAGCCCGTTCACGGGATTGAACCGCGCATccacc           >  1:8521880/1‑150 (MQ=255)
                                                                                                                                             ggCTGCCCGAGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCGTGCCCGCGAAGAGGCCGGCCCAGACACCACAAAACCACAGTTTTTGTCCGCCCGCGCAAGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac         >  2:3340352/1‑150 (MQ=255)
                                                                                                                                             ggCTGCCCGAGGGCAGCCTTCTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCcaccac         >  2:7316763/1‑150 (MQ=255)
                                                                                                                                               cTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGTGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGa       <  2:7685433/150‑1 (MQ=255)
                                                                                                                                                 gCCCTCTGGATGCCTTTTTTATTTGTCCTGTTCCGGCCTATTAGCGGCCTTGCCCGCTAATAGTCCTGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCTCTCAAGCAGCCCTTTCAGGGTATTGAACGGCTCATCCACCACTAAc     <  2:6538829/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGCTGGATCTCGCACTGGAAAGAAATGCTCTCCAGCCCGTACAAGATCGGCCGTCCGCGCCAGCTGTACACCGGCGAGCAGAAGCGCGACATCGTTGCCCTGAAAGACCGCAAGTAAGCTGCAGTAGCCGAACGCAAGAAGGCTGCCCTCGGGCAGCCTTTTTTATTTGTCCTGTGCCGGCCTCTTCGCGGCCTTGCCCGCGAAGAGGCCGGCACAGACACCACAAAACCTCAGTTTTTCTCCGCCCGCTCACGCAGCCCTTTCAGGGTATTGAACGGCGCATCCACCACGAACTTG  >  NC_002947/4745641‑4745938

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