Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,242,509 C→T 100% L26F (CTC→TTC)  ubiC → chorismate‑‑pyruvate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,242,5090CT97.1% 904.0 / 12.1 271L26F (CTC→TTC) ubiCchorismate‑‑pyruvate lyase
Reads supporting (aligned to +/- strand):  ref base C (0/8);  new base T (118/145);  total (118/153)
Fisher's exact test for biased strand distribution p-value = 1.06e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-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 271 total aligned reads displayed.

TTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTA  >  CP009273/4242363‑4242654
                                                                                                                                                  |                                                                                                                                                 
tttGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTCTTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTc                                                                                                                                                 >  2:4285249/1‑149 (MQ=255)
 ttGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCg                                                                                                                                                >  1:3538461/1‑149 (MQ=255)
   gCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGAc                                                                                                                                              <  2:4577672/149‑1 (MQ=255)
   gCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGAc                                                                                                                                              <  2:3038203/149‑1 (MQ=255)
     tttACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTg                                                                                                                                            >  1:1265156/1‑149 (MQ=255)
         cATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTg                                                                                                                                        >  2:3673119/1‑149 (MQ=255)
          aTCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGt                                                                                                                                       >  1:424935/1‑149 (MQ=255)
           tctcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGtt                                                                                                                                      >  2:2050033/1‑149 (MQ=255)
           tctcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGtt                                                                                                                                      >  1:1598977/1‑149 (MQ=255)
           tctcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGtt                                                                                                                                      >  2:1256345/1‑149 (MQ=255)
             tcCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGc                                                                                                                                    >  2:3355405/1‑149 (MQ=255)
                 gTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTgga                                                                                                                                >  2:717801/1‑149 (MQ=255)
                  tAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCACCTGTTCGACTGGCTGTTGCTggag                                                                                                                               >  1:3786911/1‑149 (MQ=255)
                   aaaCACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggagg                                                                                                                              <  1:1574245/149‑1 (MQ=255)
                    aaCACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCACCTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTggagga                                                                                                                             >  2:2890509/1‑149 (MQ=255)
                     acacTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGAt                                                                                                                            >  1:4456474/1‑149 (MQ=255)
                     acacTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGAt                                                                                                                            <  2:3780268/149‑1 (MQ=255)
                       acTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTc                                                                                                                          >  1:3699582/1‑149 (MQ=255)
                         ttttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCa                                                                                                                        >  1:4765943/1‑149 (MQ=255)
                          tttCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCAt                                                                                                                       <  2:749890/149‑1 (MQ=255)
                            tCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGa                                                                                                                     <  1:4167461/149‑1 (MQ=255)
                            tCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGa                                                                                                                     >  2:3828277/1‑149 (MQ=255)
                                 gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaaa                                                                                                                <  2:4562021/149‑1 (MQ=255)
                                 gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaaa                                                                                                                <  2:4332267/149‑1 (MQ=255)
                                 gATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACaaaa                                                                                                                <  1:2191565/149‑1 (MQ=255)
                                  aTACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAAc                                                                                                               >  1:1864360/1‑149 (MQ=255)
                                   tACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACg                                                                                                              <  1:951001/149‑1 (MQ=255)
                                      aaTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGttt                                                                                                           <  1:2841197/149‑1 (MQ=255)
                                      aaTGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGttt                                                                                                           <  1:2366135/149‑1 (MQ=255)
                                        tGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGAGTCCATGACAAAACGTTTTg                                                                                                         >  1:1090350/1‑149 (MQ=255)
                                           cTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAc                                                                                                      <  1:4589466/149‑1 (MQ=255)
                                           cTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAc                                                                                                      <  1:4306601/149‑1 (MQ=255)
                                            tttACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCTCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAAca                                                                                                     <  1:2660719/149‑1 (MQ=255)
                                                 gTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAgg                                                                                                <  2:1295386/149‑1 (MQ=255)
                                                    aTGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGaa                                                                                             >  1:3031205/1‑149 (MQ=255)
                                                         aCGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACg                                                                                        <  2:3432936/149‑1 (MQ=255)
                                                         aCGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACg                                                                                        >  1:4518083/1‑149 (MQ=255)
                                                          cGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACgg                                                                                       <  2:3876689/149‑1 (MQ=255)
                                                          cGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACgg                                                                                       >  1:1973187/1‑149 (MQ=255)
                                                          cGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACgg                                                                                       <  1:2924160/149‑1 (MQ=255)
                                                           ggAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGt                                                                                      >  1:1192394/1‑149 (MQ=255)
                                                            gagaGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTa                                                                                     >  2:4309795/1‑149 (MQ=255)
                                                            gagaGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTa                                                                                     >  2:1385512/1‑149 (MQ=255)
                                                                gTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCg                                                                                 <  2:4086260/149‑1 (MQ=255)
                                                                    ggCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGAc                                                                             <  2:3735248/149‑1 (MQ=255)
                                                                    ggCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGAc                                                                             <  1:2951579/149‑1 (MQ=255)
                                                                      cATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACga                                                                           <  2:4415011/149‑1 (MQ=255)
                                                                        tGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatg                                                                         >  1:1874499/1‑149 (MQ=255)
                                                                         gTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatga                                                                        >  2:2324150/1‑149 (MQ=255)
                                                                          tCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatgat                                                                       <  1:4227583/149‑1 (MQ=255)
                                                                          tCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACgatgat                                                                       <  1:2099425/149‑1 (MQ=255)
                                                                           cacaCCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATc                                                                      <  1:4143054/149‑1 (MQ=255)
                                                                           cacaCCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATc                                                                      <  2:2119267/149‑1 (MQ=255)
                                                                              acCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCcgc                                                                   <  1:3374989/149‑1 (MQ=255)
                                                                              acCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCcgc                                                                   >  2:2521624/1‑149 (MQ=255)
                                                                               ccccGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCcgcg                                                                  >  1:709833/1‑149 (MQ=255)
                                                                                  cgcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAg                                                               <  1:4640756/149‑1 (MQ=255)
                                                                                  cgcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAg                                                               <  1:4640741/149‑1 (MQ=255)
                                                                                   gcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAgg                                                              <  2:3701129/149‑1 (MQ=255)
                                                                                   gcgTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAgg                                                              <  1:2032377/149‑1 (MQ=255)
                                                                                      ttAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGtt                                                           <  1:3241860/149‑1 (MQ=255)
                                                                                        aaCGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTg                                                         >  1:948830/1‑149 (MQ=255)
                                                                                         aCGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGt                                                        >  2:4634312/1‑149 (MQ=255)
                                                                                           gCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCg                                                      >  2:2768569/1‑149 (MQ=255)
                                                                                            cAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGa                                                     <  2:4690583/149‑1 (MQ=255)
                                                                                             aaCTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAg                                                    >  1:1898142/1‑149 (MQ=255)
                                                                                                  cgtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGaa                                               <  2:503598/149‑1 (MQ=255)
                                                                                                  cgtgcgCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGaa                                               <  2:2711157/149‑1 (MQ=255)
                                                                                                     gcgctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGa                                            <  2:2408501/149‑1 (MQ=255)
                                                                                                        ctgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAAt                                         <  2:3244386/149‑1 (MQ=255)
                                                                                                         tgcgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATc                                        <  1:4575409/149‑1 (MQ=255)
                                                                                                           cgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATTACAAAACGTTTTGAACAGCAGGTAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATccc                                      <  1:2438723/149‑1 (MQ=255)
                                                                                                           cgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATccc                                      <  1:235527/149‑1 (MQ=255)
                                                                                                           cgctATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATccc                                      <  2:1176545/149‑1 (MQ=255)
                                                                                                              tATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCga                                   >  1:880112/1‑149 (MQ=255)
                                                                                                              tATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCga                                   >  1:1130521/1‑149 (MQ=255)
                                                                                                              tATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCga                                   >  1:1659374/1‑149 (MQ=255)
                                                                                                               aTTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaa                                  <  2:4478528/149‑1 (MQ=255)
                                                                                                                ttGTAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCgaag                                 >  2:3092786/1‑149 (MQ=255)
                                                                                                                   tAAAGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAAc                              <  2:1761176/149‑1 (MQ=255)
                                                                                                                    aaaGAGATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACt                             <  2:568643/149‑1 (MQ=255)
                                                                                                                    aaaGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACt                             <  2:31410/149‑1 (MQ=255)
                                                                                                                    aaaGAGAGCCCTGCCCGGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACt                             <  2:3296549/149‑1 (MQ=255)
                                                                                                                      agagATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGc                           >  1:3715205/1‑149 (MQ=255)
                                                                                                                      agagATCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGc                           >  1:3427177/1‑149 (MQ=255)
                                                                                                                         gaTCCCTGCCCTGGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgc                        <  2:1339026/149‑1 (MQ=255)
                                                                                                                            ccctgccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgctgc                     <  1:2994438/149‑1 (MQ=255)
                                                                                                                            ccctgccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCgctgc                     >  1:4743860/1‑149 (MQ=255)
                                                                                                                               tgccctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCtgctgc                  <  1:2418949/149‑1 (MQ=255)
                                                                                                                                   ctgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGaa              >  1:2915095/1‑149 (MQ=255)
                                                                                                                                    tgGATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGaaa             <  1:365079/149‑1 (MQ=255)
                                                                                                                                      gATCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAaga           >  1:4713223/1‑149 (MQ=255)
                                                                                                                                        tCCGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGt         >  1:1805266/1‑149 (MQ=255)
                                                                                                                                         ccGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGCACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtc        <  1:1375573/149‑1 (MQ=255)
                                                                                                                                          cGCAACTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGtct       <  1:3953367/149‑1 (MQ=255)
                                                                                                                                             aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGt    >  1:3302161/1‑149 (MQ=255)
                                                                                                                                             aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGt    <  1:2481963/149‑1 (MQ=255)
                                                                                                                                             aaCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGt    >  1:174237/1‑149 (MQ=255)
                                                                                                                                              aCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt   >  2:18485/1‑149 (MQ=255)
                                                                                                                                              aCTGTTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGtt   <  2:868588/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TTTGCTTTACATCTCCCGTAAACACTTTTCTGCGATACAATGCCTTTACGTTATGTAACGGAGAGTTCGGCATGTCACACCCCGCGTTAACGCAACTGCGTGCGCTGCGCTATTGTAAAGAGATCCCTGCCCTGGATCCGCAACTGCTCGACTGGCTGTTGCTGGAGGATTCCATGACAAAACGTTTTGAACAGCAGGGAAAAACGGTAAGCGTGACGATGATCCGCGAAGGGTTTGTCGAGCAGAATGAAATCCCCGAAGAACTGCCGCTGCTGCCGAAAGAGTCTCGTTA  >  CP009273/4242363‑4242654

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: ‑