Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,695,395 A→G 5.5% Y61C (TAC→TGC)  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,3950AG 5.5% 219.3 / 2.7 79Y61C (TAC→TGC) malYbifunctional beta‑cystathionase, PLP‑dependent/ regulator of maltose regulon
Reads supporting (aligned to +/- strand):  ref base A (30/44);  new base G (3/2);  total (33/46)
Fisher's exact test for biased strand distribution p-value = 6.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

TGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAA  >  CP009273/1695253‑1695541
                                                                                                                                              |                                                                                                                                                  
tGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCGCTGCTGACCTGTTAGCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTAccgccg                                                                                                                                              >  1:654296/1‑149 (MQ=255)
  gtgtACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCt                                                                                                                                            <  1:621051/149‑1 (MQ=255)
    gtACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                                                                                          <  2:800094/149‑1 (MQ=255)
    gtACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTgg                                                                                                                                          <  2:632811/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:406644/149‑1 (MQ=255)
                              cGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTAt                                                                                                                <  1:751550/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                                                                                                   <  2:633745/149‑1 (MQ=255)
                                           ctgACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                                                                                                   <  1:1186512/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)
                                                                                                                  cAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTc                            <  2:673915/149‑1 (MQ=255)
                                                                                                                   aGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCa                           <  1:515606/149‑1 (MQ=255)
                                                                                                                    gcgcCTGATGCACGGCGTATTTGGCTGCAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACAGCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAg                          <  2:806768/149‑1 (MQ=255)
                                                                                                                    gcgcCTGATGCACGGCGTATTTGGCTGCAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAg                          <  2:184329/149‑1 (MQ=255)
                                                                                                                    gcgcCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCCCAGACGGTGGTGTATGGCCCTTCAGTCATCCATATggaatcac                          >  2:615970/1‑143 (MQ=255)
                                                                                                                    gcgcCTGATGCACGGCGTATTTGGCTAAAGGCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTCCCGCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTTTCATCTATATGGTTTCAg                          <  2:797927/149‑1 (MQ=255)
                                                                                                                      gcCCGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGaa                        <  1:599451/149‑1 (MQ=255)
                                                                                                                           aTGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATGTATATGGTTTCAGAACTGAt                   >  1:1141296/1‑149 (MQ=255)
                                                                                                                               aCGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGt               <  2:975561/149‑1 (MQ=255)
                                                                                                                                 ggCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCa             <  1:579837/149‑1 (MQ=255)
                                                                                                                                   cGTATTTGGTTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGt           <  1:976563/149‑1 (MQ=255)
                                                                                                                                    gTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTg          <  1:1114782/149‑1 (MQ=255)
                                                                                                                                     tATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTgg         >  2:1022416/1‑149 (MQ=255)
                                                                                                                                      aTTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGACCTGATTCGTCAGTGGt        >  1:1056544/1‑149 (MQ=255)
                                                                                                                                      aTTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGGTTCAGAACTGATTCGTCAGTGGt        >  1:1027409/1‑149 (MQ=255)
                                                                                                                                       tttGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACACACCATTACACCGCCATGGATTCTTCGTCTTCTCGCAACAACCCGGCTCGCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTc       >  2:875176/1‑149 (MQ=255)
                                                                                                                                           gCTACGGCCGCGGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGaa   <  1:458956/149‑1 (MQ=255)
                                                                                                                                           gCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGaa   <  2:556233/149‑1 (MQ=255)
                                                                                                                                           gCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTACGACGGTGGCGAAGGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGaa   >  1:595279/1‑149 (MQ=255)
                                                                                                                                            cTACAGCCGCTGTAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGaaa  <  2:672253/149‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
TGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTTTCAGAACTGATTCGTCAGTGGTCTGAAA  >  CP009273/1695253‑1695541

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: ‑