Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,738,088 T→A 19.2% 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,0880TA19.2% 741.7 / 166.9 307F134Y (TTC→TAC) yiaM2,3‑diketo‑L‑gulonate TRAP transporter small permease protein
Reads supporting (aligned to +/- strand):  ref base T (125/122);  new base A (32/27);  total (158/149)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 319 total aligned reads displayed.

ATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG  >  CP009273/3737941‑3738235
                                                                                                                                                   |                                                                                                                                                   
  tCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCt                                                                                                                                                   >  2:1900185/1‑148 (MQ=255)
  tCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACt                                                                                                                                                   <  2:1824652/148‑1 (MQ=255)
    cTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCttt                                                                                                                                                 <  1:1001030/148‑1 (MQ=255)
    cTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCttt                                                                                                                                                 <  2:1595704/148‑1 (MQ=255)
      tAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGa                                                                                                                                               <  2:233087/148‑1 (MQ=255)
      tAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGa                                                                                                                                               <  2:2242252/148‑1 (MQ=255)
       aaTCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGaa                                                                                                                                              >  2:1194683/1‑148 (MQ=255)
         tCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACt                                                                                                                                            >  2:1508902/1‑148 (MQ=255)
          cTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTg                                                                                                                                           >  2:2676720/1‑148 (MQ=255)
          cTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTg                                                                                                                                           >  1:2845867/1‑148 (MQ=255)
            tatttattTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCg                                                                                                                                         >  1:1624715/1‑148 (MQ=255)
            tatttattTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCg                                                                                                                                         >  1:791544/1‑148 (MQ=255)
               ttatttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCa                                                                                                                                      >  1:905092/1‑148 (MQ=255)
                   ttGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTta                                                                                                                                  >  1:147427/1‑148 (MQ=255)
                    tgtgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTtat                                                                                                                                 <  2:2610032/148‑1 (MQ=255)
                      tgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTatat                                                                                                                               <  1:2037611/148‑1 (MQ=255)
                      tgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTatat                                                                                                                               <  2:2888329/148‑1 (MQ=255)
                      tgGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTatat                                                                                                                               <  1:2745346/148‑1 (MQ=255)
                        gcgcACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCa                                                                                                                             >  2:775419/1‑148 (MQ=255)
                           cACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACt                                                                                                                          <  1:3334620/148‑1 (MQ=255)
                            aCTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtc                                                                                                                         <  2:2450946/148‑1 (MQ=255)
                            aCTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtc                                                                                                                         <  2:3152924/148‑1 (MQ=255)
                              tGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtcat                                                                                                                       <  2:2845867/148‑1 (MQ=255)
                               ggCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtcatc                                                                                                                      <  1:3583691/148‑1 (MQ=255)
                                gCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtcatca                                                                                                                     >  2:2013160/1‑148 (MQ=255)
                                gCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACtcatca                                                                                                                     >  2:3260588/1‑148 (MQ=255)
                                    ggggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCAcgcg                                                                                                                 <  1:1575377/148‑1 (MQ=255)
                                    ggggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCAcgcg                                                                                                                 <  2:2705269/148‑1 (MQ=255)
                                    ggggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCAcgcg                                                                                                                 <  1:52100/148‑1 (MQ=255)
                                     gggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                <  1:1519141/148‑1 (MQ=255)
                                     gggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                <  2:1321665/148‑1 (MQ=255)
                                     gggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                <  1:3418596/148‑1 (MQ=255)
                                     gggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                >  2:387140/1‑148 (MQ=255)
                                     gggCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGt                                                                                                                >  2:579406/1‑148 (MQ=255)
                                         gcAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCa                                                                                                            <  1:3809027/148‑1 (MQ=255)
                                         gcAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCa                                                                                                            <  1:1920007/148‑1 (MQ=255)
                                            aCGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAAtt                                                                                                         >  2:1085335/1‑148 (MQ=255)
                                             cGCTAAAAACAATTCAGGACTGGAGCGATTAGTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATtc                                                                                                        <  2:2181184/148‑1 (MQ=255)
                                                 aaaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTctctg                                                                                                    <  1:2478883/148‑2 (MQ=255)
                                                 aaaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTa                                                                                                    <  2:2102915/148‑1 (MQ=255)
                                                 aaaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTa                                                                                                    >  1:1736155/1‑148 (MQ=255)
                                                 aaaaaCAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTa                                                                                                    >  1:110631/1‑148 (MQ=255)
                                                      cAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTc                                                                                               >  2:1829884/1‑148 (MQ=255)
                                                        aTTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGc                                                                                             >  2:258673/1‑148 (MQ=255)
                                                        aTTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGc                                                                                             >  2:517489/1‑148 (MQ=255)
                                                        aTTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGc                                                                                             >  2:2110607/1‑148 (MQ=255)
                                                        aTTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCAGCATGTCTGCCCACCAGCCTCTTTATCGAAGACTTTGAACTGCGTCATTTATATCAACTCATCACGAGTAGCAATTCTCTaaagacga                                                                                             >  2:2371595/1‑142 (MQ=255)
                                                         ttCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGcc                                                                                            >  2:484360/1‑148 (MQ=255)
                                                          tCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCa                                                                                           >  2:1690593/1‑148 (MQ=255)
                                                            aGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAAGTCTCTAACGTCGCCAcc                                                                                         >  2:3234823/1‑148 (MQ=255)
                                                             ggACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCg                                                                                        >  1:3065023/1‑148 (MQ=255)
                                                              gACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTCTCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGc                                                                                       >  2:2946375/1‑148 (MQ=255)
                                                                  ggAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAgg                                                                                   <  1:3272440/148‑1 (MQ=255)
                                                                   gAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGa                                                                                  >  2:2507426/1‑148 (MQ=255)
                                                                     gCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGc                                                                                <  1:2367804/148‑1 (MQ=255)
                                                                        attattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTa                                                                             <  2:661936/148‑1 (MQ=255)
                                                                         ttattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTaa                                                                            >  2:79359/1‑148 (MQ=255)
                                                                          tattCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAAt                                                                           <  2:1710130/148‑1 (MQ=255)
                                                                             tCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCAt                                                                        >  1:2189734/1‑148 (MQ=255)
                                                                               aCCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATgg                                                                      >  2:370999/1‑148 (MQ=255)
                                                                                ccGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGc                                                                     >  1:2904544/1‑148 (MQ=255)
                                                                                    tCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgtg                                                                 >  2:140957/1‑148 (MQ=255)
                                                                                    tCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCtgtg                                                                 >  2:1362648/1‑148 (MQ=255)
                                                                                       tCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTg                                                              <  1:1180666/148‑1 (MQ=255)
                                                                                       tCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTg                                                              >  2:3342751/1‑148 (MQ=255)
                                                                                         ggTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAt                                                            <  2:2384849/148‑1 (MQ=255)
                                                                                          gTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAtt                                                           >  1:1080292/1‑148 (MQ=255)
                                                                                           tttACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAttt                                                          >  2:940220/1‑148 (MQ=255)
                                                                                            ttACCCCTCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAtttt                                                         >  2:2984406/1‑148 (MQ=255)
                                                                                            ttACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGAtttt                                                         <  2:433670/148‑1 (MQ=255)
                                                                                                ccATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTg                                                     >  2:1894984/1‑148 (MQ=255)
                                                                                                      gCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgt                                               <  1:2720708/148‑1 (MQ=255)
                                                                                                       ccTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtc                                              <  1:2274184/148‑1 (MQ=255)
                                                                                                       ccTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGctgtc                                              >  2:3499418/1‑148 (MQ=255)
                                                                                                           aTGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGtt                                          >  2:2803215/1‑148 (MQ=255)
                                                                                                            tGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTg                                         <  2:359789/148‑1 (MQ=255)
                                                                                                             gTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTgg                                        >  1:1676772/1‑148 (MQ=255)
                                                                                                                tGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggtg                                     <  2:3480380/148‑1 (MQ=255)
                                                                                                                  ccGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGggtggt                                   <  2:754092/148‑1 (MQ=255)
                                                                                                                      ttGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTg                               <  2:327321/148‑1 (MQ=255)
                                                                                                                      ttGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTg                               <  1:2933261/148‑1 (MQ=255)
                                                                                                                       tGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGc                              >  2:2651097/1‑148 (MQ=255)
                                                                                                                       tGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGc                              <  1:3474412/148‑1 (MQ=255)
                                                                                                                       tGTCTGCCCACCAGCCTCGTTATCGCATACTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGc                              <  2:1056796/148‑1 (MQ=255)
                                                                                                                        gTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCt                             <  1:3260588/148‑1 (MQ=255)
                                                                                                                             cccACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCgg                        >  2:2312724/1‑148 (MQ=255)
                                                                                                                                accaGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGttt                     >  1:805256/1‑148 (MQ=255)
                                                                                                                                     ccTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCGGTTGGGTGGTATTGCTATCGGTTTGCCTa                >  1:2584831/1‑148 (MQ=255)
                                                                                                                                        cGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTg             >  2:374042/1‑148 (MQ=255)
                                                                                                                                           tATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCt          >  2:889649/1‑148 (MQ=255)
                                                                                                                                           tATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGGATCAATTCTCTAACGTCGACACCGCATGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCt          <  2:1525698/148‑1 (MQ=255)
                                                                                                                                             tCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTgg        >  1:825296/1‑148 (MQ=255)
                                                                                                                                             tCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTgg        >  2:1278425/1‑148 (MQ=255)
                                                                                                                                              cGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTggg       >  2:455630/1‑148 (MQ=255)
                                                                                                                                              cGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTggg       >  1:2512187/1‑148 (MQ=255)
                                                                                                                                                 aTTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCAc    >  1:1305498/1‑148 (MQ=255)
                                                                                                                                                   tCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtg  >  1:3129137/1‑148 (MQ=255)
                                                                                                                                                   tCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtg  >  1:727749/1‑148 (MQ=255)
                                                                                                                                                   tCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtg  >  1:954125/1‑148 (MQ=255)
                                                                                                                                                   tCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACtg  >  1:1421091/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
ATTCCTTAATCTTATTTATTTGTGGCGCACTGGCCTGGGGCGCAACGCTAAAAACAATTCAGGACTGGAGCGATTATTCACCGATCCTCGGTTTACCCATCGGCCTGATGTATGCCGCTTGTCTGCCCACCAGCCTCGTTATCGCATTCTTTGAACTGCGTCATTTATATCAACTCATCACGCGTAGCAATTCTCTAACGTCGCCACCGCAAGGAGCTTAATCATGGCTGTGCTGATTTTTCTGGGCTGTCTGTTGGGTGGTATTGCTATCGGTTTGCCTATTGCCTGGGCACTG  >  CP009273/3737941‑3738235

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: