Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,096,752 G→T 8.6% S127S (TCG→TCT)  mutY → adenine DNA glycosylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,096,7520GT 8.6% 1267.2 / 62.3 350S127S (TCG→TCT) mutYadenine DNA glycosylase
Reads supporting (aligned to +/- strand):  ref base G (152/167);  new base T (14/17);  total (166/184)
Fisher's exact test for biased strand distribution p-value = 8.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 357 total aligned reads displayed.

ACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGT  >  CP009273/3096603‑3096901
                                                                                                                                                     |                                                                                                                                                     
   gggCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCtt                                                                                                                                                    >  2:1343419/1‑150 (MQ=255)
    ggCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACCAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCttt                                                                                                                                                   <  2:498996/150‑1 (MQ=255)
    ggCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCttt                                                                                                                                                   <  2:3489248/150‑1 (MQ=255)
      cTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTtct                                                                                                                                                 <  1:1676703/150‑1 (MQ=255)
           cTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTggg                                                                                                                                            >  1:4893591/1‑150 (MQ=255)
             aTTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTa                                                                                                                                          >  1:5037072/1‑150 (MQ=255)
             aTTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTa                                                                                                                                          <  1:2679909/150‑1 (MQ=255)
              ttACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCTCTTTCTCTGGGTaa                                                                                                                                         <  1:2049046/150‑1 (MQ=255)
                  gCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGACGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCAc                                                                                                                                     >  2:2976744/1‑150 (MQ=255)
                       cgcgcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCTCTTTCTCTGGGTAAGCACTTTcc                                                                                                                                >  2:3657427/1‑150 (MQ=255)
                       cgcgcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTcc                                                                                                                                >  1:3869535/1‑150 (MQ=255)
                       cgcgcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTcc                                                                                                                                <  1:4156244/150‑1 (MQ=255)
                       cgcgcgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTcc                                                                                                                                <  1:3029414/150‑1 (MQ=255)
                           cgcAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTATCTCGCTTTCTCTGGGTAAGCACTTTCCGAtt                                                                                                                            <  1:1003613/150‑1 (MQ=255)
                             cAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCTCTTTCTCTGGGTAAGCACTTTCCGATtct                                                                                                                          >  1:4843145/1‑150 (MQ=255)
                              aaTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATtctc                                                                                                                         >  1:1222438/1‑150 (MQ=255)
                              aaTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTCCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATtctc                                                                                                                         >  1:4026168/1‑150 (MQ=255)
                              aaTCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCGTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATtctc                                                                                                                         >  2:2227606/1‑150 (MQ=255)
                                 cTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGAc                                                                                                                      >  2:3835751/1‑150 (MQ=255)
                                     aTAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTa                                                                                                                  >  2:3088070/1‑150 (MQ=255)
                                     aTAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTa                                                                                                                  >  2:374340/1‑150 (MQ=255)
                                      tAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTaa                                                                                                                 >  1:2872653/1‑150 (MQ=255)
                                         aGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCTCTTTCTCTGGGTAAGCACTTTCCGATTATCGACGGTAACGt                                                                                                              <  2:5183568/150‑1 (MQ=255)
                                           cGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCa                                                                                                            <  2:4379151/150‑1 (MQ=255)
                                             gCACAACAAGTGGCGACCTTACACGGCGGTAAATTACCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCaat                                                                                                          >  2:886590/1‑149 (MQ=255)
                                               acaacaAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAAcg                                                                                                        >  1:922703/1‑150 (MQ=255)
                                                 aacaaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAAcgcg                                                                                                      >  1:4766677/1‑150 (MQ=255)
                                                 aacaaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAAcgcg                                                                                                      >  1:2165534/1‑150 (MQ=255)
                                                 aacaaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCGTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCTCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAAcgcg                                                                                                      >  1:4316187/1‑150 (MQ=255)
                                                  acaaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGt                                                                                                     >  2:1682751/1‑150 (MQ=255)
                                                   caaGTGGCGACCTTACACGGCGGTAAAGTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCTCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTg                                                                                                    <  1:248265/150‑1 (MQ=255)
                                                    aaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGc                                                                                                   >  2:747653/1‑150 (MQ=255)
                                                    aaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGc                                                                                                   <  1:5110657/150‑1 (MQ=255)
                                                    aaGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGc                                                                                                   >  2:5173802/1‑150 (MQ=255)
                                                     aGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCt                                                                                                  <  2:1414998/150‑1 (MQ=255)
                                                      gTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTg                                                                                                 <  1:4428215/150‑1 (MQ=255)
                                                      gTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTg                                                                                                 >  2:1599302/1‑150 (MQ=255)
                                                           gACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGgcgcg                                                                                            <  2:2569826/150‑1 (MQ=255)
                                                              cTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCgctg                                                                                         <  1:2858921/150‑1 (MQ=255)
                                                              cTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCgctg                                                                                         <  2:1981162/150‑1 (MQ=255)
                                                                 acacGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTa                                                                                      >  2:104167/1‑150 (MQ=255)
                                                                 acacGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTa                                                                                      <  2:1105504/150‑1 (MQ=255)
                                                                   acGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATg                                                                                    <  2:2109239/150‑1 (MQ=255)
                                                                     ggcggTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCt                                                                                  <  1:4614399/150‑1 (MQ=255)
                                                                     ggcggTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCt                                                                                  <  1:3789620/150‑1 (MQ=255)
                                                                       cggTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTACGCACGTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGt                                                                                >  2:587444/1‑150 (MQ=255)
                                                                        ggTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTa                                                                               >  1:5170390/1‑150 (MQ=255)
                                                                             aTTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCgg                                                                          <  1:5165609/150‑1 (MQ=255)
                                                                              ttCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGc                                                                         <  2:2313324/150‑1 (MQ=255)
                                                                               tCCCGGAAACCTTTGAGGAAGTTGCACCACTGCAGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCTCTTTCTCTGGGTAAGCACTTTCCGATTCTCGAAGGTAAAGTCAAACGCGTGCTGTCGCTCTGCTATGATGTAAGAggat                                                                        >  2:2731783/1‑148 (MQ=255)
                                                                                    gAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCt                                                                   >  1:594751/1‑150 (MQ=255)
                                                                                         cTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGac                                                              >  2:4259358/1‑149 (MQ=255)
                                                                                         cTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGaa                                                              <  2:1263634/150‑1 (MQ=255)
                                                                                            tGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGaaaaa                                                           >  1:1856955/1‑150 (MQ=255)
                                                                                            tGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGaaaaa                                                           >  1:769636/1‑150 (MQ=255)
                                                                                             gAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGaaaaaa                                                          >  1:3830115/1‑150 (MQ=255)
                                                                                                gAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGACAAAtgag                                                       >  1:2496757/1‑150 (MQ=255)
                                                                                                gAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAagag                                                       <  2:4836149/150‑1 (MQ=255)
                                                                                                  aGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGt                                                     >  2:1576516/1‑150 (MQ=255)
                                                                                                   gTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTc                                                    >  1:3605920/1‑150 (MQ=255)
                                                                                                    ttGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCg                                                   <  1:2694316/150‑1 (MQ=255)
                                                                                                    ttGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGAGTCTCTCTCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCg                                                   <  2:2336313/150‑1 (MQ=255)
                                                                                                     tGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGACAAAAGAGGTCga                                                  >  1:4569202/1‑150 (MQ=255)
                                                                                                      gcagcaCTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCgag                                                 <  2:1184900/150‑1 (MQ=255)
                                                                                                      gcagcaCTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCgag                                                 <  2:2803623/150‑1 (MQ=255)
                                                                                                       cagcaCTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCGGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAACAAAGAGGTCgaga                                                >  1:2454562/1‑150 (MQ=255)
                                                                                                        agcaCTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCTCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGaa                                               >  1:4505275/1‑150 (MQ=255)
                                                                                                        agcaCTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGAGTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGaa                                               <  1:1868080/150‑1 (MQ=255)
                                                                                                         gcaCTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAAt                                              <  2:3933624/150‑1 (MQ=255)
                                                                                                         gcaCTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAAt                                              >  1:1399686/1‑150 (MQ=255)
                                                                                                            cTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATaaa                                           <  2:430801/150‑1 (MQ=255)
                                                                                                            cTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATaaa                                           <  2:4646808/150‑1 (MQ=255)
                                                                                                             tGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAACAAAGAGGTCGAGAATAAAt                                          >  2:4542868/1‑150 (MQ=255)
                                                                                                               ccGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGACTAAATTa                                        >  1:134169/1‑150 (MQ=255)
                                                                                                               ccGGGCGTCGGGCGTTCCACCGCAGGCGCGAGTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTa                                        <  2:3126501/150‑1 (MQ=255)
                                                                                                                 gggcgtcgggcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTCTg                                      >  1:3360190/1‑150 (MQ=255)
                                                                                                                  ggcgtcgggcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATgg                                     <  1:911335/150‑1 (MQ=255)
                                                                                                                     gtcgggcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGt                                  <  2:1526470/150‑1 (MQ=255)
                                                                                                                        gggcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTg                               <  1:2852940/150‑1 (MQ=255)
                                                                                                                          gcgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgag                             <  1:555471/150‑1 (MQ=255)
                                                                                                                          gcgtTCCACCGCAGGCGCGATTATCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgag                             <  1:3134936/150‑1 (MQ=255)
                                                                                                                           cgtTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgagc                            >  1:564822/1‑150 (MQ=255)
                                                                                                                              tCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgagcgag                         <  1:2479652/150‑1 (MQ=255)
                                                                                                                               ccACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgagcgagc                        <  2:2602714/150‑1 (MQ=255)
                                                                                                                               ccACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTgagcgagc                        <  2:445111/150‑1 (MQ=255)
                                                                                                                                     cAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGa                  <  2:2685863/150‑1 (MQ=255)
                                                                                                                                      aGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGAc                 <  1:2503297/150‑1 (MQ=255)
                                                                                                                                          gcgATTCTCTCTCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGccc             <  2:4316187/150‑1 (MQ=255)
                                                                                                                                          gcgATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGccc             <  2:4084083/150‑1 (MQ=255)
                                                                                                                                          gcgATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGccc             <  1:3401392/150‑1 (MQ=255)
                                                                                                                                            gATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCcgc           >  2:4495485/1‑150 (MQ=255)
                                                                                                                                            gATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCcgc           >  1:4933745/1‑150 (MQ=255)
                                                                                                                                              ttCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCgg         >  2:2924603/1‑150 (MQ=255)
                                                                                                                                               tctctcGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGt        >  1:4729647/1‑150 (MQ=255)
                                                                                                                                               tctctcGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGt        >  1:998800/1‑150 (MQ=255)
                                                                                                                                                ctctcGCTTTCTCTGGGTACGCACTGTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGtt       <  2:3196944/150‑1 (MQ=255)
                                                                                                                                                ctctcGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGtt       <  1:1810487/150‑1 (MQ=255)
                                                                                                                                                ctctcGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGtt       <  1:3435233/150‑1 (MQ=255)
                                                                                                                                                   tcGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGc    <  1:4127303/150‑1 (MQ=255)
                                                                                                                                                     gCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGt  <  1:236950/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCGGGCTTGGCTATTACGCCCGCGCGCGCAATCTGCATAAAGCGGCACAACAAGTGGCGACCTTACACGGCGGTAAATTCCCGGAAACCTTTGAGGAAGTTGCAGCACTGCCGGGCGTCGGGCGTTCCACCGCAGGCGCGATTCTCTCGCTTTCTCTGGGTAAGCACTTTCCGATTCTCGACGGTAACGTCAAACGCGTGCTGGCGCGCTGCTATGCTGTAAGCGGCTGGCCTGGGAAAAAAGAGGTCGAGAATAAATTATGGAGTTTGAGCGAGCAGGTGACGCCCGCGGTTGGCGT  >  CP009273/3096603‑3096901

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