Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,681,873 T→A 85.7% Q97L (CAG→CTG)  glrR ← response regulator regulating glmY sRNA in two‑component system with sensor protein GlrK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,681,8730TA85.7% 916.3 / 128.4 293Q97L (CAG→CTG) glrRresponse regulator regulating glmY sRNA in two‑component system with sensor protein GlrK
Reads supporting (aligned to +/- strand):  ref base T (19/23);  new base A (103/148);  total (122/171)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 312 total aligned reads displayed.

CAAACGCAGCATCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATCTACTTTTTC  >  CP009273/2681725‑2682021
                                                                                                                                                    |                                                                                                                                                    
cAAACGCAGCATCAGCGGGCTGCGGGTGACATTTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCctgct                                                                                                                                                      <  2:1215931/149‑1 (MQ=255)
  aaCGCAGCATCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCtgt                                                                                                                                                    >  2:2706717/1‑149 (MQ=255)
    cGCAGCATCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAgtgt                                                                                                                                                  >  1:2348582/1‑149 (MQ=255)
     gcagcaTCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGtt                                                                                                                                                 <  2:3767000/149‑1 (MQ=255)
        gcaTCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTgca                                                                                                                                              >  2:1708581/1‑149 (MQ=255)
        gcaTCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTgca                                                                                                                                              <  2:2410920/149‑1 (MQ=255)
           tCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTgcagca                                                                                                                                           >  1:164609/1‑149 (MQ=255)
            cAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTTCAGCaa                                                                                                                                          <  2:3857994/149‑1 (MQ=255)
             aGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAAc                                                                                                                                         >  2:571942/1‑149 (MQ=255)
             aGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAAc                                                                                                                                         >  2:571911/1‑149 (MQ=255)
                 ggCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCa                                                                                                                                     >  1:4671336/1‑149 (MQ=255)
                  gCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCAt                                                                                                                                    >  1:1190182/1‑149 (MQ=255)
                        ggTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGaa                                                                                                                              <  2:2050415/149‑1 (MQ=255)
                        ggTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGaa                                                                                                                              >  1:136933/1‑149 (MQ=255)
                        ggTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTATTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGaa                                                                                                                              <  1:506616/149‑1 (MQ=255)
                         gTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAAt                                                                                                                             <  1:1902110/149‑1 (MQ=255)
                             cAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATCGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGaa                                                                                                                         <  1:877203/149‑1 (MQ=255)
                             cAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGaa                                                                                                                         >  1:2294654/1‑149 (MQ=255)
                             cAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGaa                                                                                                                         <  1:1515797/149‑1 (MQ=255)
                              aaTTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAAc                                                                                                                        >  2:4275123/1‑149 (MQ=255)
                               aTTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAAcc                                                                                                                       >  2:2018948/1‑149 (MQ=255)
                                ttGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCa                                                                                                                      <  2:1378827/149‑1 (MQ=255)
                                  gCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATg                                                                                                                    >  1:245142/1‑149 (MQ=255)
                                     tCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATgcgc                                                                                                                 >  1:2277013/1‑149 (MQ=255)
                                          ccAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTaa                                                                                                            <  2:4229902/149‑1 (MQ=255)
                                           cAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAg                                                                                                           >  2:1966702/1‑149 (MQ=255)
                                             gCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaa                                                                                                         <  1:4221683/149‑1 (MQ=255)
                                             gCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaa                                                                                                         <  2:1890508/149‑1 (MQ=255)
                                               gTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaata                                                                                                       <  2:717058/149‑1 (MQ=255)
                                               gTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaata                                                                                                       >  2:4197066/1‑149 (MQ=255)
                                                ttCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaataa                                                                                                      <  1:3066259/149‑1 (MQ=255)
                                                ttCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGaataa                                                                                                      <  1:435655/149‑1 (MQ=255)
                                                  cATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAAtt                                                                                                    <  2:363782/149‑1 (MQ=255)
                                                  cATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAAtt                                                                                                    <  2:4012406/149‑1 (MQ=255)
                                                   aTCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTa                                                                                                   <  2:4612961/149‑1 (MQ=255)
                                                   aTCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTa                                                                                                   <  2:2837020/149‑1 (MQ=255)
                                                   aTCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTa                                                                                                   <  2:3326484/149‑1 (MQ=255)
                                                   aTCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTa                                                                                                   <  2:858139/149‑1 (MQ=255)
                                                    tCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTAc                                                                                                  >  1:3394216/1‑149 (MQ=255)
                                                     cGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTAcc                                                                                                 <  1:606863/149‑1 (MQ=255)
                                                     cGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAAGTGCCTGATATCGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCCTGCGCGGTAAGAATAATTAcc                                                                                                 <  2:25868/149‑1 (MQ=255)
                                                      ggTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCg                                                                                                <  1:4127614/149‑1 (MQ=255)
                                                         ggctggcGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCa                                                                                             <  2:3936499/149‑1 (MQ=255)
                                                            tggcGCGGATTGCTCCAGCGCATTGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTc                                                                                          >  2:4416992/1‑149 (MQ=255)
                                                             ggcGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTcc                                                                                         >  2:2621921/1‑149 (MQ=255)
                                                              gcgcGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTccc                                                                                        <  1:4009220/149‑1 (MQ=255)
                                                               cgcgGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCg                                                                                       <  1:2956538/149‑1 (MQ=255)
                                                                gcgGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCgg                                                                                      <  2:3989878/149‑1 (MQ=255)
                                                                  ggATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCt                                                                                    >  2:4101318/1‑149 (MQ=255)
                                                                   gATTGCTCCAGCGCATCGTCAATTGCCTGGTATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTg                                                                                   >  1:1295386/1‑149 (MQ=255)
                                                                          ccAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTc                                                                            <  1:3656394/149‑1 (MQ=255)
                                                                           cAGCGCATCGTCATTTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCt                                                                           <  1:2814862/149‑1 (MQ=255)
                                                                           cAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCt                                                                           <  1:2921825/149‑1 (MQ=255)
                                                                           cAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCt                                                                           <  1:4579736/149‑1 (MQ=255)
                                                                           cAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCt                                                                           >  2:599667/1‑149 (MQ=255)
                                                                            aGCGCATCGTCATGTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTg                                                                          <  1:3441152/149‑1 (MQ=255)
                                                                            aGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAACAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTg                                                                          <  1:3553457/149‑1 (MQ=255)
                                                                               gcATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGAt                                                                       <  1:3705531/149‑1 (MQ=255)
                                                                               gcATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGAt                                                                       >  2:1992744/1‑149 (MQ=255)
                                                                                cATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGAtt                                                                      <  2:1248258/149‑1 (MQ=255)
                                                                                cATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGAtt                                                                      <  1:1884099/149‑1 (MQ=255)
                                                                                  tCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTc                                                                    <  1:1097334/149‑1 (MQ=255)
                                                                                   cGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCa                                                                   <  2:1789247/149‑1 (MQ=255)
                                                                                    gTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAg                                                                  >  2:608537/1‑149 (MQ=255)
                                                                                    gTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAg                                                                  >  2:4664980/1‑149 (MQ=255)
                                                                                    gTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAg                                                                  >  2:2323945/1‑149 (MQ=255)
                                                                                     tCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGc                                                                 >  1:4165471/1‑149 (MQ=255)
                                                                                     tCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGc                                                                 >  2:3001991/1‑149 (MQ=255)
                                                                                        aTTGCCTGATATCGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCaaa                                                              <  2:1518663/149‑1 (MQ=255)
                                                                                        aTTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCaaa                                                              >  1:2422266/1‑149 (MQ=255)
                                                                                             cTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCt                                                         <  1:2092183/149‑1 (MQ=255)
                                                                                              tGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTg                                                        >  2:4447422/1‑149 (MQ=255)
                                                                                                atatAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCa                                                      >  2:1065234/1‑149 (MQ=255)
                                                                                                 tattGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCAt                                                     >  2:3782771/1‑149 (MQ=255)
                                                                                                  ataGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATa                                                    <  1:4642128/149‑1 (MQ=255)
                                                                                                    aGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATAcc                                                  >  1:3869987/1‑149 (MQ=255)
                                                                                                      cgcgTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGt                                                >  1:3265115/1‑149 (MQ=255)
                                                                                                        cgTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTcc                                              <  1:316951/149‑1 (MQ=255)
                                                                                                         gTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCa                                             >  1:1686880/1‑149 (MQ=255)
                                                                                                          tCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCAt                                            <  1:73788/149‑1 (MQ=255)
                                                                                                          tCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCAt                                            <  1:2148010/149‑1 (MQ=255)
                                                                                                           cTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCAtt                                           <  1:1798964/149‑1 (MQ=255)
                                                                                                            tttGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCAttt                                          >  1:3603286/1‑149 (MQ=255)
                                                                                                               gTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTcat                                       <  2:2169589/149‑1 (MQ=255)
                                                                                                                  gACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTcatcca                                    <  1:813857/149‑1 (MQ=255)
                                                                                                                      ggCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCg                                >  1:2681645/1‑149 (MQ=255)
                                                                                                                      ggCTTGGTGAGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCg                                >  1:4501338/1‑149 (MQ=255)
                                                                                                                               aGGAAACTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCt                       <  1:250577/149‑1 (MQ=255)
                                                                                                                                  aaaCTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgat                    <  1:4547788/149‑1 (MQ=255)
                                                                                                                                   aaCTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgatg                   <  1:1966624/149‑1 (MQ=255)
                                                                                                                                   aaCTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgatg                   >  1:3942385/1‑149 (MQ=255)
                                                                                                                                   aaCTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCtgatg                   >  1:125052/1‑149 (MQ=255)
                                                                                                                                     cTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGAc                 <  1:334300/149‑1 (MQ=255)
                                                                                                                                     cTAAAAACGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGAc                 >  2:4374980/1‑149 (MQ=255)
                                                                                                                                        aaaaCGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTaa              <  1:4244245/149‑1 (MQ=255)
                                                                                                                                         aaaCGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTaaa             <  2:245142/149‑1 (MQ=255)
                                                                                                                                           aCGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATc           <  1:4405773/149‑1 (MQ=255)
                                                                                                                                           aCGCCCTGCAGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATc           <  2:4523987/149‑1 (MQ=255)
                                                                                                                                                ctgcagTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATCTACtt      <  2:1295386/149‑1 (MQ=255)
                                                                                                                                                  gctgTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATCTACtttt    <  2:2294638/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAACGCAGCATCAGCGGGCTGCGGGTGACAATTGCCTCGCGCCAGCGTTCATCGGTGGCTGGCGCGGATTGCTCCAGCGCATCGTCAATTGCCTGATATAGCGCGTCTTTGTCGACAGGCTTGGTGAGGAAACTAAAAACGCCCTGCTGTGTTGCAGCAACGGCATCGGGAATAGAACCATGCGCGGTAAGAATAATTACCGGCATTCCCGGCTGCACTTTCTGGATTTCAGCAAACAGCTGCATACCGTCCATTTCATCCATCCGCAGGTCGCTGATGACTAAATCTACTTTTTC  >  CP009273/2681725‑2682021

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: