Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,921,592 T→A 20.1% intergenic (+25/+36) PP_1720 → / ← PP_1721 alcohol dehydrogenase/phosphoserine phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,921,5920TA20.1% 410.6 / 93.3 196intergenic (+25/+36)PP_1720/PP_1721alcohol dehydrogenase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (98/57);  new base A (15/22);  total (113/83)
Fisher's exact test for biased strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 202 total aligned reads displayed.

GCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCT  >  NC_002947/1921444‑1921741
                                                                                                                                                    |                                                                                                                                                     
gCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTa                                                                                                                                                      <  2:1005311/150‑1 (MQ=255)
   gTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAccc                                                                                                                                                   <  1:2863100/150‑1 (MQ=255)
   gTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAAccc                                                                                                                                                   <  1:6955839/150‑1 (MQ=255)
      gTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgcg                                                                                                                                                >  2:7471376/1‑150 (MQ=255)
      gTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcgcg                                                                                                                                                >  2:3552171/1‑150 (MQ=255)
      gTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCcgcg                                                                                                                                                <  2:7852556/150‑1 (MQ=255)
       tCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGa                                                                                                                                               <  2:5811064/150‑1 (MQ=255)
        cgccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGaa                                                                                                                                              <  2:8717878/150‑1 (MQ=255)
        cgccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGaa                                                                                                                                              <  2:704830/150‑1 (MQ=255)
        cgccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                              <  2:667388/150‑1 (MQ=255)
        cgccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                              <  2:1757909/150‑1 (MQ=255)
         gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAg                                                                                                                                             <  2:6331853/150‑1 (MQ=255)
         gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGGTTACCCGCGAAg                                                                                                                                             >  1:1010448/1‑150 (MQ=255)
          ccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAgg                                                                                                                                            <  2:5635753/150‑1 (MQ=255)
          ccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAgg                                                                                                                                            <  1:5325858/150‑1 (MQ=255)
           cgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           <  1:4269273/150‑1 (MQ=255)
           cgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           <  1:4213506/150‑1 (MQ=255)
           cgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCCGCGAAggg                                                                                                                                           <  1:2723451/150‑1 (MQ=255)
            gcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgggg                                                                                                                                          >  2:8801736/1‑150 (MQ=255)
            gcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGACGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTATCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgggg                                                                                                                                          <  1:6799215/150‑1 (MQ=255)
               gCTGGTCGACACCTATCCGGGGGAGTTTGCCCAGGGGGCGTATGCCGTGCTGGAGACCAACCAGGTGTCGGGCACGCTGGTGATGGTGTTCGACCCTAGCCTTGTTTATGCGGGGGTGGCCCCTTCGCGGGTTTACCCGCTAAGGGGCCg                                                                                                                                       <  2:7168264/150‑1 (MQ=255)
                 tGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGAAg                                                                                                                                             >  1:5711061/1‑142 (MQ=255)
                   gTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGcac                                                                                                                                   >  2:6575856/1‑150 (MQ=255)
                      gACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGc                                                                                                                                >  1:1195749/1‑150 (MQ=255)
                        cacCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGccc                                                                                                                              >  2:3063431/1‑150 (MQ=255)
                          ccTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCAt                                                                                                                            >  1:8713763/1‑150 (MQ=255)
                           cTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATa                                                                                                                           >  2:8183427/1‑150 (MQ=255)
                           cTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATa                                                                                                                           >  1:6472056/1‑150 (MQ=255)
                           cTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATa                                                                                                                           >  2:708040/1‑150 (MQ=255)
                            tATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATaa                                                                                                                          >  1:709040/1‑150 (MQ=255)
                            tATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATaa                                                                                                                          >  1:3633977/1‑150 (MQ=255)
                             aTCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAc                                                                                                                                                     <  2:6141839/122‑1 (MQ=255)
                             aTCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAc                                                                                                                                                     >  1:6141839/1‑122 (MQ=255)
                               ccGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATaaaaa                                                                                                                       >  2:4015507/1‑150 (MQ=255)
                                   tGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                   <  2:2100609/150‑1 (MQ=255)
                                    ggttttggCCCTGGCGGCGTATGCCGAGCTGGAGACCACCCTGGTGCCGGGCACGCTGGTGTTGGTGATCGACCCTCGCGTGTTGTACGCGGGGCTGGCCCCTTCGGGGTTTTCCCCGCCACGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                                  <  1:2212778/143‑1 (MQ=255)
                                    ggtgTTTGCCCTGGGGGCGTTTGCCGAGCGGTAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGGGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                                  <  1:8255828/147‑1 (MQ=255)
                                             ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccagc                                                                                                         >  2:8737155/1‑150 (MQ=255)
                                             ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAACCGCGaa                                                                                                                                              >  1:4036300/1‑113 (MQ=255)
                                             ccAGGCGGCGGATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCGGGTGTAAGCGGGGCTGGCCCCTTCGCGGGGGTACCCGCGAAGGGGCCGGCACAGCCCAGAAAAACCGTcagctccagc                                                                                                         >  2:3986089/1‑150 (MQ=255)
                                               aGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGGTTACCCGCGAAGGGGCCGGCCCAGCCCCTAAAAACCTTcagctccagctc                                                                                                       >  2:8504254/1‑150 (MQ=255)
                                                  cggcgTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCaga                                                                                                    <  2:2892891/150‑1 (MQ=255)
                                                        aTGCCGAGCTGGAGACCAACCTGGTGTCGGGAATGGTGGTTATGGTTATCGACCCTAGCGTGGTGTAAGCGGGGCTGGCCCCTTGGGGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGc                                                                                              <  2:6392606/150‑1 (MQ=255)
                                                            cGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCCGCTCCAGCTCAGAGTCGGCCAcc                                                                                          >  2:544414/1‑150 (MQ=255)
                                                            cGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGGTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAcc                                                                                          >  1:6164698/1‑150 (MQ=255)
                                                             gAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCATCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCAAAAAAAAAATCAACTCCAGCACAGAGTAGGcccccc                                                                                         >  1:3351222/1‑150 (MQ=255)
                                                                   gagaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGAt                                                                                   <  2:2743688/150‑1 (MQ=255)
                                                                     gaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGAtct                                                                                 >  2:7321189/1‑150 (MQ=255)
                                                                       ccaacctGGTGTCGGGCACGCTGGTGATGGTGATCGACCCTAGCGTGGTGTAAGGGGGGCTGGCCCCTGCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCg                                                                               <  2:3208183/150‑1 (MQ=255)
                                                                       ccaacccGGTGTGGGGCACGCGGGTGATGGTGATCGACCCTAGCCTGGTGCACGCGGGGCGGGCCCCTTCGCGGGTTTACCCGCGACGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCg                                                                               <  1:3753533/150‑1 (MQ=255)
                                                                         aaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGc                                                                             >  2:6882318/1‑150 (MQ=255)
                                                                         aaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGc                                                                             >  2:243627/1‑150 (MQ=255)
                                                                          accaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCCTAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAATGATCTCGGCa                                                                            >  1:3909739/1‑150 (MQ=255)
                                                                                 gtCGGGCATGCTGTGGATGGGGCTCGACCCTCGCCTGGTGTACGGGGGGCTGCCCCCTTGGCGGGTTTACCGGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTcgc                                                                     <  1:6330223/150‑1 (MQ=255)
                                                                                 gtCGGGCACGCTGGTGATGGTGATCGACCCTAGCCTGGTGTACGGGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTcgc                                                                     <  2:1202745/150‑1 (MQ=255)
                                                                                     ggCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCaa                                                                 >  2:866068/1‑150 (MQ=255)
                                                                                      gCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCGCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAac                                                                >  1:4616168/1‑150 (MQ=255)
                                                                                          gCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTg                                                            >  2:4236459/1‑150 (MQ=255)
                                                                                           cTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGc                                                           >  1:7321530/1‑150 (MQ=255)
                                                                                           cTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGc                                                           >  2:3908810/1‑150 (MQ=255)
                                                                                              gtgatggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAATCCTTGAGCTCCAGGTCAGAGTCGGGCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTc                                                        >  2:7771686/1‑150 (MQ=255)
                                                                                                gatggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCg                                                      >  2:3957104/1‑150 (MQ=255)
                                                                                                 atggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCACAGAGTCGGCCACCCAATTATCTCGGCATGCTCGCGCAACACTGCGtaccg                                                     >  1:9111268/1‑146 (MQ=255)
                                                                                                     tgttCGCCCCTAGCGTGGTTTACGGGGGGCGGGCCCCTTCGGGGGTTCACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTc                                                 <  1:3466411/147‑1 (MQ=255)
                                                                                                      gatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTca                                                >  1:2959556/1‑150 (MQ=255)
                                                                                                      gatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAACCGCGAAGGGGCCGGGACAGCCCATAAAAAACTTGAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTca                                                >  1:1297504/1‑150 (MQ=255)
                                                                                                      gatCGACCCTAGCCTGGTGTAAGCGGGGCTGACCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTca                                                >  1:2784430/1‑150 (MQ=255)
                                                                                                       atCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCGTCGCGGGTTTACCCGCGAAGGGGCCGGCACACCCCAGAAAAACCGTCAGCTCCCGCTCAGAGTCGGCCACCCAAAGATCTCCGCATTCTCGCGCAACACAGCGTCCGGGTTcac                                               >  1:7801614/1‑150 (MQ=255)
                                                                                                       atCGACCCTAGCCTGGGGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcac                                               <  2:8779558/150‑1 (MQ=255)
                                                                                                           aCCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTg                                                                                     >  1:5890652/1‑108 (MQ=255)
                                                                                                            cccTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCAAAAAAACCATCAGCGCCAGCTCCGAGTCGGCCACCCATTGATCTCGGCAAGCTCGCGCAACACGGCGTCCGGGGTCACCACAt                                          >  1:1973749/1‑150 (MQ=255)
                                                                                                             ccTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCCTTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATg                                         >  1:676541/1‑150 (MQ=255)
                                                                                                             ccTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATg                                         >  1:7397828/1‑150 (MQ=255)
                                                                                                                tgcgtGGTGTATGGGGGGGTGGCCCGTTGGGGGTTTTACCGGCTACGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCgg                                      <  1:6575856/146‑1 (MQ=255)
                                                                                                                  ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGt                                    >  2:2170765/1‑150 (MQ=255)
                                                                                                                    tGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGa                                  <  2:5552808/150‑1 (MQ=255)
                                                                                                                      gtgtAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATc                                <  2:6816178/150‑1 (MQ=255)
                                                                                                                           aGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCActc                           <  2:4964381/150‑1 (MQ=255)
                                                                                                                                gggtGGCCCCTTGGGGTTTTTACCCGCTACGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaac                      <  1:2433757/147‑1 (MQ=255)
                                                                                                                                ggCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaac                      <  2:6909970/150‑1 (MQ=255)
                                                                                                                                     gCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                   <  2:5007834/52‑1 (MQ=255)
                                                                                                                                     gCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                   >  1:5007834/1‑52 (MQ=255)
                                                                                                                                       cccTTCGCGGGTAAACCCGCGAAGGGGCCAGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg               >  2:5840534/1‑150 (MQ=255)
                                                                                                                                        ccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg              >  1:2938265/1‑150 (MQ=255)
                                                                                                                                        ccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg              >  2:517485/1‑150 (MQ=255)
                                                                                                                                        ccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg              <  2:1307023/150‑1 (MQ=255)
                                                                                                                                        ccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg              <  1:5206053/150‑1 (MQ=255)
                                                                                                                                         cTTCGCGGGTAAACCCGCGAAGGGGCCAGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGc             >  1:8092398/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa            >  2:2539052/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa            >  1:8856215/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa            >  1:7892331/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa            >  2:7906205/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa            >  2:4307713/1‑150 (MQ=255)
                                                                                                                                          ttCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa            >  1:3276317/1‑150 (MQ=255)
                                                                                                                                           tCGCGGGTTTACCCGCTAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCTCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACCGCaa           <  2:5923078/150‑1 (MQ=255)
                                                                                                                                            cgcgGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCa                                                                                                      <  2:5235147/58‑1 (MQ=37)
                                                                                                                                             gcgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGACATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGt         >  1:1573459/1‑150 (MQ=255)
                                                                                                                                                  tttACCCGCTAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTc    <  2:447628/150‑1 (MQ=255)
                                                                                                                                                  tttACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTc    <  2:1195749/150‑1 (MQ=255)
                                                                                                                                                   ttACCCGCGATGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTcc   <  1:8802017/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTCCT  >  NC_002947/1921444‑1921741

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: