Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 607,660 G→T 9.1% R97R (CGC→CGA fepA ← iron‑enterobactin outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273607,6600GT 9.1% 1080.3 / 63.5 286R97R (CGC→CGAfepAiron‑enterobactin outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (125/135);  new base T (12/14);  total (137/149)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

Only 100 of 302 total aligned reads displayed.

GAACTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCAG  >  CP009273/607512‑607808
                                                                                                                                                    |                                                                                                                                                    
gAACTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACcg                                                                                                                                                      >  1:1047836/1‑149 (MQ=255)
gAACTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACcg                                                                                                                                                      >  2:2317730/1‑149 (MQ=255)
 aaCTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACcgc                                                                                                                                                     >  2:3655017/1‑149 (MQ=255)
  aCTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACcgcg                                                                                                                                                    >  1:754682/1‑149 (MQ=255)
    ttCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGaa                                                                                                                                                  <  1:4042833/149‑1 (MQ=255)
     tCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAat                                                                                                                                                 <  1:263703/149‑1 (MQ=255)
     tCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAat                                                                                                                                                 >  2:2788836/1‑149 (MQ=255)
     tCAATACGTTCAATCATTTCAGGTGGAACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGACAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAat                                                                                                                                                 >  2:2823300/1‑149 (MQ=255)
      cAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAata                                                                                                                                                <  1:3310582/149‑1 (MQ=255)
       aaTACGTTCTTTCATTTCAGTTGGCACACAGGAAGTATCACCACGGGTATCTCGCTCGCCTAGCAAGCCCTGTCGCTCCGATTTACGGCTGCTTACCGGCTTGCCGTCAATCATAATCAGCGTGTTTTCCGGACCCATACCGCGAatat                                                                                                                                               <  1:1475507/149‑1 (MQ=255)
        agaCGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATc                                                                                                                                              <  2:5067654/147‑1 (MQ=255)
         tACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCa                                                                                                                                             >  1:2695101/1‑149 (MQ=255)
          aCGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCaa                                                                                                                                            >  2:1057925/1‑149 (MQ=255)
          aCGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCaa                                                                                                                                            <  1:5145519/149‑1 (MQ=255)
          aCGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCaa                                                                                                                                            <  2:4848730/149‑1 (MQ=255)
            gTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCATCGTGTTTTCCGGACCCATACCGCGAATATCAATc                                                                                                                                          >  2:3569425/1‑149 (MQ=255)
            gTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATc                                                                                                                                          >  2:1744527/1‑149 (MQ=255)
               cAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGt                                                                                                                                       >  2:4168862/1‑149 (MQ=255)
                aaTCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTc                                                                                                                                      <  2:3231465/149‑1 (MQ=255)
                aaTCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTc                                                                                                                                      <  2:3446367/149‑1 (MQ=255)
                  tCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCgg                                                                                                                                    <  2:2191720/149‑1 (MQ=255)
                    aTTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGtt                                                                                                                                  >  2:1524486/1‑149 (MQ=255)
                      ttCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTAt                                                                                                                                <  1:3123494/149‑1 (MQ=255)
                      ttCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTAt                                                                                                                                <  1:3529788/149‑1 (MQ=255)
                       tCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTAtt                                                                                                                               <  1:1357435/149‑1 (MQ=255)
                         aGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTcc                                                                                                                             <  2:488866/149‑1 (MQ=255)
                          ggTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTccc                                                                                                                            <  1:4112060/149‑1 (MQ=255)
                           gTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCa                                                                                                                           <  1:4046409/149‑1 (MQ=255)
                           gTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCa                                                                                                                           <  1:3663625/149‑1 (MQ=255)
                               cACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCt                                                                                                                       >  2:64450/1‑149 (MQ=255)
                                   cAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGATTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTAACACGCTGAcc                                                                                                                   >  1:2863229/1‑149 (MQ=255)
                                   cAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGAcc                                                                                                                   >  2:3326435/1‑149 (MQ=255)
                                      gAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACt                                                                                                                >  1:32108/1‑149 (MQ=255)
                                      gAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACt                                                                                                                >  2:462398/1‑149 (MQ=255)
                                        aGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACtgg                                                                                                              >  2:4080821/1‑149 (MQ=255)
                                        aGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACtgg                                                                                                              >  1:1585617/1‑149 (MQ=255)
                                          tATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACtggtg                                                                                                            <  1:5048261/149‑1 (MQ=255)
                                            tCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGa                                                                                                          <  1:4206779/149‑1 (MQ=255)
                                            tCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGa                                                                                                          <  2:4033839/149‑1 (MQ=255)
                                            tCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCGTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGGTATTCCCACGCTGACCACTGGGGGa                                                                                                          >  1:687115/1‑149 (MQ=255)
                                                 acGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTAc                                                                                                     <  1:3152633/149‑1 (MQ=255)
                                                 acGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTAc                                                                                                     <  1:1734453/149‑1 (MQ=255)
                                                    ggTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCgg                                                                                                  >  1:5103315/1‑149 (MQ=255)
                                                           cgcTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGtt                                                                                           <  1:3729486/149‑1 (MQ=255)
                                                             cTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTaa                                                                                         <  2:1656322/149‑1 (MQ=255)
                                                              tCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAAc                                                                                        >  2:1874642/1‑149 (MQ=255)
                                                                   acgccaGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTg                                                                                   >  2:1565255/1‑149 (MQ=255)
                                                                    cgccaGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTgg                                                                                  <  1:4423999/149‑1 (MQ=255)
                                                                     gccaGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGc                                                                                 >  1:211816/1‑149 (MQ=255)
                                                                       caGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAACTCCGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCAt                                                                               <  1:2269785/149‑1 (MQ=255)
                                                                        aGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATg                                                                              >  1:3075103/1‑149 (MQ=255)
                                                                        aGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATg                                                                              <  1:2927985/149‑1 (MQ=255)
                                                                          cccTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGt                                                                            >  1:3114924/1‑149 (MQ=255)
                                                                           ccTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTa                                                                           <  2:2438518/149‑1 (MQ=255)
                                                                           ccTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTa                                                                           <  1:1195318/149‑1 (MQ=255)
                                                                            cTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTAc                                                                          >  1:1440527/1‑149 (MQ=255)
                                                                               aCGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGga                                                                       >  1:4568812/1‑149 (MQ=255)
                                                                                cGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGgat                                                                      <  2:4739429/149‑1 (MQ=255)
                                                                                       aGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTc                                                               <  2:2259132/149‑1 (MQ=255)
                                                                                       aGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTc                                                               <  2:3078782/149‑1 (MQ=255)
                                                                                        gTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCg                                                              >  2:5146057/1‑149 (MQ=255)
                                                                                        gTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCg                                                              >  2:4090799/1‑149 (MQ=255)
                                                                                             ggCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGAcaca                                                         >  1:3010895/1‑149 (MQ=255)
                                                                                             ggCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGAcaca                                                         >  1:2573181/1‑149 (MQ=255)
                                                                                             ggCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGAcaca                                                         >  1:4133580/1‑149 (MQ=255)
                                                                                              gctgctTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACAt                                                        >  1:916202/1‑149 (MQ=255)
                                                                                              gctgctTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACAt                                                        <  2:1047836/149‑1 (MQ=255)
                                                                                              gctgctTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACAt                                                        <  2:2913355/149‑1 (MQ=255)
                                                                                                        ggCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAAc                                              >  1:4359286/1‑149 (MQ=255)
                                                                                                          cTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCg                                            <  1:2788896/149‑1 (MQ=255)
                                                                                                             gCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGt                                         <  2:2560244/149‑1 (MQ=255)
                                                                                                              ccGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGtt                                        <  2:2063402/149‑1 (MQ=255)
                                                                                                                   aatcaaAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGc                                   <  1:4582756/149‑1 (MQ=255)
                                                                                                                   aatcaaAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGc                                   <  1:1988905/149‑1 (MQ=255)
                                                                                                                   aatcaaAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGc                                   <  1:4797208/149‑1 (MQ=255)
                                                                                                                    atcaaAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCg                                  >  1:2361891/1‑149 (MQ=255)
                                                                                                                     tcaaAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCgg                                 >  2:3965391/1‑149 (MQ=255)
                                                                                                                     tcaaAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCgg                                 >  1:1599511/1‑149 (MQ=255)
                                                                                                                       aaaaTCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGAt                               <  1:4677177/149‑1 (MQ=255)
                                                                                                                       aaaaTCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGAt                               >  2:2005486/1‑149 (MQ=255)
                                                                                                                       aaaaTCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGAt                               >  1:203888/1‑149 (MQ=255)
                                                                                                                        aaaTCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGAtt                              <  2:1841367/149‑1 (MQ=255)
                                                                                                                        aaaTCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGAtt                              <  2:4436366/149‑1 (MQ=255)
                                                                                                                          aTCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTc                            <  2:2485424/149‑1 (MQ=255)
                                                                                                                           tCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCa                           >  1:2224316/1‑149 (MQ=255)
                                                                                                                           tCAGCGTGTTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCa                           >  2:1272016/1‑149 (MQ=255)
                                                                                                                                 tgtTTTCCGGACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCg                     <  1:3851822/149‑1 (MQ=255)
                                                                                                                                 tgtTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCg                     >  1:4951450/1‑149 (MQ=255)
                                                                                                                                     ttCCGGACCCATACCGCGAATATCAATCTGTCGGTTAGTCCCACGCTGACCCCGGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGa                 <  2:1205332/149‑1 (MQ=255)
                                                                                                                                       ccGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATg               <  1:664888/149‑1 (MQ=255)
                                                                                                                                         ggACCCATACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGt             <  2:26884/149‑1 (MQ=255)
                                                                                                                                         ggACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGt             <  1:5077366/149‑1 (MQ=255)
                                                                                                                                         ggACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGt             >  1:3433480/1‑149 (MQ=255)
                                                                                                                                              cATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGaaa        <  1:1837477/149‑1 (MQ=255)
                                                                                                                                               aTACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAAc       >  2:3046046/1‑149 (MQ=255)
                                                                                                                                               aTACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAAc       >  2:1722502/1‑149 (MQ=255)
                                                                                                                                                tACCTCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACg      >  2:1085453/1‑149 (MQ=255)
                                                                                                                                                 aCCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGc     >  1:3897576/1‑149 (MQ=255)
                                                                                                                                                   cgcgAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCa   >  1:165652/1‑149 (MQ=255)
                                                                                                                                                   cgcgAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCa   <  2:1403697/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GAACTTCAATACGTTCAATCATTTCAGGTGGCACCCAGGAAGTATCACCACGGGTATCGCGCTCGCCACGCCAGCCCTGACGCACCGAGTTACGGCTGCTTACCGGCTTGCCGTCAATCAAAATCAGCGTGTTTTCCGGACCCATACCGCGAATATCAATCTGTCGGTTATTCCCACGCTGACCACTGGTGGAGTTACCGGTCAGGTTAACGCCTGGCATGGTACGGATGATCTTCGACACATCGCGGGCAACCGGGTTTTTGCGGATTTCATCTGCGGTGATGGTCGAAACGCCAG  >  CP009273/607512‑607808

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: