Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,819,067 C→A 18.6% intergenic (‑132/‑124) mltB ← / → srlA membrane‑bound lytic murein transglycosylase B/glucitol/sorbitol‑specific enzyme IIC component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,819,0670CA18.6% 788.3 / 119.3 268intergenic (‑132/‑124)mltB/srlAmembrane‑bound lytic murein transglycosylase B/glucitol/sorbitol‑specific enzyme IIC component of PTS
Reads supporting (aligned to +/- strand):  ref base C (102/116);  new base A (18/32);  total (120/148)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.89e-01

Only 100 of 271 total aligned reads displayed.

TAACGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGC  >  CP009273/2818919‑2819216
                                                                                                                                                    |                                                                                                                                                     
taccgtcgCTTGAACATGAGGGGTCCATTTAACAGATTCATCCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGACGCGGTAGTCACTGCCCGATACGGACTTTACCTAACTCAACTCATTCCCCTCGCTATCCTTTTATTCa                                                                                                                                                      <  2:2683688/144‑1 (MQ=255)
tAACGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCa                                                                                                                                                      <  2:3817983/150‑1 (MQ=255)
   cgacgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAAc                                                                                                                                                   >  2:386231/1‑150 (MQ=255)
    gacgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACt                                                                                                                                                  >  2:4145293/1‑150 (MQ=255)
    gacgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACt                                                                                                                                                  >  2:4145148/1‑150 (MQ=255)
    gacgCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACt                                                                                                                                                  <  1:2230916/150‑1 (MQ=255)
       gCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTc                                                                                                                                               <  1:974941/150‑1 (MQ=255)
         ttGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCaa                                                                                                                                             <  1:785698/150‑1 (MQ=255)
         ttGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCaa                                                                                                                                             <  1:1589576/150‑1 (MQ=255)
          tGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATACCCCTCGCTATCCTTTTATTCAAACTTTCaaa                                                                                                                                            >  1:943489/1‑150 (MQ=255)
            aaCATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAAtt                                                                                                                                          <  1:4359658/150‑1 (MQ=255)
               aTGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTaaa                                                                                                                                       >  1:3654904/1‑150 (MQ=255)
                 gAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCCTTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAat                                                                                                                                     <  2:720285/150‑1 (MQ=255)
                 gAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAat                                                                                                                                     >  1:3382462/1‑150 (MQ=255)
                 gAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAat                                                                                                                                     >  1:3373132/1‑150 (MQ=255)
                  aGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAata                                                                                                                                    <  1:2624715/150‑1 (MQ=255)
                    gggTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATAtt                                                                                                                                  <  2:1287139/150‑1 (MQ=255)
                     ggTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATAttt                                                                                                                                 <  1:1064/150‑1 (MQ=255)
                      gTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTa                                                                                                                                >  2:2790568/1‑150 (MQ=255)
                          aTTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCtt                                                                                                                            >  1:4006333/1‑150 (MQ=255)
                           tttAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCttt                                                                                                                           <  2:747249/150‑1 (MQ=255)
                            ttAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTc                                                                                                                          <  1:536186/150‑1 (MQ=255)
                              aaCAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCAt                                                                                                                        <  1:2662209/150‑1 (MQ=255)
                              aaCAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCAt                                                                                                                        >  1:1943165/1‑150 (MQ=255)
                               aCAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCAtt                                                                                                                       <  1:16105/150‑1 (MQ=255)
                               aCAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCAtt                                                                                                                       <  1:3269304/150‑1 (MQ=255)
                                      cAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGAt                                                                                                                >  2:4614156/1‑150 (MQ=255)
                                      cAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGAt                                                                                                                <  1:737243/150‑1 (MQ=255)
                                          cAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCaaa                                                                                                            >  1:2150532/1‑150 (MQ=255)
                                           aGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCaaaa                                                                                                           >  1:4253836/1‑150 (MQ=255)
                                                cAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAaca                                                                                                      >  1:2593012/1‑150 (MQ=255)
                                                  aGTATGGTAAAGCATCACGCCCCGCACCAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACt                                                                                                    <  2:2572155/150‑1 (MQ=255)
                                                  aGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACt                                                                                                    <  2:3319022/150‑1 (MQ=255)
                                                   gTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACtt                                                                                                   >  1:3818998/1‑150 (MQ=255)
                                                   gTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACtt                                                                                                   >  2:1452678/1‑150 (MQ=255)
                                                   gTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACtt                                                                                                   <  1:244445/150‑1 (MQ=255)
                                                   gTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATACCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATATTTCATTTTGCGATCAAAATAACACtt                                                                                                   >  2:2784882/1‑150 (MQ=255)
                                                    tATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACttt                                                                                                  >  2:2834639/1‑150 (MQ=255)
                                                     aTGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACtttt                                                                                                 <  1:685134/150‑1 (MQ=255)
                                                     aTGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACtttt                                                                                                 >  1:2994779/1‑150 (MQ=255)
                                                      tGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAGACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTa                                                                                                <  1:615569/150‑1 (MQ=255)
                                                      tGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTa                                                                                                <  1:3702188/150‑1 (MQ=255)
                                                      tGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTa                                                                                                <  2:684804/150‑1 (MQ=255)
                                                        gTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTCTCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTaaa                                                                                              >  2:868456/1‑150 (MQ=255)
                                                         tAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAAt                                                                                             >  1:2385481/1‑150 (MQ=255)
                                                         tAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAAt                                                                                             >  2:1537476/1‑150 (MQ=255)
                                                          aaaGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATc                                                                                            >  2:2144647/1‑150 (MQ=255)
                                                            aGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCtt                                                                                          <  2:2278183/150‑1 (MQ=255)
                                                             gCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCttt                                                                                         <  2:1408857/150‑1 (MQ=255)
                                                             gCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAAATCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCttt                                                                                         <  2:2786950/150‑1 (MQ=255)
                                                              cATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTc                                                                                        >  1:4159974/1‑150 (MQ=255)
                                                               aTCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCa                                                                                       >  1:2687576/1‑150 (MQ=255)
                                                                 cACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAAt                                                                                     >  1:3918890/1‑150 (MQ=255)
                                                                 cACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTGTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAAt                                                                                     <  1:4573255/150‑1 (MQ=255)
                                                                   cGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCt                                                                                   <  1:3920078/150‑1 (MQ=255)
                                                                    gCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTg                                                                                  <  2:522035/150‑1 (MQ=255)
                                                                       ccGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTCCATAACTCACCTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTgatt                                                                               <  1:4538757/150‑1 (MQ=255)
                                                                       ccGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTgatt                                                                               <  1:2226122/150‑1 (MQ=255)
                                                                        cGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgatta                                                                              <  1:468524/150‑1 (MQ=255)
                                                                        cGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTgagta                                                                              >  2:584433/1‑150 (MQ=255)
                                                                            caAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTgattagatt                                                                          >  1:1554828/1‑150 (MQ=255)
                                                                               ggAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAgg                                                                       >  2:3560166/1‑150 (MQ=255)
                                                                                gAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGt                                                                      <  1:1475900/150‑1 (MQ=255)
                                                                                 aaGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGtt                                                                     >  1:1983968/1‑150 (MQ=255)
                                                                                 aaGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGtt                                                                     <  2:1984096/150‑1 (MQ=255)
                                                                                      gTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGt                                                                <  2:2341911/150‑1 (MQ=255)
                                                                                      gTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGt                                                                >  2:1540174/1‑150 (MQ=255)
                                                                                       tAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGtt                                                               >  1:2182868/1‑150 (MQ=255)
                                                                                       tAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTAGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGtt                                                               <  2:397724/150‑1 (MQ=255)
                                                                                        aGTCACTGCCCGATACGGACTTTACATAACTCAACTCTTTCCCCTCGCTATCCTTTTATTCAAACTTTCTAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAATTCTTTCAATCTGATTAGATTAGGGTGCCGttt                                                              <  2:1595786/150‑1 (MQ=255)
                                                                                        aGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGttt                                                              <  1:3560166/150‑1 (MQ=255)
                                                                                             cTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGtaa                                                         <  2:3147862/150‑1 (MQ=255)
                                                                                                cccGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGACGGTTGGtaataa                                                      >  2:2418299/1‑150 (MQ=255)
                                                                                                    aTACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACa                                                  <  2:1841291/150‑1 (MQ=255)
                                                                                                          aCTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAAt                                            <  1:2914809/150‑1 (MQ=255)
                                                                                                          aCTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAAt                                            >  1:1418041/1‑150 (MQ=255)
                                                                                                            tttACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATcc                                          >  2:133393/1‑150 (MQ=255)
                                                                                                               aCATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATTTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGa                                       >  2:2984105/1‑150 (MQ=255)
                                                                                                               aCATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGa                                       >  2:3540470/1‑150 (MQ=255)
                                                                                                               aCATAACTCAACTCATTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGa                                       >  1:40360/1‑150 (MQ=255)
                                                                                                                  tAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAgg                                    <  2:2593008/150‑1 (MQ=255)
                                                                                                                      tcaactcaTTCCCCTCGCTATCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGgagag                                <  2:64183/150‑1 (MQ=255)
                                                                                                                               tCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGAt                       <  2:433570/150‑1 (MQ=255)
                                                                                                                                    tCGCTATCCTTGTATTCAAACTTTCAAATTAAAATATATATCTATCATTTTGCGATCAAAATAACACTTTTAAATCATTCAAACTGATTAGATTAGGTTGCCGTTTGGTAAGAAAACAATAAATCCTGAAGGAGAGAACAATGATAGaaa                  >  2:365624/1‑150 (MQ=255)
                                                                                                                                       cTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCa               <  1:3863724/150‑1 (MQ=255)
                                                                                                                                          tCCTTTTATTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTa            <  1:2819333/150‑1 (MQ=255)
                                                                                                                                           ccTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTAc           >  1:994617/1‑150 (MQ=255)
                                                                                                                                           ccTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTAc           <  2:439853/150‑1 (MQ=255)
                                                                                                                                            cTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACt          <  2:188450/150‑1 (MQ=255)
                                                                                                                                            cTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACt          >  2:335074/1‑150 (MQ=255)
                                                                                                                                              tttATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCa        >  1:3052953/1‑150 (MQ=255)
                                                                                                                                              tttATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCa        >  2:1124795/1‑150 (MQ=255)
                                                                                                                                               ttATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCAt       >  2:945995/1‑150 (MQ=255)
                                                                                                                                                 aTTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATgg     >  1:1674485/1‑150 (MQ=255)
                                                                                                                                                 aTTAAAACTTTCAAATTAAAATATTTATCTTTCATTTTTCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATgg     >  1:2496979/1‑150 (MQ=255)
                                                                                                                                                  ttCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGt    >  1:1812010/1‑150 (MQ=255)
                                                                                                                                                  ttCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTCGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGt    <  2:3593905/150‑1 (MQ=255)
                                                                                                                                                  ttCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCCGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGt    <  1:794256/150‑1 (MQ=255)
                                                                                                                                                   tCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTg   >  1:1126412/1‑150 (MQ=255)
                                                                                                                                                    cAAACTTTCAAATTAAAATATTTATCTTTCATTATGCGAACAAAATAACACATTTAAATCTTTCAATCTGATTAGATTAGGTTGCCTTTAGGTAATAAAACAATAAATCCAGAAGGAGAGAACAATGATAGTAACCATTAATCATGGttc  >  2:4411280/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TAACGACGCTTGAACATGAGGGGTCCATTTAACAGATTCAACCAGGGGCAAGTATGGTAAAGCATCACGCCCCGCACAAGGAAGCGGTAGTCACTGCCCGATACGGACTTTACATAACTCAACTCATTCCCCTCGCTATCCTTTTATTCAAACTTTCAAATTAAAATATTTATCTTTCATTTTGCGATCAAAATAACACTTTTAAATCTTTCAATCTGATTAGATTAGGTTGCCGTTTGGTAATAAAACAATAAATCCTGAAGGAGAGAACAATGATAGAAACCATTACTCATGGTGC  >  CP009273/2818919‑2819216

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: