Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,960,118:1 +C 65.6% coding (1676/2247 nt) ptsP ← fused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,960,1181.C65.6% 147.0 / 87.6 61coding (1676/2247 nt)ptsPfused PTS enzyme: PEP‑protein phosphotransferase (enzyme I)/GAF domain containing protein
Reads supporting (aligned to +/- strand):  ref base . (9/12);  new base C (21/19);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 5.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGT  >  CP009273/2959970‑2960265
                                                                                                                                                     |                                                                                                                                                   
cAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAc                                                                                                                                                    >  2:727696/1‑150 (MQ=255)
cAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAc                                                                                                                                                    >  1:126843/1‑150 (MQ=255)
 aGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAcc                                                                                                                                                   >  2:1144642/1‑150 (MQ=255)
  gATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAccc                                                                                                                                                  <  2:1383679/150‑2 (MQ=255)
    tGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGAcccat                                                                                                                                                <  1:1336114/150‑4 (MQ=255)
      cGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCgg                                                                                                                                              <  1:350418/150‑1 (MQ=255)
         cagcaTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGc                                                                                                                                             >  2:858054/1‑149 (MQ=255)
            caTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACa                                                                                                                                         <  2:479174/150‑1 (MQ=255)
             aTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCaa                                                                                                                                           >  1:808508/1‑147 (MQ=255)
             aTAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCaa                                                                                                                                           >  1:324452/1‑147 (MQ=255)
              tAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGaa                                                                                                                                      >  1:1131494/1‑150 (MQ=255)
                aaTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGaa                                                                                                                                      >  1:546267/1‑149 (MQ=255)
                 aTACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAatat                                                                                                                                   >  1:262962/1‑150 (MQ=255)
                    ccATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATAtt                                                                                                                                  >  2:1384928/1‑148 (MQ=255)
                      aTTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAg                                                                                                                               >  1:73903/1‑150 (MQ=255)
                      aTTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAg                                                                                                                               >  2:689228/1‑150 (MQ=255)
                            ggCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGccc                                                                                                                        <  1:1300333/150‑1 (MQ=255)
                             gCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGccc                                                                                                                        <  2:408363/150‑1 (MQ=255)
                              cACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCg                                                                                                                       <  2:711517/150‑1 (MQ=255)
                                   ccAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCg                                                                                                                  >  1:350749/1‑150 (MQ=255)
                                    cAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCg                                                                                                                  <  1:1219694/149‑1 (MQ=255)
                                     aGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGc                                                                                                                <  2:1450427/150‑1 (MQ=255)
                                     aGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCa                                                                                                               >  1:667477/1‑150 (MQ=255)
                                       cATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCAt                                                                                                              >  2:972031/1‑150 (MQ=255)
                                             gCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCAc                                                                                                        >  1:489829/1‑150 (MQ=255)
                                               cGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACg                                                                                                       >  2:1376507/1‑149 (MQ=255)
                                               cGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACg                                                                                                       >  2:1017142/1‑148 (MQ=255)
                                                gATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACgca                                                                                                     <  2:269254/150‑1 (MQ=255)
                                                gATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACgca                                                                                                     <  1:1144642/150‑1 (MQ=255)
                                                    cGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCAt                                                                                                 <  1:1384928/149‑1 (MQ=255)
                                                    cGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATc                                                                                                <  1:1385025/150‑1 (MQ=255)
                                                           tGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCg                                                                                           <  2:1276156/149‑1 (MQ=255)
                                                                aTTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTgg                                                                                    <  2:977058/150‑1 (MQ=255)
                                                                 tttCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGCGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTg                                                                                     >  1:825333/1‑149 (MQ=255)
                                                                    cGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGAt                                                                                  <  2:825333/149‑1 (MQ=255)
                                                                        aCCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCaag                                                                              <  2:1080981/149‑1 (MQ=255)
                                                                        aCCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCaagaa                                                                            <  2:355784/150‑1 (MQ=255)
                                                                         ccGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCaagaag                                                                           <  2:1131494/150‑1 (MQ=255)
                                                                           gATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCaag                                                                              >  1:393698/1‑146 (MQ=255)
                                                                             tcatcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtc                                                                        >  1:569199/1‑150 (MQ=255)
                                                                             tcatcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtct                                                                       >  1:1112180/1‑150 (MQ=255)
                                                                               atcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtctc                                                                      >  1:1176/1‑150 (MQ=255)
                                                                               atcTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGAtctc                                                                      <  1:1076460/149‑1 (MQ=255)
                                                                                    ctcGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTcc                                                                     <  2:1237041/146‑1 (MQ=255)
                                                                                    ctcGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTg                                                                <  2:262962/150‑1 (MQ=255)
                                                                                             aCGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGc                                                         <  1:972031/149‑1 (MQ=255)
                                                                                                tCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTaa                                                     <  1:1419952/150‑1 (MQ=255)
                                                                                                  cGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAAt                                                    <  2:350749/149‑1 (MQ=255)
                                                                                                     cACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCg                                                 >  2:1110236/1‑149 (MQ=255)
                                                                                                           cAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATccc                                           <  2:126843/149‑1 (MQ=255)
                                                                                                           cAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCa                                          >  2:915216/1‑150 (MQ=255)
                                                                                                               cAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGc                                       >  2:262137/1‑149 (MQ=255)
                                                                                                                aGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGcc                                      >  1:991033/1‑149 (MQ=255)
                                                                                                                aGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGcc                                      <  2:393698/149‑1 (MQ=255)
                                                                                                                    ggcgTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAAcc                                  >  2:655201/1‑148 (MQ=255)
                                                                                                                       gTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAgg                              <  2:808508/150‑1 (MQ=255)
                                                                                                                       gTGCTTCATCAACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAgg                              >  1:282872/1‑150 (MQ=255)
                                                                                                                            tcatcaACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATg                          >  2:391049/1‑149 (MQ=255)
                                                                                                                             catcaACTTCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATgg                         <  2:546267/149‑1 (MQ=255)
                                                                                                                             catcaACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGAt                       >  1:344508/1‑150 (MQ=255)
                                                                                                                                    ttCATCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCtt                 >  1:304648/1‑150 (MQ=255)
                                                                                                                                       aTCGAGGCTTGTGACCCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCg               <  1:689228/149‑1 (MQ=255)
                                                                                                                                            ggCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCg        <  1:727696/150‑1 (MQ=255)
                                                                                                                                             gCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCgg       <  2:758096/150‑1 (MQ=255)
                                                                                                                                               ttGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCa     >  2:464636/1‑150 (MQ=255)
                                                                                                                                                  tgACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGt  <  1:339010/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGATGCGGCAGCATAAATACCATTGACGGCACTTCCAGCATGATGCCGATACGTGGTTTGGGAATTTCGTAACCGATCATCTCCTCGACTTCACGTCCGGCACGTTCAATCAGGCGGCGTGCTTCATCAACTTCATCGAGGCTTGTGACCATCGGCAACAGAATATTCAGGTTGCCCGTAGCGGCATTAGCACGCAGCATCGCCCGCACCTGGATCAAGAAGATCTCCGGCTGATCGAGCGTAATGCGAATCCCACGCCAACCCAGGCATGGATTCTCTTCGCTGATCGGCATGT  >  CP009273/2959970‑2960265

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: