Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,911,376 C→T 7.2% intergenic (+204/+28) barA → / ← gudD hybrid sensory histidine kinase, in two‑component regulatory system with UvrY/(D)‑glucarate dehydratase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,911,3760CT 7.2% 247.2 / 7.5 85intergenic (+204/+28)barA/gudDhybrid sensory histidine kinase, in two‑component regulatory system with UvrY/(D)‑glucarate dehydratase 1
Reads supporting (aligned to +/- strand):  ref base C (36/42);  new base T (2/4);  total (38/47)
Fisher's exact test for biased strand distribution p-value = 6.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.81e-01

ACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTT  >  CP009273/2911228‑2911519
                                                                                                                                                    |                                                                                                                                               
aCAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTcc                                                                                                                                                 >  1:1120482/1‑149 (MQ=255)
   aaGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGa                                                                                                                                              >  1:514149/1‑149 (MQ=255)
    aGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTATCATCAGTCTCAAACCGGCTCCAGAt                                                                                                                                             <  1:694684/149‑1 (MQ=255)
    aGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGAt                                                                                                                                             <  1:1058562/149‑1 (MQ=255)
    aGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGAt                                                                                                                                             <  1:1119132/149‑1 (MQ=255)
      attggttttgtaggCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATag                                                                                                                                           >  1:789498/14‑149 (MQ=255)
            tgtTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGCCTCAAACCGGCTCCAGAtgaagacgg                                                                                                                                     >  2:303019/1‑141 (MQ=255)
               ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggttt                                                                                                                                  <  2:179314/149‑1 (MQ=38)
               ttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttt                                                                                                                                  >  2:1015672/1‑149 (MQ=255)
                 attgcGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATCAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCACTCTCAAACCGGCTCTATCTGGAGCCggtttga                                                                                                                                <  2:889228/145‑3 (MQ=38)
                 aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCTATCTGGAGCCggtttga                                                                                                                                <  1:780432/149‑3 (MQ=38)
                 aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATGAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCATTCTCAAACCGGCTCTATCTGGAGCCggtttga                                                                                                                                <  2:889277/149‑3 (MQ=38)
                  tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttgg                                                                                                                               >  2:164122/1‑149 (MQ=255)
                     cGGATTCTGCGTTAACGCCTTATCCGGCCTACAAACCCAATTAAATTCAATGCATTAAAAATTGAGGTAGACCTGATAGATAGCACAACAGGCTATGTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggttt                                                                                                                            <  2:748803/149‑1 (MQ=255)
                      ggATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCggttttggtttt                                                                                                                           <  1:1148710/149‑1 (MQ=255)
                            ggcgTTAACGCCTTATCCGGCTCACAAAACCAATTAAATTCCATGCCTTCATCGCTCTTGTAGACGTGTTCCATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGTGCCGGTTTTGGTTTTCTGTCt                                                                                                                     <  1:236538/149‑1 (MQ=255)
                            ggcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATGAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGACCCGGTTCTGGCTTTCTGTCt                                                                                                                     >  2:84794/1‑149 (MQ=255)
                             gcgTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCtt                                                                                                                    <  1:566996/149‑1 (MQ=255)
                               gTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTaa                                                                                                                  <  2:1042640/149‑1 (MQ=255)
                                    cGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGCGCCGGTTTTGGTTTTCTGTCTTAACGCAc                                                                                                             <  2:11175/149‑1 (MQ=255)
                                    cGCCTAAACCGGCCAACAAAACCAAACAAAACCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTATCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCAc                                                                                                             <  2:50202/149‑1 (MQ=255)
                                        ttttCCGGCCCACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATg                                                                                                         <  2:1008315/146‑1 (MQ=255)
                                        ttATCCGGCCTACAAAACCAATTACATTCAATGAATTAAAAAAGGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATg                                                                                                         <  1:412265/149‑1 (MQ=255)
                                          aTCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGTGTTTCTGTCTTAACGCACACTGCa                                                                                                       >  1:1067235/1‑149 (MQ=255)
                                             cGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                    >  1:652081/1‑149 (MQ=255)
                                             cGGCCTAAAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGGCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACgg                                                                                                    <  1:268670/149‑1 (MQ=255)
                                              ggCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACggg                                                                                                   <  2:965310/149‑1 (MQ=255)
                                               gCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCCTGCACGGgc                                                                                                  >  2:671676/1‑149 (MQ=255)
                                               gCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAACACAACAGGCTATTTTGCGTTTGGCATAAGTCTCAAACCGACTCCAGATAGAGCCGGATTTGGTTTTCTGTCTTAACGCACCCTGCACGGgc                                                                                                  >  2:605552/1‑149 (MQ=255)
                                                 cTACAAAACCGATTAAATTCAATGAGTTAAAAGATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGgcgc                                                                                                <  2:716102/149‑1 (MQ=255)
                                                   aCAAAACCAATTAAATTCAATGACTTAAAAAATGATGTAGACCTCATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCTGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCtt                                                                                              <  2:778579/149‑1 (MQ=255)
                                                     aaaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGt                                                                                            >  2:56430/1‑149 (MQ=255)
                                                      aaaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCCCGGGCGCTTGtt                                                                                           >  1:1067108/1‑149 (MQ=255)
                                                       aaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                          >  1:603781/1‑149 (MQ=255)
                                                       aaCCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTa                                                                                          >  1:970740/1‑149 (MQ=255)
                                                        aCCAATTAAAGTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTAt                                                                                         <  1:809636/149‑1 (MQ=255)
                                                            aTTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGaa                                                                                     <  1:715974/149‑1 (MQ=255)
                                                             ttgaaTTAAAAGTATTAAAAAAGGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAAACGGCTCCAGCTAGGGCCGGTTTCGGTCTTCTGTCTTAATGCACCATTCACGGGCTCTTGTTATCGAAc                                                                                    <  2:421681/146‑1 (MQ=255)
                                                               aaaTTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATGTGACGGTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGt                                                                                  >  1:966927/1‑149 (MQ=255)
                                                                   tCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccag                                                                              <  2:370413/149‑1 (MQ=255)
                                                                    cAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagc                                                                             <  2:1067233/149‑1 (MQ=255)
                                                                     aaTGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTccagcc                                                                            >  1:290706/1‑149 (MQ=255)
                                                                       tGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCGGGCTTAACGCACCAGGCACGGGCGCTTGTTATCGAACGTcccgccag                                                                          >  1:115484/1‑149 (MQ=255)
                                                                           ttAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAAt                                                                      <  2:918419/149‑1 (MQ=255)
                                                                           ttAAAAAATGATGTAGACCGGATACATAGCACAACAGGCTATTTTGCGTTCAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAAt                                                                      <  1:286080/149‑1 (MQ=255)
                                                                               aaaaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGa                                                                  >  1:167829/1‑149 (MQ=255)
                                                                               aaaaTGACGTAGCCCTGATACATAGCACAACCGGCTATTTTGCTTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGa                                                                  <  1:734272/149‑1 (MQ=255)
                                                                                 aaTGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATa                                                                <  2:997332/149‑1 (MQ=255)
                                                                                    gatgTCGACCTGTTACATCGCCCAACAGGCTATTTTGCGTTGGGCATCAGTCTCAAACCGGCTCCAGATAGGGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTg                                                             <  2:551295/149‑1 (MQ=255)
                                                                                    gatgTAGCCTTGCTACATAGCACAACCGGCTATTTTGCGTTTCGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTg                                                             <  2:357527/149‑1 (MQ=255)
                                                                                    gatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTg                                                             <  2:1222088/149‑1 (MQ=255)
                                                                                    gatgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTg                                                             <  2:594211/149‑1 (MQ=255)
                                                                                     atgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGc                                                            <  1:15744/149‑1 (MQ=255)
                                                                                      tgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTGCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCa                                                           <  2:839259/149‑1 (MQ=255)
                                                                                      tgTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCa                                                           <  1:526987/149‑1 (MQ=255)
                                                                                      tgTAGACCAGATCCATATCACCACAGTCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTGGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCa                                                           <  2:917978/149‑1 (MQ=255)
                                                                                       gTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCag                                                          >  1:62174/1‑148 (MQ=255)
                                                                                           aCCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTccc                                                      >  2:865566/1‑149 (MQ=255)
                                                                                           aCCGGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTccc                                                      <  2:1141271/149‑1 (MQ=255)
                                                                                             cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACTGGCTCCGGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATCCTGCATTCCcat                                                    >  2:806817/1‑149 (MQ=255)
                                                                                             cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGACCGTCCAGCCCGGAATCAGATCCTGCATTCCcat                                                    >  2:1051911/1‑149 (MQ=255)
                                                                                             cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCcat                                                    <  2:1120482/149‑1 (MQ=255)
                                                                                             cTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACAGGCACCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAAAGTCCAGCCAGGAATCCGAGAAGTCATTCCcag                                                    >  2:32317/1‑148 (MQ=255)
                                                                                              tGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTccccca                                                   >  2:283169/1‑146 (MQ=255)
                                                                                                     tAGCACAACAGGCTGTTTTGCGTTTAGCATCAGTCTGAAACCGGCTCCAGGTAGAGCCAGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGGGCTTGTTATGGATGGTCCGGCGTGGGATCAGATTCTGCATTGCcatggctgcgt                                            >  2:816440/1‑144 (MQ=255)
                                                                                                     tAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCACGAATCAGATACTGCATTCCCATCGCCTCGt                                            >  2:171526/1‑149 (MQ=255)
                                                                                                        cacaACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGGACCATGCACGGGCGCTTGTTATCGAACGCCCAGTCATGAATCAGATACTGCATTCCCATCGCATCGTCAc                                         <  2:211732/149‑1 (MQ=255)
                                                                                                          caacaGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACACACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCAcgc                                       >  1:503644/1‑149 (MQ=255)
                                                                                                                 ttATTTTGCGTTTTGCATCAGGCTCAAACCGGCTCGAGATAGAGGCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCaa                                <  1:198183/148‑1 (MQ=255)
                                                                                                                 cTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTACAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCaa                                <  1:686317/149‑1 (MQ=255)
                                                                                                                  tATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAg                               <  2:136783/149‑1 (MQ=255)
                                                                                                                               gCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGa                  <  1:580667/149‑1 (MQ=255)
                                                                                                                                  tCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCCAGCCCGTGTTTCTGATAc               >  1:1197926/1‑149 (MQ=255)
                                                                                                                                  tCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCCAGCCCGGGTTTCTGATAc               >  1:314821/1‑149 (MQ=255)
                                                                                                                                  tCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATAc               >  1:641291/1‑149 (MQ=255)
                                                                                                                                  tCAGTCCCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATAc               <  1:67061/149‑1 (MQ=255)
                                                                                                                                     gTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGc            >  1:501391/1‑149 (MQ=255)
                                                                                                                                     aaccaaaaCCGGCTCTACCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCCGGAATCAGATACTGCATTCCCATCGCATCGTGCCGCGCGCCACGCCCGTGTTTCTGAtgcggc            >  2:81871/6‑144 (MQ=255)
                                                                                                                                        tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCGCGTGTTTCTGATACAGCTCa         >  1:1205978/1‑149 (MQ=255)
                                                                                                                                        tcAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACAGCATTCCCATCGCATCGTCACACGCGCCAAGCCAATGTTTCAGATACACCTCa         >  1:471541/1‑149 (MQ=255)
                                                                                                                                         cAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACCGCTCAt        >  1:314024/1‑149 (MQ=255)
                                                                                                                                          aaaCCGGCTCTATCTGGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTCATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCACATg       >  2:905254/1‑149 (MQ=255)
                                                                                                                                          aaaCCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTCATGCAGCTCATg       >  2:122271/1‑149 (MQ=255)
                                                                                                                                             ccGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGGGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCCCGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGc    >  1:43378/1‑149 (MQ=255)
                                                                                                                                              cGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCt   <  1:1015672/149‑1 (MQ=255)
                                                                                                                                               ggCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCtt  <  2:789498/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ACAAAGTGCCTGTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTACAAAACCAATTAAATTCAATGAATTAAAAAATGATGTAGACCTGATACATAGCACAACAGGCTATTTTGCGTTTAGCATCAGTCTCAAACCGGCTCCAGATAGAGCCGGTTTTGGTTTTCTGTCTTAACGCACCATGCACGGGCGCTTGTTATCGAACGTCCAGCCAGGAATCAGATACTGCATTCCCATCGCATCGTCACGCGCGCCAAGCCCGTGTTTCTGATACAGCTCATGGGCTT  >  CP009273/2911228‑2911519

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