Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,898,349 G→A 8.7% F203F (TTC→TTT bglF ← fused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,898,3490GA 8.7% 695.0 / 38.2 206F203F (TTC→TTTbglFfused beta‑glucoside‑specific PTS enzymes: IIA component/IIB component/IIC component
Reads supporting (aligned to +/- strand):  ref base G (112/76);  new base A (10/8);  total (122/84)
Fisher's exact test for biased strand distribution p-value = 8.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 210 total aligned reads displayed.

CAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAA  >  CP009273/3898202‑3898495
                                                                                                                                                   |                                                                                                                                                  
 aGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGa                                                                                                                                                   >  2:1839537/1‑148 (MQ=255)
     aTAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATc                                                                                                                                               <  1:3094993/148‑1 (MQ=255)
       aGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCa                                                                                                                                             >  2:1058939/1‑148 (MQ=255)
         cAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGc                                                                                                                                           >  2:2772352/1‑148 (MQ=255)
               tgtgAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAgc                                                                                                                                     >  1:3270610/1‑148 (MQ=255)
                gtgAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAgcg                                                                                                                                    >  2:2687644/1‑148 (MQ=255)
                   aagaaATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCAt                                                                                                                                 <  1:3178220/148‑1 (MQ=255)
                   aagaaATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCAt                                                                                                                                 >  2:359302/1‑148 (MQ=255)
                   aagaaATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCAt                                                                                                                                 >  2:359349/1‑148 (MQ=255)
                    agaaATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATc                                                                                                                                >  1:617907/1‑148 (MQ=255)
                      aaaTTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCg                                                                                                                              >  2:1771005/1‑148 (MQ=255)
                          ttttGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCttt                                                                                                                          <  1:2764708/148‑1 (MQ=255)
                           tttGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCtttt                                                                                                                         <  1:2729850/148‑1 (MQ=255)
                           tttGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCtttt                                                                                                                         <  1:3250356/148‑1 (MQ=255)
                                  gCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGt                                                                                                                  <  1:2072224/148‑1 (MQ=255)
                                    cGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTtc                                                                                                                <  2:683019/148‑1 (MQ=255)
                                    cGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTtc                                                                                                                <  1:1423725/148‑1 (MQ=255)
                                       cGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCg                                                                                                             >  2:2114486/1‑148 (MQ=255)
                                         gTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGaa                                                                                                           <  2:1418598/148‑1 (MQ=255)
                                          tAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGaaa                                                                                                          >  2:1597798/1‑148 (MQ=255)
                                             ccACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCa                                                                                                       >  1:313316/1‑148 (MQ=255)
                                             ccACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCa                                                                                                       <  2:69120/148‑1 (MQ=255)
                                              cACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAg                                                                                                      <  1:1328011/148‑1 (MQ=255)
                                               aCGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGt                                                                                                     <  1:2397183/148‑1 (MQ=255)
                                               aCGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGt                                                                                                     <  1:555648/148‑1 (MQ=255)
                                                  cATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAg                                                                                                  >  1:240415/1‑148 (MQ=255)
                                                    ttAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGaa                                                                                                <  2:611548/148‑1 (MQ=255)
                                                    ttAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGaa                                                                                                <  2:618289/148‑1 (MQ=255)
                                                    ttAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACACAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGaa                                                                                                >  2:28366/1‑148 (MQ=255)
                                                     tAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAAt                                                                                               >  1:695498/1‑148 (MQ=255)
                                                      aaGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAAtt                                                                                              <  1:183766/148‑1 (MQ=255)
                                                       aGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTa                                                                                             <  1:2486911/148‑1 (MQ=255)
                                                        gTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTaa                                                                                            <  2:2135040/148‑1 (MQ=255)
                                                        gTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTaa                                                                                            <  2:1901314/148‑1 (MQ=255)
                                                         tCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAAt                                                                                           <  1:261188/148‑1 (MQ=255)
                                                          cGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATg                                                                                          >  2:1189996/1‑148 (MQ=255)
                                                            gCGTTCCAGAATGCGGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGa                                                                                        <  2:3230870/148‑1 (MQ=255)
                                                               ttCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGc                                                                                     >  2:1168098/1‑148 (MQ=255)
                                                               ttCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGc                                                                                     >  2:2294968/1‑148 (MQ=255)
                                                                tCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCa                                                                                    <  1:523725/148‑1 (MQ=255)
                                                                   aGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTa                                                                                 >  2:489248/1‑148 (MQ=255)
                                                                    gAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCTCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTaa                                                                                >  1:777897/1‑148 (MQ=255)
                                                                    gAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTaa                                                                                <  2:1716058/148‑1 (MQ=255)
                                                                     aaTGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAAc                                                                               >  1:1918179/1‑148 (MQ=255)
                                                                         ctgcACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTc                                                                           <  2:2574661/148‑1 (MQ=255)
                                                                              caACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTccacca                                                                      >  1:1694460/1‑148 (MQ=255)
                                                                                aCCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAAt                                                                    >  2:1399848/1‑148 (MQ=255)
                                                                                  cAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCa                                                                  >  1:1232643/1‑148 (MQ=255)
                                                                                  cAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCa                                                                  <  1:2586286/148‑1 (MQ=255)
                                                                                   aGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCAc                                                                 >  1:2377850/1‑148 (MQ=255)
                                                                                   aGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCAc                                                                 <  1:1203080/148‑1 (MQ=255)
                                                                                     gCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCa                                                               <  1:1106175/148‑1 (MQ=255)
                                                                                     gCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCa                                                               >  1:906409/1‑148 (MQ=255)
                                                                                     gCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCa                                                               >  1:2559323/1‑148 (MQ=255)
                                                                                        gAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATgg                                                            >  1:2451119/1‑148 (MQ=255)
                                                                                         aaaaaaTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGc                                                           >  2:1545396/1‑148 (MQ=255)
                                                                                         aaaaaaTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGc                                                           >  2:489759/1‑148 (MQ=255)
                                                                                         aaaaaaTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGc                                                           >  2:2617889/1‑148 (MQ=255)
                                                                                            aaaTAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGt                                                        >  1:2247071/1‑148 (MQ=255)
                                                                                              ataatGTGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTaa                                                      <  2:4331/148‑1 (MQ=255)
                                                                                                 atGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATg                                                   >  2:457084/1‑148 (MQ=255)
                                                                                                  tGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATgg                                                  <  2:2825367/148‑1 (MQ=255)
                                                                                                  tGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATgg                                                  <  1:3004401/148‑1 (MQ=255)
                                                                                                   gggAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGa                                                 >  2:600043/1‑148 (MQ=255)
                                                                                                    ggAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGAt                                                <  2:3122510/148‑1 (MQ=255)
                                                                                                     gAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGAtt                                               <  1:1521614/148‑1 (MQ=255)
                                                                                                      aaTAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTa                                              <  2:884511/148‑1 (MQ=255)
                                                                                                      aaTAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTa                                              >  1:347713/1‑148 (MQ=255)
                                                                                                      aaTAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTa                                              <  1:489243/148‑1 (MQ=255)
                                                                                                      aaTAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTa                                              <  1:489661/148‑1 (MQ=255)
                                                                                                           ccGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccgcc                                         >  2:2241720/1‑148 (MQ=255)
                                                                                                           ccGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTAccgcc                                         >  2:1896721/1‑148 (MQ=255)
                                                                                                             gATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGa                                       >  1:2048073/1‑148 (MQ=255)
                                                                                                              aTGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGaa                                      >  2:1481273/1‑148 (MQ=255)
                                                                                                              aTGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGaa                                      >  2:3101343/1‑148 (MQ=255)
                                                                                                              aTGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGaa                                      >  2:3147089/1‑148 (MQ=255)
                                                                                                              aTGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGaa                                      >  2:1237531/1‑148 (MQ=255)
                                                                                                               tGACGAGTAATTCAACAATGTGACCGGAATACCCAGAAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAg                                     >  1:1053774/1‑148 (MQ=255)
                                                                                                                gACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGc                                    <  1:256627/148‑1 (MQ=255)
                                                                                                                gACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGc                                    >  1:406117/1‑148 (MQ=255)
                                                                                                                 aCGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCg                                   <  2:630263/148‑1 (MQ=255)
                                                                                                                 aCGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCg                                   <  2:1273362/148‑1 (MQ=255)
                                                                                                                      tAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTcc                              <  1:1893713/148‑1 (MQ=255)
                                                                                                                      tAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTcc                              >  1:857266/1‑148 (MQ=255)
                                                                                                                         ttCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCcg                           >  1:2161664/1‑148 (MQ=255)
                                                                                                                          tCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCcgc                          >  1:969626/1‑148 (MQ=255)
                                                                                                                            aacaaTGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCgg                        >  1:377758/1‑148 (MQ=255)
                                                                                                                              caaTGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGgtg                      <  1:2940521/148‑1 (MQ=255)
                                                                                                                                aTGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTa                    >  1:1874882/1‑148 (MQ=255)
                                                                                                                                    gACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATccc                <  2:313341/148‑1 (MQ=255)
                                                                                                                                     aCCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATcccc               >  2:115865/1‑147 (MQ=255)
                                                                                                                                       cGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAgg             <  2:2718674/148‑1 (MQ=255)
                                                                                                                                       cGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAgg             <  1:2687644/148‑1 (MQ=255)
                                                                                                                                          aaTACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGAtt          <  1:2346547/148‑1 (MQ=255)
                                                                                                                                            tACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTAt        <  2:2059868/148‑1 (MQ=255)
                                                                                                                                             aCCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATc       >  2:2327633/1‑148 (MQ=255)
                                                                                                                                              cccAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCg      >  1:1614737/1‑148 (MQ=255)
                                                                                                                                                cAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCggg    >  2:3270320/1‑148 (MQ=255)
                                                                                                                                                cAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCggg    >  2:2403232/1‑148 (MQ=255)
                                                                                                                                                  ggAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGgaa  >  2:675086/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGACATAGCAATGGTGTGAAGAAATTTTTGATTGCCGACGGTAACCACGCATTAAGTCGGCGTTCCAGAATGCTGCACAACCAGGCAGAAAAAATAATGGGAATAACCGATGACGAGTAATTCAACAATGTGACCGGAATACCCAGGAAATCCAGCCCCAGCGCATCCGCTTTTTGCCCGTTCTCGAAAGCAGTCAGAATTAATGGATGCACTAACGCTCCACCAATCACCATGGCAGTAAATGGATTACCGCCGAAGCGTTTCCCCGCGGTGTATCCCAGGATTATCGGGAA  >  CP009273/3898202‑3898495

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: