Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,505,141 G→T 100% R452S (CGT→AGT)  fecA ← ferric citrate outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,505,1410GT95.7% 627.5 / 17.9 211R452S (CGT→AGT) fecAferric citrate outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (0/9);  new base T (88/114);  total (88/123)
Fisher's exact test for biased strand distribution p-value = 1.13e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-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 211 total aligned reads displayed.

GCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTG  >  CP009273/4504997‑4505285
                                                                                                                                                |                                                                                                                                                
gCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCAt                                                                                                                                               <  1:3445618/148‑1 (MQ=255)
gCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCAt                                                                                                                                               <  1:3037241/148‑1 (MQ=255)
gCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCAt                                                                                                                                               <  1:3037224/148‑1 (MQ=255)
gCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCAt                                                                                                                                               >  1:2735337/1‑148 (MQ=255)
 cATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCAtt                                                                                                                                              >  1:1218040/1‑148 (MQ=255)
 cATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCAtt                                                                                                                                              >  2:663898/1‑148 (MQ=255)
 cATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCAtt                                                                                                                                              <  2:2361600/148‑1 (MQ=255)
  aTACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCAttt                                                                                                                                             <  2:2456789/148‑1 (MQ=255)
  aTACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCAttt                                                                                                                                             <  2:671848/148‑1 (MQ=255)
    aCCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCa                                                                                                                                           <  2:2955850/148‑1 (MQ=255)
         gCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGt                                                                                                                                      <  1:3162405/148‑1 (MQ=255)
          cGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGtt                                                                                                                                     >  1:480664/1‑148 (MQ=255)
           gTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTg                                                                                                                                    <  2:3279915/148‑1 (MQ=255)
            tGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTACGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGa                                                                                                                                   >  2:786973/1‑148 (MQ=255)
             gATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGat                                                                                                                                  <  2:1138095/148‑1 (MQ=255)
              aTGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCACGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGatt                                                                                                                                 >  1:1004542/1‑148 (MQ=255)
                 gTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGattcat                                                                                                                              >  2:647072/1‑148 (MQ=255)
                 gTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGattcat                                                                                                                              <  1:1883770/148‑1 (MQ=255)
                  tCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGattcatt                                                                                                                             <  1:570793/148‑1 (MQ=255)
                   ccAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGattcattc                                                                                                                            <  2:2290496/148‑1 (MQ=255)
                   ccAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGattcattc                                                                                                                            <  2:1534127/148‑1 (MQ=255)
                      gTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCaaa                                                                                                                         <  1:2393426/148‑1 (MQ=255)
                         gCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAg                                                                                                                      <  2:2791164/148‑1 (MQ=255)
                             aTGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGAt                                                                                                                  <  2:869283/148‑1 (MQ=255)
                             aTGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGAt                                                                                                                  <  2:2203320/148‑1 (MQ=255)
                             aTGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGAt                                                                                                                  <  2:3350219/148‑1 (MQ=255)
                              tGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATa                                                                                                                 >  2:873014/1‑148 (MQ=255)
                               gTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAg                                                                                                                >  1:1306662/1‑148 (MQ=255)
                               gTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTAAAATATCGATAg                                                                                                                >  1:653829/1‑148 (MQ=255)
                               gTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCGTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGGAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAg                                                                                                                >  1:636815/1‑148 (MQ=255)
                                 cGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcc                                                                                                              >  1:3321444/1‑148 (MQ=255)
                                  gATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgccc                                                                                                             <  2:3456437/148‑1 (MQ=255)
                                  gATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgccc                                                                                                             >  1:2662316/1‑148 (MQ=255)
                                   aTTTTGTCATCCAGATACCATGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgccca                                                                                                            >  1:804509/1‑148 (MQ=255)
                                   aTTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgccca                                                                                                            >  1:2558251/1‑148 (MQ=255)
                                       tGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgcc                                                                                                        <  2:3295462/148‑1 (MQ=255)
                                          cATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgcccac                                                                                                     >  2:271723/1‑148 (MQ=255)
                                          cATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgcccac                                                                                                     >  1:1586075/1‑148 (MQ=255)
                                          cATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgcccac                                                                                                     >  1:1913041/1‑148 (MQ=255)
                                          cATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAgcccacgcccac                                                                                                     >  2:271737/1‑148 (MQ=255)
                                           aTCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACt                                                                                                    >  1:3302986/1‑148 (MQ=255)
                                            tCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACtt                                                                                                   >  2:740013/1‑148 (MQ=255)
                                            tCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACtt                                                                                                   <  1:3341011/148‑1 (MQ=255)
                                            tCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACtt                                                                                                   <  1:1603744/148‑1 (MQ=255)
                                              cAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCTCGCCCACTTCg                                                                                                 <  2:2787801/148‑1 (MQ=255)
                                                    ccAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGcgcg                                                                                           >  1:3153962/1‑148 (MQ=255)
                                                     cAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCgg                                                                                          >  2:1233945/1‑148 (MQ=255)
                                                     cAGGCGTGCGCCTCGGTGCCGGAACGCGTATCACGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCgg                                                                                          <  2:1004542/148‑1 (MQ=255)
                                                      aGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGa                                                                                         <  2:219291/148‑1 (MQ=255)
                                                      aGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGa                                                                                         >  1:2999099/1‑148 (MQ=255)
                                                       ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGaa                                                                                        >  1:343723/1‑148 (MQ=255)
                                                       ggCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGaa                                                                                        >  1:1810458/1‑148 (MQ=255)
                                                        gCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAg                                                                                       >  2:2607057/1‑148 (MQ=255)
                                                        gCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAg                                                                                       <  2:41665/148‑1 (MQ=255)
                                                        gCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAg                                                                                       >  2:2120160/1‑148 (MQ=255)
                                                        gCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAg                                                                                       >  2:1774694/1‑148 (MQ=255)
                                                                cTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATc                                                                               >  1:711907/1‑148 (MQ=255)
                                                                cTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATc                                                                               >  1:2088041/1‑148 (MQ=255)
                                                                           aaCGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTg                                                                    <  1:1403257/148‑1 (MQ=255)
                                                                            aCGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTgg                                                                   >  1:1067067/1‑148 (MQ=255)
                                                                             cgcgTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGc                                                                  <  2:1793784/148‑1 (MQ=255)
                                                                                 tATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTa                                                              <  2:2902344/148‑1 (MQ=255)
                                                                                    cgcgGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg                                                           >  2:1002829/1‑148 (MQ=255)
                                                                                    cgcgGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCg                                                           <  2:1598085/148‑1 (MQ=255)
                                                                                        gTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGc                                                       >  1:3114322/1‑148 (MQ=255)
                                                                                             aaGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAAt                                                  >  2:3499695/1‑148 (MQ=255)
                                                                                              aGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATa                                                 <  1:352571/148‑1 (MQ=255)
                                                                                              aGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATa                                                 >  1:2347704/1‑148 (MQ=255)
                                                                                               gggCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATAc                                                >  2:1934476/1‑148 (MQ=255)
                                                                                               gggCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGATGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCGGTAGCGTGGCTCAATAc                                                >  2:136514/1‑148 (MQ=255)
                                                                                                      aGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCAc                                         <  2:2079219/148‑1 (MQ=255)
                                                                                                      aGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCAc                                         >  2:2808980/1‑148 (MQ=255)
                                                                                                       gCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCAcc                                        <  1:1690608/148‑1 (MQ=255)
                                                                                                        ccGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCAccc                                       >  2:1012384/1‑148 (MQ=255)
                                                                                                         cGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCa                                      <  1:271726/148‑1 (MQ=255)
                                                                                                         cGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCa                                      <  1:2751029/148‑1 (MQ=255)
                                                                                                             cGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCagtag                                  <  2:804509/148‑1 (MQ=255)
                                                                                                               gCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGtt                                <  2:1826454/148‑1 (MQ=255)
                                                                                                               gCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGtt                                <  1:3233496/148‑1 (MQ=255)
                                                                                                               gCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGtt                                <  2:1860948/148‑1 (MQ=255)
                                                                                                                cAACTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTa                               <  1:2195990/148‑1 (MQ=255)
                                                                                                                  aCTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTAcg                             <  2:1987067/148‑1 (MQ=255)
                                                                                                                  aCTGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTAcg                             <  2:1685418/148‑1 (MQ=255)
                                                                                                                    tGCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTAcgcg                           >  1:3051596/1‑148 (MQ=255)
                                                                                                                     gCCCGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCTTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACgcgg                          <  2:3513389/148‑1 (MQ=255)
                                                                                                                        cGCTGCTGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCga                       <  1:1266779/148‑1 (MQ=255)
                                                                                                                         gctgctGGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCgag                      >  2:3053992/1‑148 (MQ=255)
                                                                                                                          ctgctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGTGGCgaga                     <  1:1936714/148‑1 (MQ=255)
                                                                                                                             ctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAggg                  <  2:2328062/148‑1 (MQ=255)
                                                                                                                             ctgGTGGCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAggg                  <  2:1045628/148‑1 (MQ=255)
                                                                                                                               ggtggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTg                <  2:3302986/148‑1 (MQ=255)
                                                                                                                               ggtggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTg                <  1:2017257/148‑1 (MQ=255)
                                                                                                                                gtggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGa               <  2:1135166/148‑1 (MQ=255)
                                                                                                                                 tggCGGTGTAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGAt              <  2:1913079/148‑1 (MQ=255)
                                                                                                                                       tgtAGTAACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGttt        >  1:631829/1‑148 (MQ=255)
                                                                                                                                        gtagtaACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTg       <  1:1893517/148‑1 (MQ=255)
                                                                                                                                        gtagtaACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTg       <  1:2460289/148‑1 (MQ=255)
                                                                                                                                         tagtaACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGc      <  1:113226/148‑1 (MQ=255)
                                                                                                                                           gtaACTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCt    <  2:2285772/148‑1 (MQ=255)
                                                                                                                                             aaCTCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTg  <  1:2987892/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
GCATACCCGGCGTGATGGTCCAGTTGCCGATGTCGATTTTGTCATCCAGATACCAGGCGTGCGCCTCGGTGCCGGAACGCGTATCGCGGTCGTAAGGGCTTGAGCCGGACGGCAACTGCCCGCTGCTGGTGGCGGTGTAGTAACGCATTTCATGCGTTGATTCATTCAAATAGCGATAGCCCACGCCCACTTCGTGCGCGGAAGGGCCGATCATAAAGATCTGGCTGTAGCGTGGCTCAATACCGCGCACCCAGTAGTTACGCGGCGAGAGGGTGATGCGTTTGCCTTG  >  CP009273/4504997‑4505285

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: