Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,770,878 C→A 12.5% I345I (ATC→ATA ydiF → putative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,770,8780CA12.5% 1111.5 / 86.2 344I345I (ATC→ATAydiFputative acetyl‑CoA:acetoacetyl‑CoA transferase: alpha subunit/beta subunit
Reads supporting (aligned to +/- strand):  ref base C (130/171);  new base A (17/26);  total (147/197)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

Only 100 of 352 total aligned reads displayed.

TTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAA  >  CP009273/1770729‑1771027
                                                                                                                                                     |                                                                                                                                                     
     cGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCtt                                                                                                                                                  >  2:3991415/1‑150 (MQ=255)
     cGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCtt                                                                                                                                                  <  2:255056/150‑1 (MQ=255)
      gAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCttt                                                                                                                                                 >  2:3243606/1‑150 (MQ=255)
       aaaTGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTg                                                                                                                                                >  2:3166621/1‑150 (MQ=255)
          tGCGTAATGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGgcg                                                                                                                                             >  2:443106/1‑150 (MQ=255)
            cGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGcgcg                                                                                                                                           >  2:1249703/1‑150 (MQ=255)
             gTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGa                                                                                                                                          >  2:2898363/1‑150 (MQ=255)
              tAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGaa                                                                                                                                         >  1:1742646/1‑150 (MQ=255)
               aaaGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAAc                                                                                                                                        <  1:4214405/150‑1 (MQ=255)
               aaaGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGAGCTGATGAATTTATTCTGACGGTAGAAAAAGGTCCGATTGGCGGAATAACTTCACAGGGGATCGCCTATGGCGCGAAc                                                                                                                                        >  2:4761737/1‑150 (MQ=255)
                 aGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGt                                                                                                                                      <  2:2103420/150‑1 (MQ=255)
                    cgcgGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGaa                                                                                                                                   <  2:843436/150‑1 (MQ=255)
                     gcgGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAAt                                                                                                                                  <  1:92466/150‑1 (MQ=255)
                      cgGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATa                                                                                                                                 <  1:1435919/150‑1 (MQ=255)
                      cgGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATa                                                                                                                                 <  2:4765296/150‑1 (MQ=255)
                       ggTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATAc                                                                                                                                <  1:4057281/150‑1 (MQ=255)
                       ggTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATAc                                                                                                                                <  1:3654666/150‑1 (MQ=255)
                        gTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATAcc                                                                                                                               >  1:3183765/1‑150 (MQ=255)
                        gTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAGTAcc                                                                                                                               <  2:3340769/150‑1 (MQ=255)
                         tGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATAccc                                                                                                                              >  1:4997669/1‑150 (MQ=255)
                          ggggAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCg                                                                                                                             >  2:3254345/1‑150 (MQ=255)
                          ggggAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCg                                                                                                                             <  2:4659633/150‑1 (MQ=255)
                            ggAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTg                                                                                                                           >  1:2456680/1‑150 (MQ=255)
                              aaTGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGcc                                                                                                                         >  2:3035162/1‑150 (MQ=255)
                                     gCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTgg                                                                                                                  <  2:3197442/150‑1 (MQ=255)
                                      cGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCGTTGGCGCGAACGTGAATACCCGTGCCATTCTGGa                                                                                                                 >  1:2645165/1‑150 (MQ=255)
                                       gTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGat                                                                                                                >  1:3253035/1‑150 (MQ=255)
                                       gTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGat                                                                                                                >  1:1025159/1‑150 (MQ=255)
                                         cGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGTCTTTATTCTTACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGatat                                                                                                              <  2:3381734/150‑1 (MQ=255)
                                          ggTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATg                                                                                                             >  1:3784783/1‑150 (MQ=255)
                                                gCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTcc                                                                                                       <  1:33813/150‑1 (MQ=255)
                                                gCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTcc                                                                                                       >  1:3557514/1‑150 (MQ=255)
                                                 cTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTccc                                                                                                      <  1:975425/150‑1 (MQ=255)
                                                 cTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTccc                                                                                                      <  1:1329075/150‑1 (MQ=255)
                                                  tGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCa                                                                                                     <  2:1301757/150‑1 (MQ=255)
                                                   gACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAg                                                                                                    >  1:5083970/1‑150 (MQ=255)
                                                   gACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAg                                                                                                    <  2:1526053/150‑1 (MQ=255)
                                                      ggCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGttt                                                                                                 <  2:5181102/150‑1 (MQ=255)
                                                          ttGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGAtt                                                                                             <  1:3102660/150‑1 (MQ=255)
                                                          ttGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGAtt                                                                                             <  1:4926601/150‑1 (MQ=255)
                                                          ttGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGAtt                                                                                             <  2:2929949/150‑1 (MQ=255)
                                                            ggCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGAtttt                                                                                           >  2:1166617/1‑150 (MQ=255)
                                                             gCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGAttttt                                                                                          <  1:1524053/150‑1 (MQ=255)
                                                             gCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGAttttt                                                                                          <  2:1763829/150‑1 (MQ=255)
                                                             gCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGAttttt                                                                                          <  2:4076795/150‑1 (MQ=255)
                                                                 ggTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCa                                                                                      >  2:3113479/1‑150 (MQ=255)
                                                                 ggTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCa                                                                                      >  2:1661902/1‑150 (MQ=255)
                                                                 ggTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCa                                                                                      <  2:1227651/150‑1 (MQ=255)
                                                                   tCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACg                                                                                    <  1:5156352/150‑1 (MQ=255)
                                                                   tCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACg                                                                                    >  1:2844581/1‑150 (MQ=255)
                                                                        cGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGc                                                                               >  1:631743/1‑150 (MQ=255)
                                                                         gagaAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCg                                                                              <  1:210903/150‑1 (MQ=255)
                                                                         gagaAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCg                                                                              >  1:3527599/1‑150 (MQ=255)
                                                                            aagaagGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTc                                                                           >  1:4142103/1‑150 (MQ=255)
                                                                                agGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGa                                                                       >  2:724244/1‑150 (MQ=255)
                                                                                  gTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGAtg                                                                     <  1:3817224/150‑1 (MQ=255)
                                                                                     gtgCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGAtgttt                                                                  <  2:1929717/150‑1 (MQ=255)
                                                                                          gatgaCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTAt                                                             <  1:4464361/150‑1 (MQ=255)
                                                                                            tgaCTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTAttt                                                           <  1:1122928/150‑1 (MQ=255)
                                                                                                tttATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCATTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGt                                                       <  2:5304068/150‑1 (MQ=255)
                                                                                                 ttATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGtt                                                      <  1:4816211/150‑1 (MQ=255)
                                                                                                 ttATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGtt                                                      >  2:5092244/1‑150 (MQ=255)
                                                                                                 ttAGTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGtt                                                      <  2:3604286/150‑1 (MQ=255)
                                                                                                   aTTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGtttt                                                    <  2:4508252/150‑1 (MQ=255)
                                                                                                   aTTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGtttt                                                    <  2:5012699/150‑1 (MQ=255)
                                                                                                    ttCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTg                                                   <  2:1340670/150‑1 (MQ=255)
                                                                                                    ttCTGACGGTAGAAACAGGTAAGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACTTGAATAACAGTGCCATTCTGGTTATGACGTACCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTg                                                   <  1:1780011/150‑1 (MQ=255)
                                                                                                      cTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCt                                                 <  2:1619985/150‑1 (MQ=255)
                                                                                                      cTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCt                                                 <  1:505032/150‑1 (MQ=255)
                                                                                                      cTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCt                                                 >  2:432268/1‑150 (MQ=255)
                                                                                                       tGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTg                                                <  2:458427/150‑1 (MQ=255)
                                                                                                        gACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGa                                               >  2:2377219/1‑150 (MQ=255)
                                                                                                         aCGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGaa                                              <  2:2263585/150‑1 (MQ=255)
                                                                                                         aCGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGaa                                              <  1:3339570/150‑1 (MQ=255)
                                                                                                          cGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAg                                             <  1:4582055/150‑1 (MQ=255)
                                                                                                           ggTAGATACTGGTACGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGAGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTCTCACGGTGGCGGTCTGGTTGTTTGTTATTTGAGTTTTGCTGAAGt                                            <  1:3354899/150‑1 (MQ=255)
                                                                                                           ggTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGt                                            >  2:3039877/1‑150 (MQ=255)
                                                                                                            gTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTc                                           <  2:4244047/150‑1 (MQ=255)
                                                                                                            gTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTc                                           >  1:3084317/1‑150 (MQ=255)
                                                                                                            gTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTc                                           <  1:4687498/150‑1 (MQ=255)
                                                                                                              aGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGa                                         >  1:2000478/1‑150 (MQ=255)
                                                                                                              aGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGAGGTTTGTTATTTGAGTTTTGCTGAAGTCGa                                         >  2:2240899/1‑150 (MQ=255)
                                                                                                                aaaCAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGAcc                                       >  2:1430906/1‑150 (MQ=255)
                                                                                                                  aCAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAg                                     <  1:2384670/150‑1 (MQ=255)
                                                                                                                    aGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCa                                   >  1:5249505/1‑150 (MQ=255)
                                                                                                                       tCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACgg                                >  1:1608914/1‑150 (MQ=255)
                                                                                                                         cGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTa                              <  2:5136708/150‑1 (MQ=255)
                                                                                                                           aTTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAAc                            <  1:432268/150‑1 (MQ=255)
                                                                                                                           aTTGGCGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAAc                            <  1:689882/150‑1 (MQ=255)
                                                                                                                               gCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCg                        >  1:4291758/1‑150 (MQ=255)
                                                                                                                                cGGAATTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCgg                       >  1:2268539/1‑150 (MQ=255)
                                                                                                                                   aaTTACTTCACAGGGGATAGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGt                    <  2:1025159/150‑1 (MQ=255)
                                                                                                                                     ttACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGc                  >  1:984654/1‑150 (MQ=255)
                                                                                                                                     ttACTTCACAGGGGATCGCCCTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGc                  >  2:3458266/1‑150 (MQ=255)
                                                                                                                                        cTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATa               >  1:2012431/1‑150 (MQ=255)
                                                                                                                                         ttCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATaa              <  1:3705724/150‑1 (MQ=255)
                                                                                                                                           cacaGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACATCGGCGTGCATAAAt            <  1:4344729/150‑1 (MQ=255)
                                                                                                                                           cacaGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACATCGGCGTGCATAAAt            <  1:1776481/150‑1 (MQ=255)
                                                                                                                                                  gATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATgg     <  1:1002010/150‑1 (MQ=255)
                                                                                                                                                   aTCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGt    <  2:1479351/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTATTCGAAATGCGTAAAGGCGCGGTGGGGAATGTCGGCGTCGGTATTGCTGACGGCATTGGCCTGGTCGCCCGAGAAGAAGGTTGTGCTGATGACTTTATTCTGACGGTAGAAACAGGTCCGATTGGCGGAATTACTTCACAGGGGATCGCCTTTGGCGCGAACGTGAATACCCGTGCCATTCTGGATATGACGTCCCAGTTTGATTTTTATCACGGTGGCGGTCTGGATGTTTGTTATTTGAGTTTTGCTGAAGTCGACCAGCACGGTAACGTCGGCGTGCATAAATTCAATGGTAA  >  CP009273/1770729‑1771027

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: