Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,354,610 A→T 9.2% W84R (TGG→AGG)  dsbD ← fused thiol:disulfide interchange protein: activator of DsbC/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,354,6100AT 9.2% 549.6 / 36.1 182W84R (TGG→AGG) dsbDfused thiol:disulfide interchange protein: activator of DsbC/conserved protein
Reads supporting (aligned to +/- strand):  ref base A (65/100);  new base T (11/6);  total (76/106)
Fisher's exact test for biased strand distribution p-value = 6.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

TTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTG  >  CP009273/4354463‑4354757
                                                                                                                                                   |                                                                                                                                                   
 ttCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAg                                                                                                                                                    <  1:2957049/148‑1 (MQ=255)
  tCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGa                                                                                                                                                   <  2:1651795/148‑1 (MQ=255)
     ggCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGc                                                                                                                                                <  2:1210351/148‑1 (MQ=255)
     ggCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGc                                                                                                                                                <  2:2205219/148‑1 (MQ=255)
      gCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGcc                                                                                                                                               <  1:1690900/148‑1 (MQ=255)
      gCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGcc                                                                                                                                               <  1:1559991/148‑1 (MQ=255)
       cGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCt                                                                                                                                              <  2:633614/148‑1 (MQ=255)
            aaCAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCgg                                                                                                                                         >  1:1010124/1‑148 (MQ=255)
               aGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCATAAGCTCGCGGTAAATAGAGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAg                                                                                                                                      <  1:3249126/148‑1 (MQ=255)
               aGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAg                                                                                                                                      <  1:2713856/148‑1 (MQ=255)
               aGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAg                                                                                                                                      >  2:3278187/1‑148 (MQ=255)
                gAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGACGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGc                                                                                                                                     >  2:2424516/1‑148 (MQ=255)
                gAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGc                                                                                                                                     <  1:311950/148‑1 (MQ=255)
                 aaaCCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCTGACGCCTTGCGGCAGCt                                                                                                                                    >  2:1095340/1‑148 (MQ=255)
                   aCCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGc                                                                                                                                  >  1:1664143/1‑148 (MQ=255)
                    ccGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCa                                                                                                                                 <  2:585650/148‑1 (MQ=255)
                    ccGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCa                                                                                                                                 <  2:1495820/148‑1 (MQ=255)
                    ccGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCa                                                                                                                                 <  2:585689/148‑1 (MQ=255)
                    ccGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCa                                                                                                                                 >  1:540721/1‑148 (MQ=255)
                     cGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCAc                                                                                                                                >  1:2464365/1‑148 (MQ=255)
                        cctCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTc                                                                                                                             <  2:558864/146‑1 (MQ=255)
                           cAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGc                                                                                                                          >  2:59505/1‑148 (MQ=255)
                              cacaGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAAt                                                                                                                       <  1:3206898/148‑1 (MQ=255)
                                caGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAAttt                                                                                                                     >  1:1106370/1‑148 (MQ=255)
                                caGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAAttt                                                                                                                     >  1:677546/1‑148 (MQ=255)
                                   cccTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTcg                                                                                                                  >  1:2228338/1‑148 (MQ=255)
                                       ggTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCTGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCgtg                                                                                                              <  1:1180805/148‑1 (MQ=255)
                                        gTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCgtgt                                                                                                             <  1:1675974/148‑1 (MQ=255)
                                           ggTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTcc                                                                                                          >  2:1123477/1‑148 (MQ=255)
                                           ggTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTcc                                                                                                          >  2:2775521/1‑148 (MQ=255)
                                            gTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCg                                                                                                         >  1:72911/1‑148 (MQ=255)
                                             tGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTCAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCgg                                                                                                        >  2:3230308/1‑148 (MQ=255)
                                               aCAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCg                                                                                                      >  1:906488/1‑148 (MQ=255)
                                                     aaCGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGc                                                                                                <  1:1310/148‑1 (MQ=255)
                                                      aCGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCTGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCg                                                                                               >  2:2019270/1‑148 (MQ=255)
                                                       cGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCgg                                                                                              <  2:2141454/148‑1 (MQ=255)
                                                       cGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCgg                                                                                              >  1:1001698/1‑148 (MQ=255)
                                                       cGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCgg                                                                                              >  1:1001728/1‑148 (MQ=255)
                                                         tCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGAt                                                                                            <  2:3042277/148‑1 (MQ=255)
                                                           gCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCt                                                                                          <  1:162951/148‑1 (MQ=255)
                                                              cccGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGtt                                                                                       <  2:2955977/148‑1 (MQ=255)
                                                              cccGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGtt                                                                                       <  1:1436077/148‑1 (MQ=255)
                                                                cGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTa                                                                                     <  1:967161/148‑1 (MQ=255)
                                                                cGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTa                                                                                     <  2:1003376/148‑1 (MQ=255)
                                                                cGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTa                                                                                     <  1:544148/148‑1 (MQ=255)
                                                                 gCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTAc                                                                                    >  1:76247/1‑148 (MQ=255)
                                                                 gCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTAc                                                                                    <  1:1275876/148‑1 (MQ=255)
                                                                    cTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGt                                                                                 >  2:1705286/1‑148 (MQ=255)
                                                                     tCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTa                                                                                >  1:1864818/1‑148 (MQ=255)
                                                                      cGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTag                                                                               <  2:1106370/148‑1 (MQ=255)
                                                                      cGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTag                                                                               >  2:1760688/1‑148 (MQ=255)
                                                                          tGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCTGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGgt                                                                           >  2:1582058/1‑148 (MQ=255)
                                                                           ggTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGgta                                                                          >  2:3181310/1‑148 (MQ=255)
                                                                             ttGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCTGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGgtagt                                                                        <  2:4406/148‑1 (MQ=255)
                                                                              tGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGgtagta                                                                       <  1:850611/148‑1 (MQ=255)
                                                                                aTGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAAc                                                                     >  1:2576230/1‑148 (MQ=255)
                                                                                   gTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGt                                                                  >  2:168600/1‑148 (MQ=255)
                                                                                    tGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTc                                                                 <  2:1666485/148‑1 (MQ=255)
                                                                                     gACGGGAAGCTTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCt                                                                <  2:1347182/148‑1 (MQ=255)
                                                                                     gACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCt                                                                <  1:59505/148‑1 (MQ=255)
                                                                                      aCGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCtt                                                               <  1:1763056/148‑1 (MQ=255)
                                                                                       cGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCttt                                                              <  2:2994336/148‑1 (MQ=255)
                                                                                        gggAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTg                                                             <  2:2353966/148‑1 (MQ=255)
                                                                                          gAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGAt                                                           <  1:2089826/148‑1 (MQ=255)
                                                                                            aGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCt                                                         >  2:612650/1‑148 (MQ=255)
                                                                                              cGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCTGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGc                                                       <  2:1167336/148‑1 (MQ=255)
                                                                                              cGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGc                                                       <  1:360887/148‑1 (MQ=255)
                                                                                                tCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCa                                                     <  1:2200993/148‑1 (MQ=255)
                                                                                                   gCCGATCGCGGTAACTCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGt                                                  <  2:3067123/148‑1 (MQ=255)
                                                                                                       aTCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGa                                              <  1:1496543/148‑1 (MQ=255)
                                                                                                        tCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGAt                                             >  1:2818846/1‑148 (MQ=255)
                                                                                                        tCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGAt                                             >  2:3254910/1‑148 (MQ=255)
                                                                                                          gcgGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGAGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTATCCGTCTTTGATCTGCCAGGTCAGATTa                                           <  1:3049511/148‑1 (MQ=255)
                                                                                                             gTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAgg                                        >  1:2645354/1‑148 (MQ=255)
                                                                                                             gTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAgg                                        <  1:613792/148‑1 (MQ=255)
                                                                                                               aaaTCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTc                                      <  1:2922673/148‑1 (MQ=255)
                                                                                                                 aTCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCAt                                    <  2:2391079/148‑1 (MQ=255)
                                                                                                                  tctcGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATg                                   >  1:1134234/1‑148 (MQ=255)
                                                                                                                   ctcGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGt                                  <  2:2936367/148‑1 (MQ=255)
                                                                                                                     cGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTg                                <  2:2251775/148‑1 (MQ=255)
                                                                                                                       cTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGt                              >  1:2020933/1‑148 (MQ=255)
                                                                                                                       cTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGt                              <  1:2775521/148‑1 (MQ=255)
                                                                                                                           tGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTg                          <  1:2522233/148‑1 (MQ=255)
                                                                                                                           tGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTg                          <  2:1430871/148‑1 (MQ=255)
                                                                                                                           tGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTg                          <  1:264358/148‑1 (MQ=255)
                                                                                                                              cGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTg                       >  1:2305854/1‑148 (MQ=255)
                                                                                                                                 aaaaCTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGaaa                    <  2:1470939/148‑1 (MQ=255)
                                                                                                                                  aaaCTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGaaaa                   >  2:3148605/1‑148 (MQ=255)
                                                                                                                                    aCTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATc                 >  2:807834/1‑148 (MQ=255)
                                                                                                                                     cTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCa                <  2:1864213/148‑1 (MQ=255)
                                                                                                                                     cTCATCTTCATGCCAGACGCCTTGCGGCAGCTGACCGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCa                <  2:2621923/148‑1 (MQ=255)
                                                                                                                                        aTCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCaaaa             <  2:1664143/148‑1 (MQ=255)
                                                                                                                                            tCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGAGTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCaa         >  1:771278/1‑148 (MQ=255)
                                                                                                                                             cATGCCTGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCaaa        >  2:2287797/1‑148 (MQ=255)
                                                                                                                                             cATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCaaa        <  2:1019430/148‑1 (MQ=255)
                                                                                                                                               tGCCTGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAgg      <  1:1595496/148‑1 (MQ=255)
                                                                                                                                                gCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGc     <  2:3008316/148‑1 (MQ=255)
                                                                                                                                                 ccAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCt    <  1:3262312/148‑1 (MQ=255)
                                                                                                                                                  cAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCtt   <  2:72911/148‑1 (MQ=255)
                                                                                                                                                   aGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTg  >  1:2257906/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TTTCTGGCGGATAACAGAAACCGGCATCAGCACAGCCCTGGTAGGTGACAGTTAACGTCGCTCCCGCACTCGCCTGGTTGATGGTGACGGGAAGCGTCAGCCGATCGCGGTAAATCTCGCTTTTGCCGTAAAACTCATCTTCATGCCAGACGCCTTGCGGCAGCTGCACGTCGGCAATTTTCGCGTGTTCCGGCGTAATGCGGATCTGTTTACGGTAGAGGTAGTAACCGTCTTTGATCTGCCAGGTCAGATTAAGGTCATGTTGGTTTTGCTGAAAATCAAAAGCAAAGGCTTG  >  CP009273/4354463‑4354757

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: