Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 3,961,904 C→T T416T (ACC→ACT RPPX_17635 → damage‑inducible protein CinA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099743,961,9040CT99.6% 977.5 / NA 265T416T (ACC→ACTRPPX_17635damage‑inducible protein CinA
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (113/151);  total (113/152)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

CCCGGTGGAAATCGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGC  >  CP009974/3961757‑3962050
                                                                                                                                                   |                                                                                                                                                  
cccGGTGGAAATCGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACttct                                                                                                                                                 <  2:4184965/151‑1 (MQ=255)
  cGGTGGAAATCGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACttctct                                                                                                                                               >  1:2306924/1‑151 (MQ=255)
    gTGGAAATCGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTgg                                                                                                                                             <  2:3559145/151‑1 (MQ=255)
    gTGGAAATCGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTgg                                                                                                                                             <  1:766797/151‑1 (MQ=255)
         aaTCGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATc                                                                                                                                        <  1:247634/151‑1 (MQ=255)
          aTCGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCa                                                                                                                                       <  1:4427702/151‑1 (MQ=255)
            cGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCAcg                                                                                                                                     >  1:3030158/1‑150 (MQ=255)
            cGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACt                                                                                                                                     <  1:2750015/151‑1 (MQ=255)
            cGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACt                                                                                                                                     <  1:1387319/151‑1 (MQ=255)
             gCCCGCGCCTCGGTGCTGGGGCAGCCTTTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTg                                                                                                                                    >  1:285937/1‑151 (MQ=255)
             gCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTg                                                                                                                                    <  1:4025854/151‑1 (MQ=255)
              cccGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTgc                                                                                                                                   <  2:4301256/151‑1 (MQ=255)
              cccGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTgc                                                                                                                                   <  2:4550408/151‑1 (MQ=255)
              cccGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTgc                                                                                                                                   >  1:2584488/1‑151 (MQ=255)
               ccGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTgcg                                                                                                                                  >  1:1071168/1‑151 (MQ=255)
               ccGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTgcg                                                                                                                                  <  1:842771/151‑1 (MQ=255)
               ccGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTgcg                                                                                                                                  <  1:5140147/151‑1 (MQ=255)
                cgcgCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTgcgc                                                                                                                                 >  1:5181701/1‑151 (MQ=255)
                cgcgCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTgcgc                                                                                                                                 >  1:4507806/1‑151 (MQ=255)
                   gcCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTgg                                                                                                                              >  1:1111342/1‑151 (MQ=255)
                   gcCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTgg                                                                                                                              >  2:5002506/1‑151 (MQ=255)
                      tCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCAc                                                                                                                           >  1:2658077/1‑151 (MQ=255)
                      tCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCAc                                                                                                                           <  2:3734285/151‑1 (MQ=255)
                        ggTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCActg                                                                                                                         <  2:607272/151‑1 (MQ=255)
                        ggTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCActg                                                                                                                         >  2:4083381/1‑151 (MQ=255)
                          tGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCActgct                                                                                                                       <  1:4736091/151‑1 (MQ=255)
                             tGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGAc                                                                                                                    <  1:269535/151‑1 (MQ=255)
                             tGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGAc                                                                                                                    <  1:2825966/151‑1 (MQ=255)
                             tGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGAc                                                                                                                    <  1:269612/151‑1 (MQ=255)
                              ggggCAGCCTGGGCCCCGGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACg                                                                                                                   <  2:4682074/151‑1 (MQ=255)
                                 gCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGCCGGgc                                                                                                                >  2:2864367/1‑151 (MQ=255)
                                 gCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGgc                                                                                                                >  1:874600/1‑151 (MQ=255)
                                  cAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTGCACTGGTGATCACTGCGCTGGCACTGCTGACGGgcg                                                                                                               >  2:3022474/1‑151 (MQ=255)
                                     ccTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCt                                                                                                            >  2:2702317/1‑151 (MQ=255)
                                     ccTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCt                                                                                                            <  2:1299479/151‑1 (MQ=255)
                                      cTGTGCACCTGATGAGCCCGCTGGTGGCTTCCTCCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCtt                                                                                                           <  1:524487/151‑1 (MQ=255)
                                           cACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCt                                                                                                      <  1:1636118/151‑1 (MQ=255)
                                                gatgaGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATg                                                                                                 <  2:3333483/151‑1 (MQ=255)
                                                gatgaGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATg                                                                                                 <  2:4172227/151‑1 (MQ=255)
                                                gatgaGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATg                                                                                                 >  1:6902/1‑151 (MQ=255)
                                                 atgaGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATgg                                                                                                <  1:1670180/151‑1 (MQ=255)
                                                     gCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTATGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTg                                                                                            >  2:4615779/1‑151 (MQ=255)
                                                       ccGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGaa                                                                                          >  2:816830/1‑151 (MQ=255)
                                                       ccGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGaa                                                                                          >  1:58249/1‑151 (MQ=255)
                                                        ggCTGGTGGGTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAg                                                                                         <  2:4499965/150‑1 (MQ=255)
                                                        cGCTGGTGGCTTCCGCCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAg                                                                                         <  1:3283165/151‑1 (MQ=255)
                                                        cGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAg                                                                                         >  2:5233563/1‑151 (MQ=255)
                                                           tggtggCTTCCGCCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCAc                                                                                      <  1:3845298/151‑1 (MQ=255)
                                                           tggtggCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCAc                                                                                      >  2:3921556/1‑151 (MQ=255)
                                                             gtggCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTc                                                                                    >  1:1430388/1‑151 (MQ=255)
                                                             gtggCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTc                                                                                    <  1:3748468/151‑1 (MQ=255)
                                                               ggCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAg                                                                                  >  1:1550089/1‑151 (MQ=255)
                                                                 cTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGGTCACTGCGCTGGCACTGGTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAggg                                                                                <  2:3614809/151‑1 (MQ=255)
                                                                   tCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGgcg                                                                              <  2:3494363/151‑1 (MQ=25)
                                                                      aCCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGcgccg                                                                           <  1:3495176/151‑1 (MQ=21)
                                                                         ttgttgCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTc                                                                        <  1:4838530/151‑1 (MQ=18)
                                                                              gCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGgcgc                                                                   >  1:3179139/1‑151 (MQ=18)
                                                                              gCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGgcgc                                                                   <  1:181431/151‑1 (MQ=18)
                                                                               cTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCCCGGCGCCGCTCGGCGcc                                                                  >  2:5145490/1‑151 (MQ=255)
                                                                                 ggtggGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGcccc                                                                <  1:19308/151‑1 (MQ=25)
                                                                                  gtggGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCg                                                                       <  1:2328215/143‑1 (MQ=14)
                                                                                      gCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGACCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCACCTCTGa                                                           <  2:2332345/151‑1 (MQ=255)
                                                                                        aTGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACt                                                         >  1:2147725/1‑151 (MQ=255)
                                                                                         tGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCCTGGTTTGTAGCACTCAGGGCGCCCCTCGGCGCCCCTCTGACTc                                                        >  2:4431536/1‑151 (MQ=255)
                                                                                          ggTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCa                                                       >  1:4074886/1‑151 (MQ=255)
                                                                                           gTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAGGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCAc                                                      <  2:4005111/151‑1 (MQ=255)
                                                                                           gTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCAc                                                      >  1:4262260/1‑151 (MQ=255)
                                                                                             cGACAGGGATAGCGGTGAGTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCAccc                                                    <  2:1516363/151‑1 (MQ=255)
                                                                                              gACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCa                                                   >  2:1080149/1‑151 (MQ=255)
                                                                                              gACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCa                                                   >  2:1085239/1‑151 (MQ=255)
                                                                                               aCAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCaa                                                  >  1:1279939/1‑151 (MQ=255)
                                                                                                cAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAg                                                 <  1:5002506/151‑1 (MQ=255)
                                                                                                cAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAg                                                 >  2:2083140/1‑151 (MQ=255)
                                                                                                 aGGGATATCGGTGATTTCCAGAAAGCCGCCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGa                                                <  2:4820647/151‑1 (MQ=255)
                                                                                                 aGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGa                                                >  2:3031502/1‑151 (MQ=255)
                                                                                                 aGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGa                                                >  2:2403102/1‑151 (MQ=255)
                                                                                                 aGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGa                                                <  1:2080243/151‑1 (MQ=255)
                                                                                                  gggATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGaa                                               <  1:2978110/151‑1 (MQ=255)
                                                                                                      tatCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGaacaac                                           <  2:3139721/151‑1 (MQ=255)
                                                                                                        tCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGaacaacaa                                         <  1:2794671/151‑1 (MQ=255)
                                                                                                         cGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGaaaaacaac                                        >  2:2112182/1‑151 (MQ=255)
                                                                                                                ttCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTaa                                 >  1:5109088/1‑151 (MQ=255)
                                                                                                                 tCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAAc                                >  1:3405695/1‑151 (MQ=255)
                                                                                                                     gAAAGCCACGTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACcatc                            <  2:4932374/151‑1 (MQ=255)
                                                                                                                     gAAAGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACcatc                            <  2:169463/151‑1 (MQ=255)
                                                                                                                       aaGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACcatcat                          <  1:3535000/151‑1 (MQ=255)
                                                                                                                        aGCCACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGCGCAACCATCATg                         >  2:1737603/1‑151 (MQ=255)
                                                                                                                          ccACCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTc                       >  1:4779546/1‑151 (MQ=255)
                                                                                                                            cccTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTccc                     <  2:2607292/150‑1 (MQ=255)
                                                                                                                            aCCTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTccc                     <  2:3858683/151‑1 (MQ=255)
                                                                                                                             ccgtCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCg                    <  1:3954792/148‑1 (MQ=255)
                                                                                                                             ccTTCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCg                    <  2:2275602/151‑1 (MQ=255)
                                                                                                                               ttCAAGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGAt                  <  2:477455/151‑1 (MQ=255)
                                                                                                                                   aGTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCACCATGTCCGGATTTCt              >  2:3204385/1‑151 (MQ=255)
                                                                                                                                    gTGGGCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCtt             <  2:1653234/151‑1 (MQ=255)
                                                                                                                                       ggCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTc          <  1:2032104/151‑1 (MQ=255)
                                                                                                                                        gCCGTGCTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTcc         >  2:656053/1‑151 (MQ=255)
                                                                                                                                              cTCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGTGCCTTGTCCTTCATGGTTTGAAGCACTCTGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAAAAACAGGGTAAGCATCATGTACCGGTCTCTTCTCCGTgggg   >  1:1952617/1‑151 (MQ=21)
                                                                                                                                               tCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGc  >  1:4294881/1‑151 (MQ=255)
                                                                                                                                               tCACTTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGc  <  2:2497116/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCCGGTGGAAATCGCCCGCGCCTCGGTGCTGGGGCAGCCTGTGCACCTGATGAGCCCGCTGGTGGCTTCCACCTTGTTGCTGGTGGGCATGGTCGACAGGGATATCGGTGATTTCCAGAAAGCCACCTTCAAGTGGGCCGTGCTCACCTCTCTGGTGATCACTGCGCTGGCACTGCTGACGGGCGCCTTGTCCTTCATGGTTTGAAGCACTCAGGGCGCCGCTCGGCGCCCCTCTGACTCACCCAAGAACAACAACAGAGTAACCATCATGTCCCGATTTCTTTTCCGTGGGGC  >  CP009974/3961757‑3962050

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