Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,695,365 A→T 9.3% N51I (AAT→ATT)  malY → bifunctional beta‑cystathionase, PLP‑dependent/ regulator of maltose regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,695,3650AT 9.3% 210.4 / 11.7 74N51I (AAT→ATT) malYbifunctional beta‑cystathionase, PLP‑dependent/ regulator of maltose regulon
Reads supporting (aligned to +/- strand):  ref base A (31/35);  new base T (2/5);  total (34/40)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.50e-01

TTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTAT  >  CP009273/1695217‑1695507
                                                                                                                                                    |                                                                                                                                              
ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                                                                >  1:30837/1‑149 (MQ=255)
 tCGATTTTTCAAAGGTCGTGGATCGTCATGACACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAAt                                                                                                                                               >  1:977227/1‑149 (MQ=255)
    aTTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTAGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTATGCGACTGCCCCCTGCATTATCGAGGCGCTGAATccg                                                                                                                                            >  1:537021/1‑147 (MQ=255)
         tCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCt                                                                                                                                       <  1:506445/149‑1 (MQ=255)
          cAAAGGTCGTGGATCGTCATGGCACATGGTGTACACGGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtc                                                                                                                                      <  2:859622/149‑2 (MQ=255)
            aaGGTGGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTCCATTCTCGCGGCGCTGAATCAGCGCCTGAt                                                                                                                                    >  2:867591/1‑149 (MQ=255)
             aGGTCGTGGATCCGCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgat                                                                                                                                   <  1:1009553/149‑5 (MQ=255)
             aGGTCGCGGATCCTAAAGGCACATTGTGTACACAGTAGGATTACGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtcgat                                                                                                                                   <  1:14986/149‑5 (MQ=255)
              ggTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGGGGCGCTGAATCAGCGCCCGAtac                                                                                                                                  >  2:599451/1‑147 (MQ=255)
                tCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGGGGCGCTGAATCAGCGCCTGATGCAc                                                                                                                                >  1:276525/1‑149 (MQ=255)
                tCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCCGGTGCAc                                                                                                                                >  2:130362/1‑149 (MQ=255)
                   tGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGc                                                                                                                             <  2:316745/149‑1 (MQ=255)
                          tCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGTCCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCACCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTg                                                                                                                      >  2:812027/1‑149 (MQ=255)
                          tCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCGGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTg                                                                                                                      >  1:406191/1‑149 (MQ=255)
                           cATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTgg                                                                                                                     <  1:867879/149‑1 (MQ=255)
                           cATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCACCTGATGCACGGCGTATTTgg                                                                                                                     <  2:1214766/149‑1 (MQ=255)
                                    tGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCGCTGCTGACCTGTTAGCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTAccgccg                                                                                                            >  1:654296/1‑149 (MQ=255)
                                      gtgtACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCt                                                                                                          <  1:621051/149‑1 (MQ=255)
                                        gtACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                                                        <  2:632811/149‑1 (MQ=255)
                                        gtACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                                                        <  2:800094/149‑1 (MQ=255)
                                             caGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATAAAAATTCCATGCACGGCGTATTTGGCTACAGCCGCTGGaaaaa                                                                                                   <  1:913117/149‑1 (MQ=255)
                                             caGTGGGATTATGTCGCTGACAGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaaaa                                                                                                   <  1:861385/149‑1 (MQ=255)
                                             caGTGGGATTATGGCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaaaa                                                                                                   <  2:957691/149‑1 (MQ=255)
                                                 ggtagtATCTCGCTGGTCGTGTCGGCCTTGCTTACCTGTTGCCGTGCACGATTTCAGACATGGAGTTTGCCACTGCCCCCTGCATTAGCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGAt                                                                                               <  2:63475/144‑1 (MQ=255)
                                                 gggATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTCGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATGTGGCTACAGCCGCTGGAAAAACGAt                                                                                               >  1:764134/1‑149 (MQ=255)
                                                 gggATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTAGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATCTGGCTACAGCCGCTGGAAAAACGAt                                                                                               >  1:764007/1‑149 (MQ=255)
                                                   gATTATGTCGCTGACCGTTTCGGCACTGGAGACCTGTTACCGTTCACGATTTCAGACATGGATCTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTGCAGCCGCTGGAAAAACGATGa                                                                                             >  1:806768/1‑149 (MQ=255)
                                                   gATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTGCAGCCGCTGGAAAAACGATGa                                                                                             >  1:184329/1‑149 (MQ=255)
                                                         gTCGCTGACCGTTTCGGCCCTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTtct                                                                                       <  1:808256/149‑1 (MQ=255)
                                                         ctcccttatcGTTTCCGCACTGCTGACCTTTTACCCGCCCCGATTTCAGACAATGATTTCGCCACTGCCCCCTTCATTATCCAGGCGCTGAATCAGCGCCTGCTGCACGGCGCATTTGGCTACAGCCACTGGAAAAACGATGAGTTtct                                                                                       <  2:606249/140‑1 (MQ=255)
                                                          tCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACCACCGCTGGAAAAACGATGAGtgtttc                                                                                      >  1:1094738/1‑144 (MQ=255)
                                                          tCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGCAAAACGATGAGTTtctc                                                                                      >  1:720190/1‑149 (MQ=255)
                                                          tCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAACGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGCAAAAAGATGAGtctctc                                                                                      >  2:377547/1‑149 (MQ=255)
                                                            gCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATGAGTGCGTGATTCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGCGTTTGTcgc                                                                                    >  1:286792/1‑149 (MQ=255)
                                                                aCCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCt                                                                                >  2:579837/1‑149 (MQ=255)
                                                                  cGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTAt                                                                              <  1:751550/149‑1 (MQ=255)
                                                                  cGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTAt                                                                              <  1:406644/149‑1 (MQ=255)
                                                                     ttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGc                                                                           <  2:1131256/149‑1 (MQ=255)
                                                                     gtgggCACTCCGGACCTGTTACCGTGCACGATTGCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGc                                                                           <  2:1131250/146‑1 (MQ=255)
                                                                     gtCGGCACTGCTGGCCTGTTACCTTTCACGAGTTCAGACATGGATTTTGCCACGGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGc                                                                           <  1:867591/148‑1 (MQ=255)
                                                                        ggCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTTCCCa                                                                        >  1:115197/1‑149 (MQ=255)
                                                                           aCTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTg                                                                     <  2:977227/149‑1 (MQ=255)
                                                                            ctgctgACCTGTTACCGTTAACGATTTCAGACATGGAGTTTGCCACTTCCCCCTGCATTATCGAGGCGCTGAATCCGCGCCTGATGCACGGCGTATTTGGGTGCAGCCGCTGGAAAAACGATACGTTTCTCGCGGCTGTTGcaccatgg                                                                    >  1:211357/1‑142 (MQ=255)
                                                                               ctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                                                                 <  1:1186512/149‑1 (MQ=255)
                                                                               ctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                                                                 <  2:633745/149‑1 (MQ=255)
                                                                                tgATCTGTTACCGTTCACGTTTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGCATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtttt                                                                <  1:664580/149‑1 (MQ=255)
                                                                                tgACCTGTTACCGTTCACGATTTCAAACATGGATTTTGCCACTGCCCACTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGtttt                                                                >  1:694444/1‑149 (MQ=255)
                                                                                   ccTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCTCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAATGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCa                                                             <  1:800940/149‑1 (MQ=255)
                                                                                      gTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAccc                                                          <  2:30837/149‑1 (MQ=255)
                                                                                       ttACCGTTCATGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTGTCTCGCGGCTATTGCCCACGGGTTTTCCACCCa                                                         >  1:205222/1‑149 (MQ=255)
                                                                                       ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAACCGCTGAGTTTCCCAAGGCTATTGCCCACTGGTTTTCCACCCa                                                         >  1:165262/1‑149 (MQ=255)
                                                                                       ttACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa                                                         >  2:1171789/1‑149 (MQ=255)
                                                                                        tACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGATACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAg                                                        <  1:1164660/149‑1 (MQ=255)
                                                                                          ccGTCCACGATGTCAGACATGGATTTTGCCACTGCCCCCTGCATTAGCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCa                                                      <  2:523484/149‑1 (MQ=255)
                                                                                            gTTCACGCTTGCAGACATGGATTTTTCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAtt                                                    <  1:1205135/149‑1 (MQ=255)
                                                                                            gTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTAGCGGCTATTGCCCACTGGTTTTCCACCCAGCAtt                                                    <  1:1226416/149‑1 (MQ=255)
                                                                                             ttCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTa                                                   <  2:276525/149‑1 (MQ=255)
                                                                                                aCGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTacac                                                >  1:673915/1‑149 (MQ=255)
                                                                                                  gATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCg                                              >  2:467409/1‑149 (MQ=255)
                                                                                                  gATTTCAGACATGGATTTAGCAACTACCCCCTGCATTAGGGAGGCGCTGAATCAGCGCCTGAGGCCCGGAGTATTTGTCTACAGCCGCGGGAAAAAGAATGAGTTTCTCGCGGCTATGGCCCCAAGTACTTCCACCCAGCCTTACACCg                                              >  2:369275/1‑149 (MQ=255)
                                                                                                       cAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATc                                         <  2:1056514/149‑1 (MQ=255)
                                                                                                                      ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACCGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtg                          <  1:289654/149‑1 (MQ=255)
                                                                                                                      ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtg                          >  2:1183130/1‑149 (MQ=255)
                                                                                                                      ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtg                          <  1:469942/149‑1 (MQ=255)
                                                                                                                        aCTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAGCGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtggt                        <  2:246479/149‑1 (MQ=255)
                                                                                                                        aCTGCCCCCTGCATTATCGACGCGCTGAGTCCGCGCCTGATGCACGGCGTATTTGGCTACCGCCGCTGGAAAAAGGATGAGTTTCTCGAGGGTATTGCCCAATGGTTTTCTACCCAGCCTTACACAGCCATCGATTCTCAGACgggggt                        >  2:746887/1‑149 (MQ=255)
                                                                                                                         cTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGgtggtg                       >  1:14762/1‑149 (MQ=255)
                                                                                                                            cccccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTAt                    <  1:151614/149‑1 (MQ=255)
                                                                                                                                cTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGcc                >  1:1126901/1‑149 (MQ=255)
                                                                                                                                     ttgtCGAGGCGCTGAATCAGGGCCTGTTGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGTTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCt           <  2:809499/146‑1 (MQ=255)
                                                                                                                                     ttATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGATGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCt           <  2:1127419/149‑1 (MQ=255)
                                                                                                                                     tgcatcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCt           >  1:118778/7‑149 (MQ=255)
                                                                                                                                         cGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCa       <  1:130362/149‑1 (MQ=255)
                                                                                                                                          caGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCAt      >  2:1114782/2‑149 (MQ=255)
                                                                                                                                              cgcTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCtat  >  2:1173493/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
TTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTAT  >  CP009273/1695217‑1695507

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