Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,344,121 A→C 5.1% C8G (TGC→GGC)  PP_1168 ← TRAP dicarboxylate transporter subunit DctQ

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,344,1210AC 5.1% 453.6 / 10.9 138C8G (TGC→GGC) PP_1168TRAP dicarboxylate transporter subunit DctQ
Reads supporting (aligned to +/- strand):  ref base A (69/62);  new base C (2/5);  total (71/67)
Fisher's exact test for biased strand distribution p-value = 2.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 139 total aligned reads displayed.

GCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGA  >  NC_002947/1343972‑1344264
                                                                                                                                                     |                                                                                                                                               
  gggCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCt                                                                                                                                               <  2:2345281/150‑1 (MQ=255)
  gggCGAAGATGCCCCACGGGATGATCAGCGTCAGGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCt                                                                                                                                               <  1:651522/150‑1 (MQ=255)
    gCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGc                                                                                                                                             <  1:411605/150‑1 (MQ=255)
     cGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGTGCGGCGCCCGCCGCCCTGcc                                                                                                                                            <  2:826574/150‑1 (MQ=255)
      gAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCg                                                                                                                                           >  1:2925505/1‑150 (MQ=255)
          aTGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGt                                                                                                                                       >  2:132312/1‑150 (MQ=255)
          aTGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGt                                                                                                                                       >  2:132186/1‑150 (MQ=255)
              cccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATCGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAt                                                                                                                                   >  2:2583865/1‑150 (MQ=255)
               ccACGGGTTGATCAGCGTCTTGGCCAGGATCGACAGGCCCGCGATCCTGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGAtt                                                                                                                                  <  2:2158024/150‑1 (MQ=255)
               ccACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGGCGAtt                                                                                                                                  >  1:3176973/1‑150 (MQ=255)
                 agggtttgatCAGGGTCATGTCCAGGATCGACTGGCCCGCTATCAATATGCAACTGGGGTACAGCTTTTCGTTCACGCTCTGGAAAAGCGTTTTCATTGGCATCACAAAAGCGGCTGGGCGGCGGGGGCCGCACTGCCGAGGTCGATTgg                                                                                                                                <  1:4079807/144‑1 (MQ=255)
                 aCGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTgg                                                                                                                                <  1:2124012/150‑1 (MQ=255)
                    ggATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCa                                                                                                                             >  2:3355065/1‑150 (MQ=255)
                       ttatCAGCGTCTTGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCTGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAgg                                                                                                                          <  1:1643118/148‑1 (MQ=255)
                        gatCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAggg                                                                                                                         >  1:3768105/1‑150 (MQ=255)
                          tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcg                                                                                                                       >  2:305588/1‑150 (MQ=255)
                          tCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcg                                                                                                                       >  2:3298175/1‑150 (MQ=255)
                           cAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcgg                                                                                                                      <  2:4966588/150‑1 (MQ=255)
                            aGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACCGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATCGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcggc                                                                                                                     <  2:3623239/150‑1 (MQ=255)
                            aGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcggc                                                                                                                     <  2:3239880/150‑1 (MQ=255)
                            aGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcggc                                                                                                                     >  1:2686875/1‑150 (MQ=255)
                            aGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGggcggc                                                                                                                     >  2:1212408/1‑150 (MQ=255)
                              cGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTg                                                                                                                   <  2:2067395/150‑1 (MQ=255)
                                   gggCCCGGACCGACTGGCCCGCGAGCCAGTTGCACCGGCGGTCCCGGGTGTCGTTCACGCTCAGCAAAAGCTATTTCATCGGCATCACCAACGCGGCCGGGCGGCGGGCGCCGCACGGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACt                                                                                                              <  1:4746802/149‑1 (MQ=255)
                                     gCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTgg                                                                                                            >  2:140693/1‑150 (MQ=255)
                                       cAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGAc                                                                                                          >  2:2129204/1‑150 (MQ=255)
                                        aGGATCGACAGGCCCGCGATCCAGATGCACCTGCGGTACAGCGTGTCGTTCACGCTCAGGAACAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACg                                                                                                         <  2:4605076/150‑1 (MQ=255)
                                            tCGACCGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTGGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAACGCGGCTGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGt                                                                                                     <  1:4035926/150‑1 (MQ=255)
                                             cGACAGGCCCGCTATCCAGATGCAACTGCGTTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCACCACCACAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTc                                                                                                    <  1:2665980/150‑1 (MQ=255)
                                              gACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCg                                                                                                   >  2:2236761/1‑150 (MQ=255)
                                              gACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATGTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACGGGACGGCGTCg                                                                                                   >  2:2236541/1‑150 (MQ=255)
                                                 aGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATa                                                                                                <  2:2012062/150‑1 (MQ=255)
                                                   gCCCGCTATCCCGTTGCAACGGCGGTACCGCGTGTCGTCCCCGCTCAGTAAAAGCGATTACATCGGCATCACAAACGCGGCCGGGGGGCGGGCGCCGCACGGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACg                                                                                              <  1:995506/150‑1 (MQ=255)
                                                    cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAACGCGGCCGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                             <  1:4903401/150‑1 (MQ=255)
                                                    cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                             >  2:4667925/1‑150 (MQ=255)
                                                    cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                             <  2:3224822/150‑1 (MQ=255)
                                                    cccGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGc                                                                                             >  2:802759/1‑150 (MQ=255)
                                                      cgcgATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCtt                                                                                           >  2:3961192/1‑150 (MQ=255)
                                                       gcgATCCTGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTtc                                                                                          <  2:1845705/150‑1 (MQ=255)
                                                          aaccAGATGCAACTGCGGTACCGCGTGTCGTTCACGCTCAGGAAATGCGATTTCATAGGAATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTtcatc                                                                                       <  1:262250/148‑1 (MQ=255)
                                                            ccAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAg                                                                                     >  1:3145481/1‑150 (MQ=255)
                                                              aGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGt                                                                                   >  1:3690739/1‑150 (MQ=255)
                                                               taTGCAACTGCGTTACAGCGTGTCGTTCACGCTCAGTAAAAGCTATTTCATCGGAATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTc                                                                                  <  1:3618307/149‑1 (MQ=255)
                                                               gATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGTAAAAGCTATTTCATAGGCATCACCAAAGCGGATGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTc                                                                                  <  1:3618334/150‑1 (MQ=255)
                                                                aTGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTGCATCAGGTCg                                                                                 >  1:1327468/1‑150 (MQ=255)
                                                                  gCAACTGCGGTACTGCGTGTCTTTCACGCTCCGGAAAAGCGATTTCTCAGCAATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGc                                                                               <  1:2201974/150‑1 (MQ=255)
                                                                    aaCTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGt                                                                             >  1:705517/1‑150 (MQ=255)
                                                                           gTACAGCGTGTCGTTCACGCTCCGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgc                                                                      <  2:640735/150‑1 (MQ=255)
                                                                           gTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgc                                                                      <  2:3774744/150‑1 (MQ=255)
                                                                            tACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTgcgcg                                                                     >  2:3533986/1‑150 (MQ=255)
                                                                                      cGTCCACGCTCAGGAAAAGCGAGTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTAttttt                                                           <  1:5016717/150‑1 (MQ=255)
                                                                                       gTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTc                                                          >  1:4913564/1‑150 (MQ=255)
                                                                                       gTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTc                                                          <  2:4863978/150‑1 (MQ=255)
                                                                                        ttCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcg                                                         >  1:3798405/1‑150 (MQ=255)
                                                                                         tCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTcgc                                                        <  1:2415050/150‑1 (MQ=255)
                                                                                              cTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTAcc                                                   <  1:2880551/150‑1 (MQ=255)
                                                                                              cTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTAcc                                                   <  1:322244/150‑1 (MQ=255)
                                                                                               tCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCg                                                  >  1:1016712/1‑150 (MQ=255)
                                                                                                 aGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGaa                                                >  2:5103367/1‑150 (MQ=255)
                                                                                                   gAAAAGCGATTTCATAGGCATCACCAAGGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGc                                              <  2:3490997/150‑1 (MQ=255)
                                                                                                   gAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTgg                                                                                                                                >  1:4225737/1‑68 (MQ=255)
                                                                                                    aaaaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGcc                                             <  2:4351024/150‑1 (MQ=255)
                                                                                                      aaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGt                                           >  2:128808/1‑150 (MQ=255)
                                                                                                      aaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGt                                           >  1:1712290/1‑150 (MQ=255)
                                                                                                      aaGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGt                                           >  2:2618355/1‑150 (MQ=255)
                                                                                                       aGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTg                                          >  1:4569804/1‑150 (MQ=255)
                                                                                                          gATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcg                                       >  2:4952213/1‑150 (MQ=255)
                                                                                                          gATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcg                                       >  1:1605780/1‑150 (MQ=255)
                                                                                                           aTTTCATCGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgt                                      <  1:682262/150‑1 (MQ=255)
                                                                                                           aTTTCATAGGCATCACCAACGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCTTCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCAGTGGcgt                                      <  2:3643434/150‑1 (MQ=255)
                                                                                                            tttCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtc                                     >  1:2320796/1‑150 (MQ=255)
                                                                                                            tttCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtc                                     >  1:4891603/1‑150 (MQ=255)
                                                                                                            tttCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtc                                     <  2:3907802/150‑1 (MQ=255)
                                                                                                              tCATAGGCATCACAAAAGCGGCTGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGcgtcgt                                   <  1:2315972/150‑1 (MQ=255)
                                                                                                               cATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTa                                  <  1:2488010/150‑1 (MQ=255)
                                                                                                                aTAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAg                                 >  1:1291821/1‑150 (MQ=255)
                                                                                                                aTAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAg                                 >  1:4449076/1‑150 (MQ=255)
                                                                                                                   ggCATCACCAACGCGGCCGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAAc                              <  1:1433825/150‑1 (MQ=255)
                                                                                                                    gCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACg                             >  2:1641880/1‑150 (MQ=255)
                                                                                                                    gCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACg                             >  2:58957/1‑150 (MQ=255)
                                                                                                                         aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttc                        <  2:1605780/150‑1 (MQ=255)
                                                                                                                         aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttc                        <  1:2935889/150‑1 (MQ=255)
                                                                                                                         aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttc                        >  2:3733934/1‑150 (MQ=255)
                                                                                                                         aCCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttc                        <  1:2936048/150‑1 (MQ=255)
                                                                                                                           cAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                      <  1:3875029/150‑1 (MQ=255)
                                                                                                                           cAAAGCGGCAGGGCGGCGGGCGCCCCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGcttctt                      <  1:1467804/150‑1 (MQ=255)
                                                                                                                            aaaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTg                     <  2:4303953/150‑1 (MQ=255)
                                                                                                                             aaGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGt                    <  2:3207209/150‑1 (MQ=255)
                                                                                                                              ctcGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTc                   >  1:3903155/3‑150 (MQ=255)
                                                                                                                              aGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTc                   <  2:1503638/150‑1 (MQ=255)
                                                                                                                               tcGGCAGTGCGGCGCCCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTGCGTGTCg                  >  1:4986965/2‑150 (MQ=255)
                                                                                                                                  gCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACg               <  2:143976/150‑1 (MQ=255)
                                                                                                                                  gCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACg               >  1:116601/1‑150 (MQ=255)
                                                                                                                                   cAGTGCGGCGCCCGCCGCCCTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACgg              >  1:3928349/1‑150 (MQ=255)
                                                                                                                                    aGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGt             <  1:1705502/150‑1 (MQ=255)
                                                                                                                                     gggCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTa            <  1:1786295/150‑1 (MQ=255)
                                                                                                                                          gcgggcgCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTTATGaa       <  2:4209434/150‑1 (MQ=255)
                                                                                                                                               cgccgcACTGCCGAGGTCGATTAGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGa  <  1:2583865/150‑1 (MQ=255)
                                                                                                                                                gccgcACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGct                          <  2:2844217/125‑1 (MQ=255)
                                                                                                                                                gccgcACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGct                          >  1:2844217/1‑125 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
GCGGGCGAAGATGCCCCACGGGATGATCAGCGTCATGGCCAGGATCGACAGGCCCGCGATCCAGATGCAACTGCGGTACAGCGTGTCGTTCACGCTCAGGAAAAGCGATTTCATAGGCATCACCAAAGCGGCAGGGCGGCGGGCGCCGCACTGCCGAGGTCGATTGGGCAAGGGCGGCTGTTACTGGACGGCGTCGATACGCTTCATCAGGTCGGCGTACTGCGCGCCGTATTTTTCGCGTACCGAAGCCGTGGCGTCGTAGAACGGCTTCTTGTCGACGGTAATGAACTCGA  >  NC_002947/1343972‑1344264

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: ‑