Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,527,522 G→T 7.3% R449L (CGT→CTT)  pck → phosphoenolpyruvate carboxykinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,527,5220GT 7.3% 1141.7 / 47.1 299R449L (CGT→CTT) pckphosphoenolpyruvate carboxykinase
Reads supporting (aligned to +/- strand):  ref base G (137/140);  new base T (4/18);  total (141/158)
Fisher's exact test for biased strand distribution p-value = 6.54e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 307 total aligned reads displayed.

GAACCGACGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAAC  >  CP009273/3527374‑3527670
                                                                                                                                                    |                                                                                                                                                    
gaaccgaCGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACg                                                                                                                                                      >  1:2716919/1‑149 (MQ=255)
    cgaCGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTAtc                                                                                                                                                  >  2:3684795/1‑149 (MQ=255)
      aCGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGtctctc                                                                                                                                                >  2:595976/1‑149 (MQ=255)
      aCGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTAtctc                                                                                                                                                <  1:3290204/149‑1 (MQ=255)
      aCGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTAtctc                                                                                                                                                >  1:4448821/1‑149 (MQ=255)
        gCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGa                                                                                                                                              >  2:4601764/1‑149 (MQ=255)
         ccAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGAt                                                                                                                                             <  1:4289535/149‑1 (MQ=255)
         ccAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGAt                                                                                                                                             <  2:4290250/149‑1 (MQ=255)
         ccAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGAt                                                                                                                                             >  2:1920451/1‑149 (MQ=255)
         ccAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGAt                                                                                                                                             >  2:146848/1‑149 (MQ=255)
         ccAACCTTCTCCGCGTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGAt                                                                                                                                             <  2:361647/149‑1 (MQ=255)
              cTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAg                                                                                                                                        <  2:4308384/149‑1 (MQ=255)
              cTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAg                                                                                                                                        >  1:1196726/1‑149 (MQ=255)
               ttCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGa                                                                                                                                       <  1:2119070/149‑1 (MQ=255)
                tctcCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGAt                                                                                                                                      >  2:3305698/1‑149 (MQ=255)
                  tcCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATAc                                                                                                                                    >  1:3952691/1‑149 (MQ=255)
                   ccGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATAcc                                                                                                                                   <  2:3746957/149‑1 (MQ=255)
                   ccGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACGGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATAcc                                                                                                                                   <  2:3040726/149‑1 (MQ=255)
                       ttgcttCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCcgcg                                                                                                                               >  1:1634001/1‑149 (MQ=255)
                          cttCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCa                                                                                                                            <  1:757784/149‑1 (MQ=255)
                             cGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTa                                                                                                                         <  1:1625774/149‑1 (MQ=255)
                             cGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGAGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATAACCGCGCCATTa                                                                                                                         >  2:1430300/1‑149 (MQ=255)
                                 gcgGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGa                                                                                                                     <  1:2992154/149‑1 (MQ=255)
                                  cgGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGAc                                                                                                                    <  2:1196630/149‑1 (MQ=255)
                                    gCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGc                                                                                                                  >  2:1952432/1‑149 (MQ=255)
                                     cATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGcc                                                                                                                 <  2:1343143/149‑1 (MQ=255)
                                       ttCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCAt                                                                                                               >  1:4043883/1‑149 (MQ=255)
                                          cTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACTTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCt                                                                                                            >  1:3933683/1‑149 (MQ=255)
                                          cTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACTTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCt                                                                                                            <  2:2742047/149‑1 (MQ=255)
                                          cTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCt                                                                                                            >  1:3022984/1‑149 (MQ=255)
                                          cTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCt                                                                                                            <  2:4069407/149‑1 (MQ=255)
                                            gTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCa                                                                                                          <  2:882070/149‑1 (MQ=255)
                                             tCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCaa                                                                                                         <  2:4158977/149‑1 (MQ=255)
                                              cGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAAc                                                                                                        <  1:3194372/149‑1 (MQ=255)
                                                 tGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGt                                                                                                     <  2:3968577/149‑1 (MQ=255)
                                                 tGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGt                                                                                                     <  2:1667434/149‑1 (MQ=255)
                                                 tGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGt                                                                                                     <  2:2495481/149‑1 (MQ=255)
                                                   cACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTc                                                                                                   >  1:2539613/1‑149 (MQ=255)
                                                    aCCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCg                                                                                                  >  2:723895/1‑149 (MQ=255)
                                                      ccGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCt                                                                                                <  1:3096207/149‑1 (MQ=255)
                                                       cGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTg                                                                                               >  2:2198700/1‑149 (MQ=255)
                                                       cGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTg                                                                                               <  2:1071576/149‑1 (MQ=255)
                                                        gACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGTCAATCGTATCTCGATTAATGATTCCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTgg                                                                                              <  1:3078220/149‑1 (MQ=255)
                                                         aCTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACTTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGa                                                                                             <  2:3875134/149‑1 (MQ=255)
                                                         aCTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGa                                                                                             >  1:1769697/1‑149 (MQ=255)
                                                             aGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAAt                                                                                         <  1:1918691/149‑1 (MQ=255)
                                                             aGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAAt                                                                                         <  1:1368098/149‑1 (MQ=255)
                                                             aGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAAt                                                                                         <  1:3438046/149‑1 (MQ=255)
                                                             aGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCAGCCTCAACGGTTCGCTGGATAAt                                                                                         >  1:4396038/1‑149 (MQ=255)
                                                               tACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGc                                                                                       <  2:3143464/149‑1 (MQ=255)
                                                                aCGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCa                                                                                      <  2:907562/149‑1 (MQ=255)
                                                                   cAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACTTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGaa                                                                                   <  2:2862472/149‑1 (MQ=255)
                                                                     gAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAAc                                                                                 <  2:4473394/149‑1 (MQ=255)
                                                                      aaGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACTTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAAcc                                                                                <  2:2341385/149‑1 (MQ=255)
                                                                       aGTGCTGGAGAAACGTATGCAGGCGGCGGGCGCGCAGGCTGATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTAGCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGCATAATGCAGAAACCt                                                                               >  2:4677882/1‑149 (MQ=255)
                                                                         tGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTc                                                                             >  1:4534575/1‑149 (MQ=255)
                                                                            tGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCAct                                                                          >  1:2653952/1‑149 (MQ=255)
                                                                             ggTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCActc                                                                         <  2:4193355/149‑1 (MQ=255)
                                                                               tGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTAACTCTg                                                                       >  2:3255789/1‑149 (MQ=255)
                                                                                gAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGc                                                                      <  2:4448789/149‑1 (MQ=255)
                                                                                 aaaCGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGcc                                                                     >  1:3261418/1‑149 (MQ=255)
                                                                                 aaaCGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGcc                                                                     >  2:2135643/1‑149 (MQ=255)
                                                                                 aaaCGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGcc                                                                     <  2:2602184/149‑1 (MQ=255)
                                                                                  aaCGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCg                                                                    >  1:383140/1‑149 (MQ=255)
                                                                                  aaCGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCg                                                                    >  1:4211599/1‑149 (MQ=255)
                                                                                    cGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGAt                                                                  <  2:2872527/149‑1 (MQ=255)
                                                                                            ggcggcGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAAcc                                                          >  1:712271/1‑149 (MQ=255)
                                                                                               ggcgggcgCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTgg                                                       <  1:4670703/149‑1 (MQ=255)
                                                                                               ggcgggcgCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTgg                                                       <  1:149638/149‑1 (MQ=255)
                                                                                                gcgggcgCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGc                                                      <  1:2010168/149‑1 (MQ=255)
                                                                                                     cgcgcAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATcc                                                 >  2:3801336/1‑149 (MQ=255)
                                                                                                        gcAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCaa                                              <  2:2339188/149‑1 (MQ=255)
                                                                                                          aGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAAcc                                            <  2:4172351/149‑1 (MQ=255)
                                                                                                           ggCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCg                                           >  1:1852418/1‑149 (MQ=255)
                                                                                                            gCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGa                                          >  1:353203/1‑149 (MQ=255)
                                                                                                                  cTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGcc                                    <  1:555924/149‑1 (MQ=255)
                                                                                                                    ggTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCgg                                  >  1:2268640/1‑149 (MQ=255)
                                                                                                                     gTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCggg                                 >  1:991134/1‑149 (MQ=255)
                                                                                                                      ttAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGc                                >  1:3688981/1‑149 (MQ=255)
                                                                                                                        aaCACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGt                              >  1:338030/1‑149 (MQ=255)
                                                                                                                        aaCACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGt                              >  1:2761474/1‑149 (MQ=255)
                                                                                                                           acTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGa                           <  1:1709027/149‑1 (MQ=255)
                                                                                                                             tggctggTACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGaca                         <  1:3488525/149‑1 (MQ=255)
                                                                                                                             tggctggAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGaca                         <  1:1225717/149‑1 (MQ=255)
                                                                                                                              ggctggAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGacac                        >  2:1641500/1‑149 (MQ=255)
                                                                                                                               gctggAACGGCACTGGCAAACGTCTCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACg                       >  2:1378409/1‑149 (MQ=255)
                                                                                                                                ctggAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGa                      >  2:2123336/1‑149 (MQ=255)
                                                                                                                                ctggAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGa                      >  1:1406882/1‑149 (MQ=255)
                                                                                                                                   gAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGa                   >  2:849185/1‑149 (MQ=255)
                                                                                                                                     aCGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGAtt                 <  2:3508012/149‑1 (MQ=255)
                                                                                                                                     aCGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGAtt                 >  2:4271467/1‑149 (MQ=255)
                                                                                                                                     aCGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGAtt                 <  2:4011780/149‑1 (MQ=255)
                                                                                                                                       ggCACTGGCAAACTTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATtct               <  2:1269575/149‑1 (MQ=255)
                                                                                                                                       ggCACTGGCAAACTTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATtct               <  2:1584585/149‑1 (MQ=255)
                                                                                                                                          aCTGGCAAACTTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGa            >  1:1012952/1‑149 (MQ=255)
                                                                                                                                           cTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGAt           <  2:1833745/149‑1 (MQ=255)
                                                                                                                                            tGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATc          <  2:701810/149‑1 (MQ=255)
                                                                                                                                             ggCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATcc         >  1:3595526/1‑149 (MQ=255)
                                                                                                                                                    gTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAac  >  1:1838431/1‑149 (MQ=255)
                                                                                                                                                    gTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAac  >  1:4479760/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAACCGACGCCAACCTTCTCCGCTTGCTTCGGCGCGGCATTCCTGTCGCTGCACCCGACTCAGTACGCAGAAGTGCTGGTGAAACGTATGCAGGCGGCGGGCGCGCAGGCTTATCTGGTTAACACTGGCTGGAACGGCACTGGCAAACGTATCTCGATTAAAGATACCCGCGCCATTATCGACGCCATCCTCAACGGTTCGCTGGATAATGCAGAAACCTTCACTCTGCCGATGTTTAACCTGGCGATCCCAACCGAACTGCCGGGCGTAGACACGAAGATTCTCGATCCGCGTAAC  >  CP009273/3527374‑3527670

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: