Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,738,088 T→A 17.9% F134Y (TTC→TAC)  yiaM → 2,3‑diketo‑L‑gulonate TRAP transporter small permease protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,738,0880TA17.9% 610.8 / 118.8 245F134Y (TTC→TAC) yiaM2,3‑diketo‑L‑gulonate TRAP transporter small permease protein
Reads supporting (aligned to +/- strand):  ref base T (100/101);  new base A (15/29);  total (115/130)
Fisher's exact test for biased strand distribution p-value = 6.76e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 250 total aligned reads displayed.

CATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGT  >  CP009273/3737940‑3738236
                                                                                                                                                    |                                                                                                                                                    
cATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAtt                                                                                                                                                      <  2:3929762/149‑1 (MQ=255)
cATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAtt                                                                                                                                                      >  1:142246/1‑149 (MQ=255)
 aTTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTc                                                                                                                                                     <  1:1650806/149‑1 (MQ=255)
  ttCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCt                                                                                                                                                    <  2:1719160/149‑1 (MQ=255)
   tCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCtt                                                                                                                                                   <  1:2679435/149‑1 (MQ=255)
     cTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTg                                                                                                                                                 >  1:798873/1‑149 (MQ=255)
      ttAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGa                                                                                                                                                >  1:3995664/1‑149 (MQ=255)
       tAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGaa                                                                                                                                               <  1:3391171/149‑1 (MQ=255)
         aTCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACt                                                                                                                                             <  1:146463/149‑1 (MQ=255)
         aTCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCACTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACt                                                                                                                                             >  1:2742311/1‑149 (MQ=255)
           cTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGc                                                                                                                                           <  1:1919477/149‑1 (MQ=255)
              atttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTc                                                                                                                                        <  1:2384968/149‑1 (MQ=255)
              atttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTc                                                                                                                                        >  1:2780606/1‑149 (MQ=255)
               tttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCa                                                                                                                                       >  2:3090381/1‑149 (MQ=255)
                ttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCAt                                                                                                                                      <  2:404439/149‑1 (MQ=255)
                ttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCAt                                                                                                                                      <  1:683969/149‑1 (MQ=255)
                  atttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCAttt                                                                                                                                    >  2:2241092/1‑149 (MQ=255)
                   tttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTta                                                                                                                                   <  1:2253938/149‑1 (MQ=255)
                   tttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTta                                                                                                                                   >  2:1327645/1‑149 (MQ=255)
                      gtgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTatat                                                                                                                                <  1:2618972/149‑1 (MQ=255)
                            cACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtc                                                                                                                          <  1:3372003/149‑1 (MQ=255)
                              cTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCa                                                                                                                                       <  1:1132398/134‑1 (MQ=255)
                              cTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtcat                                                                                                                        >  1:2288120/1‑149 (MQ=255)
                                ggCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACtcatca                                                                                                                      >  2:1062269/1‑149 (MQ=255)
                                   cTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCAcgc                                                                                                                   <  1:3091039/149‑1 (MQ=255)
                                     ggggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                 <  2:2178789/149‑1 (MQ=255)
                                     ggggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                 >  2:460627/1‑149 (MQ=255)
                                      gggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTa                                                                                                                >  2:901023/1‑149 (MQ=255)
                                           cAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAAt                                                                                                           >  1:2318126/1‑149 (MQ=255)
                                             aCGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtc                                                                                                         >  1:2329706/1‑149 (MQ=255)
                                             aCGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtc                                                                                                         <  2:4067435/149‑1 (MQ=255)
                                              cGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtct                                                                                                        <  2:2278559/149‑1 (MQ=255)
                                               gCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtctc                                                                                                       <  2:684251/149‑1 (MQ=255)
                                                   aaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAAc                                                                                                   >  1:150989/1‑149 (MQ=255)
                                                     aaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGt                                                                                                 >  2:2961045/1‑149 (MQ=255)
                                                      aCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTc                                                                                                <  1:2561603/149‑1 (MQ=255)
                                                      aCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTc                                                                                                >  2:1117373/1‑149 (MQ=255)
                                                          ttCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCa                                                                                            >  2:2532774/1‑149 (MQ=255)
                                                          ttCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCa                                                                                            >  1:1454611/1‑149 (MQ=255)
                                                              ggACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGc                                                                                        <  1:2612268/149‑1 (MQ=255)
                                                                 cTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAg                                                                                     <  1:3934944/149‑1 (MQ=255)
                                                                    gAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAg                                                                                  >  1:397474/1‑149 (MQ=255)
                                                                    gAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCAAGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAg                                                                                  >  1:2143101/1‑149 (MQ=255)
                                                                      gCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCt                                                                                <  1:1979572/149‑1 (MQ=255)
                                                                        gATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTa                                                                              >  1:2725561/1‑149 (MQ=255)
                                                                          ttattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAAt                                                                            <  2:1431512/149‑1 (MQ=255)
                                                                          ttattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAAt                                                                            <  1:901023/149‑1 (MQ=255)
                                                                          ttattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATAAACTCATCACGCGTAGCACTTCTCTAACGTCGCCACCGCAAGGAGCTTAAt                                                                            >  1:1064641/1‑149 (MQ=255)
                                                                           tattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATc                                                                           >  1:1341343/1‑149 (MQ=255)
                                                                           tattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATc                                                                           <  2:3517264/149‑1 (MQ=255)
                                                                           tattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATc                                                                           <  2:1029693/149‑1 (MQ=255)
                                                                           tattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATc                                                                           <  2:1029511/149‑1 (MQ=255)
                                                                             ttCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCAt                                                                         >  1:445150/1‑149 (MQ=255)
                                                                                aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGc                                                                      <  1:1150481/149‑1 (MQ=255)
                                                                                aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGc                                                                      >  1:2258609/1‑149 (MQ=255)
                                                                                aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCAGTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGc                                                                      >  2:3341714/1‑149 (MQ=255)
                                                                                   gATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgt                                                                   >  1:2526987/1‑149 (MQ=255)
                                                                                   gATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgt                                                                   <  1:3238056/149‑1 (MQ=255)
                                                                                     tCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGc                                                                 <  1:3559220/149‑1 (MQ=255)
                                                                                      ccTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCt                                                                >  1:1889183/1‑149 (MQ=255)
                                                                                        tCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGa                                                              >  2:1554679/1‑149 (MQ=255)
                                                                                        tCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGa                                                              >  1:2838839/1‑149 (MQ=255)
                                                                                          ggTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTACCGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAtt                                                            <  2:3476130/149‑1 (MQ=255)
                                                                                           gTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAttt                                                           >  2:3694/1‑149 (MQ=255)
                                                                                           gTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAttt                                                           >  2:3738/1‑149 (MQ=255)
                                                                                           gTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGGGCTGAttt                                                           >  1:3183727/1‑149 (MQ=255)
                                                                                           gTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAttt                                                           >  1:3842015/1‑149 (MQ=255)
                                                                                             ttACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAttttt                                                         <  1:3195045/149‑1 (MQ=255)
                                                                                                cccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTg                                                      >  2:1192595/1‑149 (MQ=255)
                                                                                                 ccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTgg                                                     >  1:347134/1‑149 (MQ=255)
                                                                                                   aTCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGc                                                   <  2:2873025/149‑1 (MQ=255)
                                                                                                   aTCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGc                                                   <  2:3211286/149‑1 (MQ=255)
                                                                                                     cGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctg                                                 <  2:677690/149‑1 (MQ=255)
                                                                                                       gCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtc                                               >  1:2592905/1‑149 (MQ=255)
                                                                                                        ccTGATGTATGCCGCGTGTCTGCCCACCAGCCTGGTTATCGCCTTCTTGGAACTGCTTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtct                                              <  1:299851/149‑1 (MQ=255)
                                                                                                         cTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtctg                                             >  2:839356/1‑149 (MQ=255)
                                                                                                          tGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtctgt                                            <  1:2198794/149‑1 (MQ=255)
                                                                                                           gATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGtt                                           <  2:1336911/149‑1 (MQ=255)
                                                                                                             tGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTgg                                         <  1:1715981/149‑1 (MQ=255)
                                                                                                              gTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTggg                                        >  1:3880719/1‑149 (MQ=255)
                                                                                                              gTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTggg                                        >  2:557214/1‑149 (MQ=255)
                                                                                                               tATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggt                                       <  2:2892903/149‑1 (MQ=255)
                                                                                                                aTGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggtg                                      <  1:901820/149‑1 (MQ=255)
                                                                                                                    cGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTAt                                  <  2:3269218/149‑1 (MQ=255)
                                                                                                                      cTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTg                                >  2:2037317/1‑149 (MQ=255)
                                                                                                                        tGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCt                              <  2:2463945/149‑1 (MQ=255)
                                                                                                                          tCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTAt                            <  2:2338815/149‑1 (MQ=255)
                                                                                                                               ccaccaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGtt                       <  2:1950400/149‑1 (MQ=255)
                                                                                                                                caccaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGttt                      <  1:1005359/149‑1 (MQ=255)
                                                                                                                                caccaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGttt                      >  1:525335/1‑149 (MQ=255)
                                                                                                                                   caGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGcc                   <  1:1085938/149‑1 (MQ=255)
                                                                                                                                    aGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCt                  <  2:3850167/149‑1 (MQ=255)
                                                                                                                                       cTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTAtt               <  1:1805902/149‑1 (MQ=255)
                                                                                                                                         cGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGc             <  1:931772/149‑1 (MQ=255)
                                                                                                                                         cGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGc             <  1:2232791/149‑1 (MQ=255)
                                                                                                                                            tATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTg          >  1:1070438/1‑149 (MQ=255)
                                                                                                                                              tCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTCATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTggg        >  2:2988479/1‑149 (MQ=255)
                                                                                                                                                gCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCa      >  2:148551/1‑149 (MQ=255)
                                                                                                                                                 cATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAc     <  2:2596687/149‑1 (MQ=255)
                                                                                                                                                    tCTTTGAACTGCGTCATTTATATCAACTCATCAAGCTTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTATGGGCTGTCTGTTGGGTGGTATTGATATCGGTTTGCCTATTGCCTGGGCACtgt  >  1:2956655/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTGT  >  CP009273/3737940‑3738236

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: