Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 376,412 G→A 9.0% pseudogene (304/416 nt) yaiX ← pseudogene, interrupted by IS2A, acetyltransferase homolog; nonfunctional; interruped by IS2; putative transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273376,4120GA 9.0% 974.1 / 59.8 287pseudogene (304/416 nt)yaiXpseudogene, interrupted by IS2A, acetyltransferase homolog; nonfunctional; interruped by IS2; putative transferase
Reads supporting (aligned to +/- strand):  ref base G (142/119);  new base A (12/14);  total (154/133)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 298 total aligned reads displayed.

AAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCT  >  CP009273/376263‑376561
                                                                                                                                                     |                                                                                                                                                     
aaaaaTGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAg                                                                                                                                                       <  1:3117523/150‑1 (MQ=255)
aaaaaTGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAg                                                                                                                                                       >  2:2227088/1‑150 (MQ=255)
          gccgcTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAaataat                                                                                                                                             <  1:1150975/150‑1 (MQ=255)
          gccgcTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAaataat                                                                                                                                             <  1:2006627/150‑1 (MQ=255)
            cgcTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAAtaataa                                                                                                                                           <  2:74375/150‑1 (MQ=255)
             gcTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAAc                                                                                                                                          <  2:3242109/150‑1 (MQ=255)
             gcTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAAc                                                                                                                                          >  1:659617/1‑150 (MQ=255)
               taataaCGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACtt                                                                                                                                        <  1:3450858/150‑1 (MQ=255)
                 ataaCGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGt                                                                                                                                      >  2:2623282/1‑150 (MQ=255)
                     cGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGc                                                                                                                                  >  2:2041087/1‑150 (MQ=255)
                      gTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGcc                                                                                                                                 >  1:2266656/1‑150 (MQ=255)
                        ccGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCaa                                                                                                                               <  2:4168662/150‑1 (MQ=255)
                         cGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAg                                                                                                                              <  2:3012673/150‑1 (MQ=255)
                         cGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAg                                                                                                                              >  1:14060/1‑150 (MQ=255)
                          gTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAgg                                                                                                                             >  2:580559/1‑150 (MQ=255)
                           ttttATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGc                                                                                                                            >  1:2875222/1‑150 (MQ=255)
                               aTCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGa                                                                                                                        <  1:2213139/150‑1 (MQ=255)
                                       aaTCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGcc                                                                                                                <  2:1084276/150‑1 (MQ=255)
                                           tcCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGat                                                                                                            <  1:3276479/150‑1 (MQ=255)
                                            ccTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGata                                                                                                           <  1:766559/150‑1 (MQ=255)
                                             cTGTACGTATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGatat                                                                                                          >  1:966465/1‑150 (MQ=255)
                                             cTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGatat                                                                                                          <  2:2585011/150‑1 (MQ=255)
                                               gTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATaa                                                                                                        <  1:1314335/150‑1 (MQ=255)
                                                  cGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATaacaa                                                                                                     >  1:3997183/1‑150 (MQ=255)
                                                   ggATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATaacaac                                                                                                    <  2:3950128/150‑1 (MQ=255)
                                                   ggATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATaacaac                                                                                                    >  1:1671124/1‑150 (MQ=255)
                                                    gATAAGTTCTGGTCGGTGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAAcc                                                                                                   <  2:1672243/150‑1 (MQ=255)
                                                      tAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTa                                                                                                 >  1:582209/1‑150 (MQ=255)
                                                         gTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAAtt                                                                                              >  2:2664348/1‑150 (MQ=255)
                                                               gTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGc                                                                                        >  2:3858190/1‑150 (MQ=255)
                                                                tCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCa                                                                                       >  1:3567482/1‑150 (MQ=255)
                                                                  ggAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATc                                                                                     >  1:3520807/1‑150 (MQ=255)
                                                                  ggAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATc                                                                                     >  2:3573213/1‑150 (MQ=255)
                                                                  ggAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATc                                                                                     <  2:1782953/150‑1 (MQ=255)
                                                                   taGTGAATAAGTTCCGGTAGGTAAATTACTTTCTACAATCACGCGCGGGACAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATATCAACCTATTTTATCGCATcc                                                                                    <  2:3657722/149‑1 (MQ=255)
                                                                    aGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCg                                                                                   <  1:783946/150‑1 (MQ=255)
                                                                    aGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCg                                                                                   <  1:2754953/150‑1 (MQ=255)
                                                                      tGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGt                                                                                 >  2:2580697/1‑150 (MQ=255)
                                                                      tGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGt                                                                                 >  2:2514517/1‑150 (MQ=255)
                                                                       gAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTa                                                                                >  1:713253/1‑150 (MQ=255)
                                                                        aaTAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAg                                                                               >  2:1062392/1‑150 (MQ=255)
                                                                           aaGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGa                                                                            >  2:139082/1‑150 (MQ=255)
                                                                            aGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGAt                                                                           <  2:3029966/150‑1 (MQ=255)
                                                                             gTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATa                                                                          <  1:3240053/150‑1 (MQ=255)
                                                                               tCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAAt                                                                        <  2:2722074/150‑1 (MQ=255)
                                                                                ccGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAAtt                                                                       >  1:135938/1‑150 (MQ=255)
                                                                                   gtaggtaAATTACGTTCTACAATCACGCGCGGGCCAAGGTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTccc                                                                    <  1:1286017/150‑1 (MQ=255)
                                                                                    taggtaAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCct                                                                   <  1:894916/150‑1 (MQ=255)
                                                                                    taggtaAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCct                                                                   <  2:2783041/150‑1 (MQ=255)
                                                                                      ggtaAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCctct                                                                 >  2:3226465/1‑150 (MQ=255)
                                                                                      ggtaAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCctct                                                                 >  2:217679/1‑150 (MQ=255)
                                                                                       gnaAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTg                                                                >  2:2525949/1‑150 (MQ=255)
                                                                                         aaaTTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGa                                                              >  2:2429696/1‑150 (MQ=255)
                                                                                          aaTTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAg                                                             >  1:3700434/1‑150 (MQ=255)
                                                                                           aTTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGt                                                            <  1:3750895/150‑1 (MQ=255)
                                                                                             tACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTc                                                          >  1:2323878/1‑150 (MQ=255)
                                                                                               cGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGa                                                        >  1:1483251/1‑150 (MQ=255)
                                                                                               cGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGa                                                        >  2:2265659/1‑150 (MQ=255)
                                                                                                gTTCTACAATCCCGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGaa                                                       <  2:3353088/150‑1 (MQ=255)
                                                                                                   cTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAg                                                    <  2:624980/150‑1 (MQ=255)
                                                                                                    tACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGa                                                   >  1:1878831/1‑150 (MQ=255)
                                                                                                    tACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGa                                                   >  2:1623687/1‑150 (MQ=255)
                                                                                                    tACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGa                                                   <  2:1897724/150‑1 (MQ=255)
                                                                                                      cactctcgcgcgGGCCAAGGTGCGTGTTCGGGGAAATTATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGaca                                                 <  2:2822424/144‑1 (MQ=255)
                                                                                                         tCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACgg                                              >  2:906999/1‑150 (MQ=255)
                                                                                                         tCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACgg                                              >  1:3896371/1‑150 (MQ=255)
                                                                                                            cgcgcgGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCt                                           >  1:1100746/1‑150 (MQ=255)
                                                                                                            cgcgcgGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCt                                           >  1:1100717/1‑150 (MQ=255)
                                                                                                              cgcgGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGt                                         <  1:169345/150‑1 (MQ=255)
                                                                                                              cgcgGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGt                                         >  1:385236/1‑150 (MQ=255)
                                                                                                                cgGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTc                                       <  1:2727961/150‑1 (MQ=255)
                                                                                                                cgGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTc                                       <  2:586350/150‑1 (MQ=255)
                                                                                                                 gggCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCa                                      >  1:2943486/1‑150 (MQ=255)
                                                                                                                 gggCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCa                                      <  2:3891193/150‑1 (MQ=255)
                                                                                                                     cAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCa                                  <  1:608856/150‑1 (MQ=255)
                                                                                                                      aaGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAg                                 <  1:231023/150‑1 (MQ=255)
                                                                                                                            cGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAAGAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACgatg                           <  2:3197842/150‑1 (MQ=255)
                                                                                                                            cGTGTTCGGGGAAATAATTCGCCCAAGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACgatg                           >  1:2080344/1‑150 (MQ=255)
                                                                                                                             gtgtTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACgatga                          >  2:138204/1‑150 (MQ=255)
                                                                                                                               gtTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGAtt                        >  2:573244/1‑150 (MQ=255)
                                                                                                                                  cGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACt                     >  2:851747/1‑150 (MQ=255)
                                                                                                                                  cGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACt                     <  2:3135347/150‑1 (MQ=255)
                                                                                                                                  cGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACt                     <  2:3785687/150‑1 (MQ=255)
                                                                                                                                    gggAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTgg                   >  2:1335159/1‑150 (MQ=255)
                                                                                                                                    gggAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTgg                   >  2:41052/1‑150 (MQ=255)
                                                                                                                                     ggAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGgt                  <  1:513922/150‑1 (MQ=255)
                                                                                                                                      gAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGgta                 <  1:3797469/150‑1 (MQ=255)
                                                                                                                                       aaaTAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGgtac                <  2:2337098/150‑1 (MQ=255)
                                                                                                                                        aataatTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTATCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGgtacg               >  2:646325/1‑150 (MQ=255)
                                                                                                                                          taatTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGgtacgta             <  1:3314687/150‑1 (MQ=255)
                                                                                                                                           aatTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGgtacgtac            >  1:3136835/1‑150 (MQ=255)
                                                                                                                                              tCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTtg         <  2:3065833/150‑1 (MQ=255)
                                                                                                                                               cGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTtgt        >  2:872953/1‑150 (MQ=255)
                                                                                                                                                gCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTtggg       <  2:1497241/150‑3 (MQ=255)
                                                                                                                                                 cccAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTgc      <  1:1822472/150‑1 (MQ=255)
                                                                                                                                                  ccAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTgcg     <  1:451369/150‑1 (MQ=255)
                                                                                                                                                   cAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTgcgc    >  2:2939835/1‑150 (MQ=255)
                                                                                                                                                    aGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGcc   >  2:1708392/1‑150 (MQ=255)
                                                                                                                                                     ggCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCt  >  2:518037/1‑150 (MQ=255)
                                                                                                                                                     ggCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCt  <  1:2494273/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAAAATGGCCGCCGCTAATAACGTCCGTTTTATCATTTTAATCTCCTGTACGGATAAGTTCTTGTCGGAGTGAATAAGTTCCGGTAGGTAAATTACGTTCTACAATCACGCGCGGGCCAAGTTGCGTGTTCGGGGAAATAATTCGCCCAGGCAAAATAATAACTTGTACGCCAAGGCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCT  >  CP009273/376263‑376561

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: