Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,458,948 G→T 88.5% D503Y (GAT→TAT)  mgtA → magnesium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,458,9480GT88.5% 673.3 / 84.3 265D503Y (GAT→TAT) mgtAmagnesium transporter
Reads supporting (aligned to +/- strand):  ref base G (16/15);  new base T (108/126);  total (124/141)
Fisher's exact test for biased strand distribution p-value = 5.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 268 total aligned reads displayed.

TGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAAC  >  CP009273/4458802‑4459093
                                                                                                                                                  |                                                                                                                                                 
 ggtAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTtata                                                                                                                                                <  1:4082164/149‑1 (MQ=255)
  gtAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATAc                                                                                                                                               >  2:925142/1‑149 (MQ=255)
     gTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCt                                                                                                                                            <  1:717901/149‑1 (MQ=255)
     gTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCt                                                                                                                                            <  1:1827478/149‑1 (MQ=255)
       ggCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGa                                                                                                                                          <  2:447235/149‑1 (MQ=255)
       ggCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGa                                                                                                                                          >  1:4695239/1‑149 (MQ=255)
         cAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAAt                                                                                                                                        <  1:4645319/149‑1 (MQ=255)
          aGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATc                                                                                                                                       >  2:3291850/1‑149 (MQ=255)
          aGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATc                                                                                                                                       >  1:2937027/1‑149 (MQ=255)
           gAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCg                                                                                                                                      <  1:3565994/149‑1 (MQ=255)
           gAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCg                                                                                                                                      <  1:2678505/149‑1 (MQ=255)
             aaaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTc                                                                                                                                    <  1:4589283/149‑1 (MQ=255)
              aaTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCa                                                                                                                                   >  1:1567283/1‑149 (MQ=255)
               aTACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAg                                                                                                                                  <  1:1806801/149‑1 (MQ=255)
                tACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgg                                                                                                                                 <  1:2784941/149‑1 (MQ=255)
                  ccGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAgggg                                                                                                                               <  1:4684903/149‑1 (MQ=255)
                     aGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTg                                                                                                                            <  1:1861971/149‑1 (MQ=255)
                      gCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTgc                                                                                                                           <  2:2650075/149‑1 (MQ=255)
                            tCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTgg                                                                                                                     >  2:1472044/1‑149 (MQ=255)
                            tCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTgg                                                                                                                     >  1:552981/1‑149 (MQ=255)
                                cTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGc                                                                                                                 <  1:467489/149‑1 (MQ=255)
                                    tttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                             <  2:4347812/149‑1 (MQ=255)
                                    tttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                             >  2:2365271/1‑149 (MQ=255)
                                    tttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                             >  1:366515/1‑149 (MQ=255)
                                    tttGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTg                                                                                                             <  2:2627/149‑1 (MQ=255)
                                    attGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCATGTGAGTCACAATGGCGATATTGTGCCTCTCGATGAAATCATTCTGAGTAAGATTAAGCGTGTTACTTATAAGCTTAATCGTCATGGTCTGCGCGTGGTTGCGGTg                                                                                                             <  2:4651486/148‑1 (MQ=255)
                                        cAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGcga                                                                                                         >  2:647434/1‑149 (MQ=255)
                                             gTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGaaa                                                                                                    >  1:3727758/1‑149 (MQ=255)
                                               gCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATa                                                                                                  <  2:1685267/149‑1 (MQ=255)
                                                 aTTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAcc                                                                                                >  1:1585095/1‑149 (MQ=255)
                                                 aTTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAcc                                                                                                >  1:1584424/1‑149 (MQ=255)
                                                     cAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGcc                                                                                            <  2:2201247/149‑1 (MQ=255)
                                                            tCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGt                                                                                     <  1:4122502/149‑1 (MQ=255)
                                                                  aTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGga                                                                               <  1:622766/149‑1 (MQ=255)
                                                                  aTGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGga                                                                               >  2:13419/1‑149 (MQ=255)
                                                                     tgtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAt                                                                            <  1:1043140/149‑1 (MQ=255)
                                                                     tgtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAt                                                                            <  1:2619128/149‑1 (MQ=255)
                                                                     tgtggTCGCAGGTGCGTCACAATGGCGAGAGTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACCTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAt                                                                            <  2:3237119/149‑1 (MQ=255)
                                                                      gtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAtt                                                                           <  1:1397711/149‑1 (MQ=255)
                                                                      gtgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAtt                                                                           <  2:2900972/149‑1 (MQ=255)
                                                                       tgtTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTa                                                                          <  2:3214394/149‑1 (MQ=255)
                                                                          tCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCa                                                                       <  2:4734015/149‑1 (MQ=255)
                                                                          tCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCa                                                                       <  1:798702/149‑1 (MQ=255)
                                                                           cGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAg                                                                      <  1:319354/149‑1 (MQ=255)
                                                                            gCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgc                                                                     <  1:222671/149‑1 (MQ=255)
                                                                             cAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcg                                                                    >  1:16748/1‑149 (MQ=255)
                                                                              aGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAgcgg                                                                   >  2:3455537/1‑149 (MQ=255)
                                                                                   cGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGa                                                              <  2:1791461/149‑1 (MQ=255)
                                                                                     tCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATg                                                            <  2:1437593/149‑1 (MQ=255)
                                                                                       acaATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGaa                                                          >  1:2928491/1‑149 (MQ=255)
                                                                                       acaATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGaa                                                          <  2:1629126/149‑1 (MQ=255)
                                                                                        caATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAAt                                                         <  1:1871718/149‑1 (MQ=255)
                                                                                        caATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAAt                                                         <  1:2448334/149‑1 (MQ=255)
                                                                                         aaTGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATAACTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATc                                                        >  1:4332918/1‑149 (MQ=255)
                                                                                          aTGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATcc                                                       <  1:1487743/149‑1 (MQ=255)
                                                                                          aTGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATcc                                                       >  2:866748/1‑149 (MQ=255)
                                                                                             gCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGAc                                                    <  2:4173292/149‑1 (MQ=255)
                                                                                              cGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGAcc                                                   <  2:1957822/149‑1 (MQ=255)
                                                                                                agaTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTg                                                 >  2:3102436/1‑149 (MQ=255)
                                                                                                  aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGAt                                               <  1:1885726/149‑1 (MQ=255)
                                                                                                  aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGAt                                               <  1:2165777/149‑1 (MQ=255)
                                                                                                  aTTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGAt                                               <  2:3515968/149‑1 (MQ=255)
                                                                                                   ttGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCAGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATc                                              >  1:3010179/1‑149 (MQ=255)
                                                                                                     gtgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCt                                            <  1:2156867/149‑1 (MQ=255)
                                                                                                      tgCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTc                                           <  2:3820414/149‑1 (MQ=255)
                                                                                                       gCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGGTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCg                                          >  1:3362629/1‑149 (MQ=255)
                                                                                                         cGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGaa                                        <  1:443055/149‑1 (MQ=255)
                                                                                                          gCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAg                                       >  2:4518894/1‑149 (MQ=255)
                                                                                                           cTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAgg                                      <  2:2969320/149‑1 (MQ=255)
                                                                                                               aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatat                                  >  2:1652672/1‑149 (MQ=255)
                                                                                                               aTGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGGGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatat                                  >  2:4039443/1‑149 (MQ=255)
                                                                                                                tGACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatata                                 >  1:4479710/1‑149 (MQ=255)
                                                                                                                 gACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatatat                                <  2:1152163/149‑1 (MQ=255)
                                                                                                                 gACATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGatatat                                >  2:100960/1‑149 (MQ=255)
                                                                                                                  aCATCATGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATAtt                               <  1:228214/149‑1 (MQ=255)
                                                                                                                  aCATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATAtt                               >  2:2235936/1‑149 (MQ=255)
                                                                                                                     tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCt                            >  1:685652/1‑149 (MQ=255)
                                                                                                                     tcatGCTGCGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTCGTTGCGGAGGCGACGAAATACCAGCCAGCGCGAGAAGGAGATAACCAGCGGGCGGATGACTCCGAACTGATCCTCGAAGCCTCTATTGCt                            >  2:4638333/1‑149 (MQ=255)
                                                                                                                          ctgcGTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCt                       <  1:803003/149‑1 (MQ=255)
                                                                                                                             cGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGa                    <  2:4599314/149‑1 (MQ=255)
                                                                                                                              gTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   >  2:970079/1‑149 (MQ=255)
                                                                                                                              gTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   >  2:1876483/1‑149 (MQ=255)
                                                                                                                              gTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   <  1:756846/149‑1 (MQ=255)
                                                                                                                              gTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   >  1:2000606/1‑149 (MQ=255)
                                                                                                                              gTAAGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   >  1:3107632/1‑149 (MQ=255)
                                                                                                                              gTAAGATTAAGCGGGTTACTGAGACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGAt                   <  1:3307743/149‑1 (MQ=255)
                                                                                                                                aaGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                 <  2:1417900/149‑1 (MQ=255)
                                                                                                                                aaGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGGTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATcc                 <  2:2699608/149‑1 (MQ=255)
                                                                                                                                 aGATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCg                >  1:3797101/1‑149 (MQ=255)
                                                                                                                                  gATTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGc               >  2:2795577/1‑149 (MQ=255)
                                                                                                                                   aTTAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGcc              >  1:1236151/1‑149 (MQ=255)
                                                                                                                                    ttAAGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCa             >  1:3658784/1‑149 (MQ=255)
                                                                                                                                      aaGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCaaa           >  1:2340001/1‑149 (MQ=255)
                                                                                                                                       aGCGGGTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCaaaa          >  2:4084603/1‑149 (MQ=255)
                                                                                                                                           ggTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAgaga      <  2:4558512/149‑1 (MQ=255)
                                                                                                                                           ggTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGAGTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAgaga      <  1:948842/149‑1 (MQ=255)
                                                                                                                                            gTTACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGAc     <  1:2907447/149‑1 (MQ=255)
                                                                                                                                             ttACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACa    >  2:763598/1‑149 (MQ=255)
                                                                                                                                              tACTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACaa   >  1:4746507/1‑149 (MQ=255)
                                                                                                                                               aCTTATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAAc  <  1:1143469/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
TGGTAGTGGCAGAAAATACCGAGCACCATCAGCTGGTTTGCAAAGGTGCATTGCAGGAAATCCTCAATGTGTGTTCGCAGGTGCGTCACAATGGCGAGATTGTGCCGCTCGATGACATCATGCTGCGTAAGATTAAGCGGGTTACTGATACGCTGAATCGTCAGGGGCTGCGCGTGGTTGCGGTGGCGACGAAATACCTGCCAGCGCGTGAAGGAGATTACCAGCGGGCGGATGAATCCGACCTGATCCTCGAAGGATATATTGCTTTTCTTGATCCGCCAAAAGAGACAAC  >  CP009273/4458802‑4459093

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