Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,574,820 C→A 100% A609E (GCG→GAG)  tktB → transketolase 2, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,574,8200CA98.2% 681.9 / 3.1 220A609E (GCG→GAG) tktBtransketolase 2, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (92/124);  total (93/127)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 220 total aligned reads displayed.

GAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAA  >  CP009273/2574675‑2574965
                                                                                                                                                 |                                                                                                                                                 
gAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGt                                                                                                                                                <  2:2687199/149‑1 (MQ=255)
gAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGt                                                                                                                                                <  2:3514304/149‑1 (MQ=255)
  aaTTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTgg                                                                                                                                              >  1:3904325/1‑149 (MQ=255)
    ttACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGaa                                                                                                                                            >  1:515915/1‑149 (MQ=255)
    ttACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGaa                                                                                                                                            <  2:4112932/149‑1 (MQ=255)
       cccTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCa                                                                                                                                         <  1:1175563/149‑1 (MQ=255)
       cccTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCa                                                                                                                                         >  1:947795/1‑149 (MQ=255)
        ccTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAg                                                                                                                                        >  2:2618537/1‑149 (MQ=255)
           gCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTa                                                                                                                                     >  2:4068999/1‑149 (MQ=255)
            cAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTAt                                                                                                                                    <  1:2270699/149‑1 (MQ=255)
             aaGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTAtt                                                                                                                                   <  2:1246431/149‑1 (MQ=255)
               gCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGc                                                                                                                                 <  1:3639386/149‑1 (MQ=255)
                  gCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGa                                                                                                                              <  2:3989685/149‑1 (MQ=255)
                    agagAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGAtt                                                                                                                            <  1:3877686/149‑1 (MQ=255)
                    agagAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGAtt                                                                                                                            >  2:1916785/1‑149 (MQ=255)
                     gagaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTa                                                                                                                           >  2:3191746/1‑149 (MQ=255)
                      agaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTAc                                                                                                                          <  1:3984121/149‑1 (MQ=255)
                      agaAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTAc                                                                                                                          >  2:3102179/1‑149 (MQ=255)
                          aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGt                                                                                                                      >  1:539167/1‑149 (MQ=255)
                          aTTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGt                                                                                                                      >  2:3412065/1‑149 (MQ=255)
                           ttAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTa                                                                                                                     >  1:3029900/1‑149 (MQ=255)
                              gCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACaa                                                                                                                  <  1:576952/149‑1 (MQ=255)
                                 ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                                                                                                               <  1:1647198/149‑1 (MQ=255)
                                 ggAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTa                                                                                                               <  1:336904/149‑1 (MQ=255)
                                  gagaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTAt                                                                                                              >  1:274239/1‑149 (MQ=255)
                                    gaAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGt                                                                                                            <  2:755351/149‑1 (MQ=255)
                                      aGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTg                                                                                                          >  1:646972/1‑149 (MQ=255)
                                       ggTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTgg                                                                                                         >  1:2425123/1‑149 (MQ=255)
                                         tCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc                                                                                                       <  1:450692/149‑1 (MQ=255)
                                         tCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTc                                                                                                       <  1:3278267/149‑1 (MQ=255)
                                           gCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTg                                                                                                     >  1:2095174/1‑149 (MQ=255)
                                             aaTGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGaa                                                                                                   >  2:280382/1‑149 (MQ=255)
                                               tGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTTCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTAATGGTACAAGTATGTTGGTCTGAAAg                                                                                                 >  1:3442464/1‑149 (MQ=255)
                                                  aCGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAgggg                                                                                              <  1:3498221/149‑1 (MQ=255)
                                                      gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGCGGAATATCGGGAGTCGGTGTTGCCGTCTAACGGTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                          <  2:2039398/149‑1 (MQ=255)
                                                      gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                          <  2:1082022/149‑1 (MQ=255)
                                                      gTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAt                                                                                          <  1:1489269/149‑1 (MQ=255)
                                                       tAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAAtt                                                                                         <  2:3835198/149‑1 (MQ=255)
                                                         gggtCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGt                                                                                       <  1:3673283/146‑1 (MQ=255)
                                                         gTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGt                                                                                       <  1:498487/149‑1 (MQ=255)
                                                          tttCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTc                                                                                      <  2:2854527/149‑1 (MQ=255)
                                                              ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGa                                                                                  <  1:2304099/149‑1 (MQ=255)
                                                              ccTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGa                                                                                  >  1:1022643/1‑149 (MQ=255)
                                                                 gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGa                                                                               <  1:651369/149‑1 (MQ=255)
                                                                 gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGa                                                                               <  1:1386124/149‑1 (MQ=255)
                                                                 gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGa                                                                               <  1:3695826/149‑1 (MQ=255)
                                                                 gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGa                                                                               <  2:3441779/149‑1 (MQ=255)
                                                                 gCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTGCAAACGTTTCGGCTCGCGTGGAGGTGGAAGCAGGTATTTCCGAATACATGTACAAGAATGTTGGTAAGAAAGGGGCAATTGTCGGGAAGa                                                                               >  1:1071378/1‑149 (MQ=255)
                                                                    cTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACgg                                                                            >  2:1514804/1‑149 (MQ=255)
                                                                        aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTa                                                                        >  2:536208/1‑149 (MQ=255)
                                                                        aCCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTa                                                                        >  2:3345389/1‑149 (MQ=255)
                                                                          cGATATTTTCGACGCCCAGGATGAGGTATTTCGGGTGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTTTTGCCGATTTGTGGTACAAGTATGTGGGTCTGAAAGGGGCAATTGTCGGTATGACGGGTTACg                                                                      <  2:894138/149‑1 (MQ=255)
                                                                          cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                                      <  2:2954971/149‑1 (MQ=255)
                                                                          cGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACg                                                                      >  2:1396783/1‑149 (MQ=255)
                                                                           gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgg                                                                     <  2:2461850/149‑1 (MQ=255)
                                                                           gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgg                                                                     >  1:21950/1‑149 (MQ=255)
                                                                           gATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACgg                                                                     <  2:1044889/149‑1 (MQ=255)
                                                                               ttttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGaa                                                                 <  1:942112/149‑1 (MQ=255)
                                                                               ttttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGaa                                                                 >  1:3295791/1‑149 (MQ=255)
                                                                                tttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAAt                                                                >  1:4017494/1‑149 (MQ=255)
                                                                                 ttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATc                                                               <  1:1284117/149‑1 (MQ=255)
                                                                                 ttCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATc                                                               <  1:2460419/149‑1 (MQ=255)
                                                                                    gACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGc                                                            >  2:2733381/1‑149 (MQ=255)
                                                                                    gACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGc                                                            <  1:708076/149‑1 (MQ=255)
                                                                                     aCGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCt                                                           <  1:1916772/149‑1 (MQ=255)
                                                                                      cGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTTGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTc                                                          <  1:1557276/149‑1 (MQ=255)
                                                                                      cGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTc                                                          <  1:2000025/149‑1 (MQ=255)
                                                                                      cGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTc                                                          <  1:3665996/149‑1 (MQ=255)
                                                                                      cGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTc                                                          <  1:3974811/149‑1 (MQ=255)
                                                                                          cAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggc                                                      >  2:4042294/1‑149 (MQ=255)
                                                                                            ggATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggcgg                                                    >  2:3974294/1‑149 (MQ=255)
                                                                                            ggATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggcgg                                                    >  1:1999495/1‑149 (MQ=255)
                                                                                            ggATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCcggcgg                                                    <  1:2613526/149‑1 (MQ=255)
                                                                                             gATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGa                                                   <  1:3536406/149‑1 (MQ=255)
                                                                                               tGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATa                                                 >  1:172007/1‑149 (MQ=255)
                                                                                               tGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATa                                                 <  1:2740832/149‑1 (MQ=255)
                                                                                                gAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATaa                                                >  2:3325074/1‑149 (MQ=255)
                                                                                                    aaTATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTg                                            <  2:515936/149‑1 (MQ=255)
                                                                                                    aaTATCGGGAGTCGGTGTTGCCTTCTAACGGTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTg                                            <  2:515890/149‑1 (MQ=255)
                                                                                                        tCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTcc                                        <  1:292278/149‑1 (MQ=255)
                                                                                                           ggAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGt                                     <  2:908048/149‑1 (MQ=255)
                                                                                                              gTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                                  <  2:1633237/149‑1 (MQ=255)
                                                                                                              gTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCt                                  <  2:4013984/149‑1 (MQ=255)
                                                                                                                  gtgtTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACACGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTgg                              >  1:50581/1‑149 (MQ=255)
                                                                                                                  gtgtTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTgg                              >  1:3093025/1‑149 (MQ=255)
                                                                                                                    gtTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCt                            <  2:274113/149‑1 (MQ=255)
                                                                                                                    gtTGCCTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCt                            <  1:1589479/149‑1 (MQ=255)
                                                                                                                        ccTTCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAc                        >  1:888960/1‑149 (MQ=255)
                                                                                                                          ttCTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTAAGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccg                      >  2:1215664/1‑149 (MQ=255)
                                                                                                                            cTAACGTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTAccgcc                    >  2:185935/1‑149 (MQ=255)
                                                                                                                                 gTTGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGaa               <  1:3126126/149‑1 (MQ=255)
                                                                                                                                  ttGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAat              >  2:564640/1‑149 (MQ=255)
                                                                                                                                   tGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAata             <  1:2859670/149‑1 (MQ=255)
                                                                                                                                   tGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAata             <  1:800861/149‑1 (MQ=255)
                                                                                                                                   tGCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAata             <  1:1144954/149‑1 (MQ=255)
                                                                                                                                    gCGGCTCGCGTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAatat            <  2:233949/149‑1 (MQ=255)
                                                                                                                                           gcgTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCa     >  1:2839148/1‑149 (MQ=255)
                                                                                                                                            cgTGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCaa    <  1:2536861/149‑1 (MQ=255)
                                                                                                                                             gTGGAGGTGGAATCAGTTATTGCCGATTACTTGTAAAAGTATGTTGGTCTTAAAGGGGCAATTGTCGTGATGACGGGTTACGGTGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCaaa   <  2:3803550/149‑1 (MQ=255)
                                                                                                                                              tGGAGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCaaaa  <  2:21950/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                 
GAAATTACCCTGCAAGCGGCAGAGAAATTAGCAGGAGAAGGTCGCAATGTACGCGTAGTTTCCCTGCCCTCGACCGATATTTTCGACGCCCAGGATGAGGAATATCGGGAGTCGGTGTTGCCTTCTAACGTTGCGGCTCGCGTGGCGGTGGAAGCAGGTATTGCCGATTACTGGTACAAGTATGTTGGTCTGAAAGGGGCAATTGTCGGGATGACGGGTTACGGGGAATCTGCTCCGGCGGATAAGCTGTTCCCGTTCTTTGGCTTTACCGCCGAGAATATTGTGGCAAAA  >  CP009273/2574675‑2574965

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: