Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,921,591 T→A 11.2% intergenic (+24/+37) 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,5910TA11.2% 269.1 / 19.4 96intergenic (+24/+37)PP_1720/PP_1721alcohol dehydrogenase/phosphoserine phosphatase
Reads supporting (aligned to +/- strand):  ref base T (41/44);  new base A (3/8);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTC  >  NC_002947/1921446‑1921739
                                                                                                                                                 |                                                                                                                                                    
ggTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAcc                                                                                                                                                  <  2:2929618/150‑1 (MQ=255)
 gTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAccc                                                                                                                                                 <  2:270043/150‑1 (MQ=255)
  ttGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCcg                                                                                                                                                <  1:4600094/150‑1 (MQ=255)
    gTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCcgcg                                                                                                                                              <  2:2882589/150‑1 (MQ=255)
     tCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGa                                                                                                                                             <  1:661768/150‑1 (MQ=255)
     tCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGa                                                                                                                                             <  1:661634/150‑1 (MQ=255)
      cgccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCTGTGCCGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                            <  2:211910/150‑1 (MQ=255)
      cgccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGaa                                                                                                                                            <  1:5311848/150‑1 (MQ=255)
      cgccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTAAACCCGCGaa                                                                                                                                            <  2:1429237/150‑1 (MQ=255)
       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAg                                                                                                                                           >  2:2568497/1‑150 (MQ=255)
       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAg                                                                                                                                           <  1:576356/150‑1 (MQ=255)
       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCTGGTAAACCCGCGAAg                                                                                                                                           <  2:2101261/150‑1 (MQ=255)
       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCCGCGAAg                                                                                                                                           <  2:87782/150‑1 (MQ=255)
       gccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCCGCGAAg                                                                                                                                           <  2:356302/150‑1 (MQ=255)
        ccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgg                                                                                                                                          <  2:1375277/150‑1 (MQ=255)
        ccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgg                                                                                                                                          <  2:1375262/150‑1 (MQ=255)
        ccgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCCGCGAAgg                                                                                                                                          <  1:5093621/150‑1 (MQ=255)
         cgccgcGGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAggg                                                                                                                                         <  1:1927294/150‑1 (MQ=255)
         cgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAggg                                                                                                                                         >  2:144335/1‑150 (MQ=255)
         cgcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCCGGCCCCTTCGCGGGTAAACCCGCGAAggg                                                                                                                                         <  1:460950/150‑1 (MQ=255)
          gcAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAgggg                                                                                                                                        <  1:2332260/150‑1 (MQ=255)
           cAGCTGGTCGACACCTCTCGGGTGGTGTTTGCCCCGGCGGCGTATGCCGTGCTGGAGACCAACCTGGTGTCGGGCATGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGGGGGGCTGGCCCCTTGGCGGGTTTACCCGCGAAGGGGc                                                                                                                                       <  2:2975247/150‑1 (MQ=255)
               tGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGc                                                                                                                                   >  1:4517583/1‑150 (MQ=255)
                    gACACCTATCGGGTGTATTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCTGGTGTCGGGCAAGCTGGTTTTGGTGATCGACCCTAGCCTGGTTTAAGCGGGGCTGGCCCCTTCGCGGTTTTACCCGCGAAGGGGCCGGCACAGc                                                                                                                              <  2:4204326/150‑1 (MQ=255)
                      cacCTATCGGGTGGAGTTTGCCCAGGCGGCTTTTGCCGAGCTGGAGACCAACCTGGTGTCGGGCACGGTGGTGTTGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTGGGGGGTTTACCCGCGAAGGGGCCGGCACAGccc                                                                                                                            <  2:4828457/150‑1 (MQ=255)
                       cccTATCCGGTGGATTTTGCCCTGGCGGCTTATGCCGCGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGAGGGTGATCGCCCCTACCCTGGTGTACGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCa                                                                                                                           <  2:3874228/149‑1 (MQ=255)
                       acCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTTCCCGCGAAGGGGCCGGCACACCCCa                                                                                                                           >  1:1955321/1‑150 (MQ=255)
                       acCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCa                                                                                                                           >  2:2086967/1‑150 (MQ=255)
                         cTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCaga                                                                                                                         >  1:2451966/1‑148 (MQ=255)
                             ccGGTGGAGTTGGCCCCGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTACGGGGGGCTGGCCCCTTCGGGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATaaaaa                                                                                                                     <  2:1636723/150‑1 (MQ=255)
                                 tgttttttGCCCTGGCGGCTTAGGCGGAGCTGGAGACCACCCCGGTGTCGGGCATGGGGGTGTTGGTGATCGACCCTCGCCTGTTTTACGGGGGGGTGGCCCCTTGGGGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 <  1:1749473/145‑1 (MQ=255)
                                 tGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTTATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 <  2:4943572/150‑1 (MQ=255)
                                 tGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 <  2:831423/150‑1 (MQ=255)
                                 tGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACAAACCAGGTGTGGGGCAAGCTTGTGATGTTGATCGACCCTCGCCTGGTTTAAGGGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCtt                                                                                                                 <  2:831521/150‑1 (MQ=255)
                                  ggtttttGCCCTGGCGGCGTATGCCGCGCTGGAGACCAACCTGGTGTGGGGCATGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGGGGGGCTGGCCCCTTCGCGGTTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                                <  2:2898383/146‑1 (MQ=255)
                                  ggAGTTTGCCCAGGCGGCGTATGCCTAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAGGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTc                                                                                                                <  2:378069/150‑1 (MQ=255)
                                   gAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTca                                                                                                               >  1:2692890/1‑150 (MQ=255)
                                       ttGCCCCGGCGGGGTATGCCGCGCTGGAGACAAACCAGGGGTCGGGCATGGTGGTGAGGGTTATCGACCCTAGCGTGGTGTACGGGGGGCGGGCCCCTTCGCGGTTTTACCGGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctc                                                                                                           <  2:4598859/150‑1 (MQ=255)
                                          cccTGGCGGCGTATGCGGCGCTGGATACCACCCCGGTGTCGGGCATGCTGGTGATGTTGATGGACCCTTGCCTGGTTTACGCGGGGCTGGCCCCTTCGGGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccag                                                                                                        <  1:4082601/150‑1 (MQ=255)
                                           ccAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTcagctccagc                                                                                                       >  2:4622241/1‑150 (MQ=255)
                                                      ttGCCGAGCTGGAGACCAACCTGGTGTGGGGAAAGCTGGTGTTGGTGATCGACCCTAGCGTGGTTTAAGGGGGGCTGGCCCCTTCGCGGTTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGc                                                                                            <  2:626989/149‑1 (MQ=255)
                                                       tGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATACAAACCTTCAGCTCCAGCTCAGAGTCGGcc                                                                                           >  2:5253854/1‑150 (MQ=255)
                                                        gCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGc                                                                                                                              <  2:451383/114‑1 (MQ=255)
                                                        gCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGc                                                                                                                              >  1:451383/1‑114 (MQ=255)
                                                          cGAGCTGGAGACCAACCTGGTGTCGGGCAAGCTGGTGATGGTTATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAcc                                                                                        <  2:4653940/150‑1 (MQ=255)
                                                          cGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAcc                                                                                        >  2:3961839/1‑150 (MQ=255)
                                                          cGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAcc                                                                                        >  2:2316353/1‑150 (MQ=255)
                                                          cGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGACCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCAcc                                                                                        >  2:2822236/1‑150 (MQ=255)
                                                            tgCTGGAGACCACCCTGGTGTGGGGCAAGCTGGTTCTGGTGATCGACCCTAGCCTGGTGTAAGGGGGGCTGCCCCCTTCGCGGTTTAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCa                                                                                      <  1:810900/149‑1 (MQ=255)
                                                             gCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGTTGGTGATCGACCCTAGCCGGGTGTAAGCGGGGCTGGCCCCTTCGGGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAt                                                                                     <  1:4414822/150‑1 (MQ=255)
                                                                ggAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGa                                                                                  >  1:4871291/1‑150 (MQ=255)
                                                                   gaCCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGAtct                                                                               >  2:2471062/1‑150 (MQ=255)
                                                                    accaaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAACCGCGAAGGGGCCGGCCCAGCCCCAAAAAAACATCAGCTCCAGCACAGAGTAGGCCACCCAATGAtatc                                                                              >  1:244203/1‑150 (MQ=255)
                                                                    accaaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAACCCCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCGCAGAGTCGGCCACCCAATGAtcac                                                                              >  1:244251/1‑148 (MQ=255)
                                                                    accaaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCGTCGCGGGTTTACCCGCGCAGGGGCCGGCACAGCCCATAAAAACCATCAGCTCCAGCGCACTGTCGGCCACCCATAGAtatc                                                                              >  1:243982/1‑150 (MQ=255)
                                                                     ccaaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGTTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCg                                                                             <  2:299575/150‑1 (MQ=255)
                                                                      caaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCCTAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAATGATCTCgg                                                                            >  1:196756/1‑150 (MQ=255)
                                                                      caaccaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAATGATCACgg                                                                            >  1:4337540/1‑150 (MQ=255)
                                                                        accaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCa                                                                          >  1:4870286/1‑150 (MQ=255)
                                                                        accaGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTAAGCTCCAGCTCAGAGTAGGCCACCCATTGTTATCGGCa                                                                          >  1:1115670/1‑150 (MQ=255)
                                                                               gtCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGTTCTGGGCATGCtagg                                                                   >  2:1352283/1‑147 (MQ=255)
                                                                               ggcGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATACAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTcgc                                                                   <  2:2860560/148‑1 (MQ=255)
                                                                                  gggCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAAACTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCa                                                                >  1:4828951/1‑150 (MQ=255)
                                                                                        gCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTg                                                          >  1:4843190/1‑150 (MQ=255)
                                                                                         cTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGc                                                         >  1:3181245/1‑150 (MQ=255)
                                                                                             tgatggttatCCACCCTAGCCTGGTGTAAGCGGGGCTGGCCCATTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTcc                                                     <  2:872008/150‑1 (MQ=255)
                                                                                             tgatggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTcc                                                     <  2:3433225/150‑1 (MQ=255)
                                                                                             tgatggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTcc                                                     <  2:3613645/150‑1 (MQ=255)
                                                                                               ttggtgttCGACCTTTGCGTGGTGTTTGCGGGGCTGGGCCGTTGGGGGGTTTTGCGGCGAAGGGGCCGGGACAGCCCTTTAAAAGGTTGTGTTGGAGCTCAGAGTCGGCCAGCCTTTGATCTCGGCATGCTCGCGCAACACTGCGTCCgg                                                   <  2:3628037/149‑1 (MQ=255)
                                                                                               atggtgatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTAACCGCGAAGGGGCCGGCACAGCCCATAAAAAACTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTAGGCATTCTCGCGCAACACTGCGTCCgg                                                   >  1:4778766/1‑150 (MQ=255)
                                                                                                    gatCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTCCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTca                                              >  2:2542756/1‑150 (MQ=255)
                                                                                                      tCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacc                                            <  2:2979075/150‑1 (MQ=255)
                                                                                                       cGACCCTAGCCTGGTGTAAGCGGGGCGGGCCCCTTCGCGGGTTTCCCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTcacca                                           <  2:4607083/150‑1 (MQ=255)
                                                                                                             tAGCCTCGTGTAAGCGGGGCTGGCCCCTACGCGGGTTAACCCGCGAATCGGCCGGCACAGCCCATAAAAACATGAATCACCAGATAACAGTCCGCCAACCAATGAACAAGGCAAGCAAGCGCAACACATCGACCCGGTTCACCACATGCg                                     >  2:457261/1‑150 (MQ=255)
                                                                                                                ccTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCCACACTGCGTCCGGGTTCACCACATGCGGGt                                  >  1:3037485/1‑150 (MQ=255)
                                                                                                                 cTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTACCCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTg                                 <  2:2930465/150‑1 (MQ=255)
                                                                                                                    gtgtAAGCGGGGCTGGCCCCTTCGCGGGTTTCCCCGCGACGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATc                              <  2:203651/150‑1 (MQ=255)
                                                                                                                       tAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCCACACTGCGTCCGGGTTCCCCACATGCGGGTGATCCAc                           >  2:1509396/1‑150 (MQ=255)
                                                                                                                        aaGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCAct                          <  2:3017867/150‑1 (MQ=255)
                                                                                                                          gCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCGTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCActct                        >  2:4220971/1‑150 (MQ=255)
                                                                                                                              ggCGGGCCCCTTCTGGGGTTTCCCCGCGATGGTGCCGGCAGAGCCCATCTAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGTCATGCTCGCGCAAGAATGCGTCCGGGTTGACCACTTGGGGGTGATCCACTCTcaac                    <  2:349884/150‑1 (MQ=255)
                                                                                                                                  ggCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaacaaca                >  2:1665250/1‑150 (MQ=255)
                                                                                                                                  ggCCCCTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTcaacaaca                >  2:3717194/1‑150 (MQ=255)
                                                                                                                                     cccTTCGGGGTTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg             <  2:686077/150‑1 (MQ=255)
                                                                                                                                     cccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACATCCAATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTACTTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg             >  2:3322456/1‑150 (MQ=255)
                                                                                                                                     cccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACg             >  2:4122903/1‑150 (MQ=255)
                                                                                                                                     cccTTCGCGGGTATACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCAATGATCGCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGAGATCCACTCTCAACAACAACg             >  2:1332275/1‑150 (MQ=255)
                                                                                                                                      ccTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACgg            >  2:4798865/1‑150 (MQ=255)
                                                                                                                                       cTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGc           <  2:3801661/150‑1 (MQ=255)
                                                                                                                                       cTTCGCGGGTAAACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGc           >  2:4557559/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa          >  1:4777640/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa          >  1:4777483/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa          >  2:4228417/1‑150 (MQ=255)
                                                                                                                                        ttCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCa          >  1:5087023/1‑150 (MQ=255)
                                                                                                                                           gcgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGt       >  1:638287/1‑150 (MQ=255)
                                                                                                                                           gcgGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGt       <  1:2086967/150‑1 (MQ=255)
                                                                                                                                               ttttACCCGCTATGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCAtt   <  1:2568497/149‑1 (MQ=255)
                                                                                                                                                tttACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTc  >  1:561954/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GGTTGTCGCCGCAGCTGGTCGACACCTATCCGGTGGAGTTTGCCCAGGCGGCGTATGCCGAGCTGGAGACCAACCAGGTGTCGGGCAAGCTGGTGATGGTGATCGACCCTAGCCTGGTGTAAGCGGGGCTGGCCCCTTCGCGGGTTTACCCGCGAAGGGGCCGGCACAGCCCATAAAAACCTTCAGCTCCAGCTCAGAGTCGGCCACCCATTGATCTCGGCATGCTCGCGCAACACTGCGTCCGGGTTCACCACATGCGGGTGATCCACTCTCAACAACAACGGCAAGTCATTC  >  NC_002947/1921446‑1921739

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