Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,403,195 G→T 8.8% S546* (TCG→TAG)  yfbS ← putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,403,1950GT 8.8% 1202.6 / 66.0 317S546* (TCG→TAG) yfbSputative transporter
Reads supporting (aligned to +/- strand):  ref base G (118/171);  new base T (15/13);  total (133/184)
Fisher's exact test for biased strand distribution p-value = 2.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 320 total aligned reads displayed.

GATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCAT  >  CP009273/2403047‑2403343
                                                                                                                                                    |                                                                                                                                                    
 aTAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGa                                                                                                                                                     >  1:2604189/1‑149 (MQ=255)
  tAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGac                                                                                                                                                    >  1:4253510/1‑149 (MQ=255)
  tAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGac                                                                                                                                                    >  1:3395774/1‑149 (MQ=255)
   aaTGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGaca                                                                                                                                                   >  1:233805/1‑149 (MQ=255)
        tGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACAccccc                                                                                                                                              >  1:447405/1‑149 (MQ=255)
        tGAACGGTACCCCCAACTTAACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACAccccc                                                                                                                                              <  1:1078787/149‑1 (MQ=255)
           aCGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATc                                                                                                                                           <  2:555741/149‑1 (MQ=255)
              gTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCTACACCCCCATCGtt                                                                                                                                        <  2:4765864/149‑1 (MQ=255)
              gTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGtt                                                                                                                                        <  1:422670/149‑1 (MQ=255)
               tACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGttt                                                                                                                                       >  2:3997962/1‑149 (MQ=255)
                aCCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGTTAAGGCTACACCCCCATCGtttt                                                                                                                                      <  1:2836407/149‑1 (MQ=255)
                aCCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGtttt                                                                                                                                      >  1:4074611/1‑149 (MQ=255)
                  ccccAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCTACACCCCCATCGTTTTgg                                                                                                                                    <  1:89375/149‑1 (MQ=255)
                   cccAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGgc                                                                                                                                   >  2:780060/1‑149 (MQ=255)
                   cccAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGgc                                                                                                                                   <  2:2353693/149‑1 (MQ=255)
                      aaCTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGgcagc                                                                                                                                >  2:1936063/1‑149 (MQ=255)
                       aCTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGgcagca                                                                                                                               <  2:2988897/149‑1 (MQ=255)
                       aCTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGgcagca                                                                                                                               <  2:3133410/149‑1 (MQ=255)
                         ttCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCagcagc                                                                                                                             >  2:4092198/1‑149 (MQ=255)
                          tCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGcagcc                                                                                                                            <  1:1590990/149‑1 (MQ=255)
                           cacaAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGcagcca                                                                                                                           <  2:1982339/149‑1 (MQ=255)
                            acaAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCTACACCCCCATCGTTTTGGCAGcagccag                                                                                                                          >  2:1476108/1‑149 (MQ=255)
                                 gTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCaa                                                                                                                     <  1:98757/149‑1 (MQ=255)
                                 gTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCaa                                                                                                                     <  2:4004322/149‑1 (MQ=255)
                                  tCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAAt                                                                                                                    <  2:1137735/149‑1 (MQ=255)
                                   cACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATc                                                                                                                   <  2:4150871/149‑1 (MQ=255)
                                   cACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATc                                                                                                                   <  2:1571631/149‑1 (MQ=255)
                                     cTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCgg                                                                                                                 >  2:4634487/1‑149 (MQ=255)
                                       gAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCTACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAg                                                                                                               >  2:5057174/1‑149 (MQ=255)
                                        aaGCTGTAATTTTCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGc                                                                                                              >  1:2192681/1‑149 (MQ=255)
                                            tGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATc                                                                                                          >  1:2540838/1‑149 (MQ=255)
                                            tGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATc                                                                                                          >  1:3954024/1‑149 (MQ=255)
                                            tGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATc                                                                                                          <  1:2952777/149‑1 (MQ=255)
                                              tAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCaa                                                                                                        >  1:1388646/1‑149 (MQ=255)
                                              tAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCaa                                                                                                        <  2:218290/149‑1 (MQ=255)
                                                aTTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCTACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAaca                                                                                                      <  2:4110265/149‑1 (MQ=255)
                                                  ttCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACAcc                                                                                                    <  2:1920704/149‑1 (MQ=255)
                                                         aCCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTcg                                                                                             >  2:1246301/1‑149 (MQ=255)
                                                            tAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCgg                                                                                          <  1:3220444/149‑1 (MQ=255)
                                                              aaaCCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTa                                                                                        <  2:1378992/149‑1 (MQ=255)
                                                              aaaCCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTa                                                                                        <  1:1565443/149‑1 (MQ=255)
                                                              aaaCCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTa                                                                                        <  1:2984535/149‑1 (MQ=255)
                                                               aaCCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATACAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTAt                                                                                       <  2:3532689/149‑1 (MQ=255)
                                                                aCCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTAtt                                                                                      <  2:2623934/149‑1 (MQ=255)
                                                                   aGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTaga                                                                                   <  1:2617210/149‑1 (MQ=255)
                                                                      gtgtTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGTGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGAt                                                                                <  1:4944806/149‑1 (MQ=255)
                                                                       tgtTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATa                                                                               >  1:631516/1‑149 (MQ=255)
                                                                       tgtTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATa                                                                               <  2:4977537/149‑1 (MQ=255)
                                                                         ttAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATaaa                                                                             <  2:3173870/149‑1 (MQ=255)
                                                                            aCAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAg                                                                          >  2:2046796/1‑149 (MQ=255)
                                                                             cAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGc                                                                         <  1:3722804/149‑1 (MQ=255)
                                                                              aGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTa                                                                                        <  2:2789808/133‑1 (MQ=255)
                                                                                   aagaaACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATa                                                                   <  2:2550101/149‑1 (MQ=255)
                                                                                   aagaaACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATa                                                                   <  2:3739881/149‑1 (MQ=255)
                                                                                   aagaaACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATa                                                                   <  2:3870394/149‑1 (MQ=255)
                                                                                    agaaACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATaa                                                                  <  1:2401004/149‑1 (MQ=255)
                                                                                       aaCCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAAccg                                                               >  1:694778/1‑149 (MQ=255)
                                                                                       aaCCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAAccg                                                               >  2:2449145/1‑149 (MQ=255)
                                                                                       aaCCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAAccg                                                               <  1:3156960/149‑1 (MQ=255)
                                                                                        aCCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTCTTAGAGATAAATAGCCCAATAAccgg                                                              >  1:2259829/1‑148 (MQ=255)
                                                                                        aCCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAAAAAccgc                                                              >  2:4848822/1‑149 (MQ=255)
                                                                                          cGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAAccgccg                                                            <  2:1266293/149‑1 (MQ=255)
                                                                                                  tAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCTACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCgacaagaca                                                    >  1:3750795/1‑149 (MQ=255)
                                                                                                     aGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACACGacaaac                                                 >  1:448999/1‑149 (MQ=255)
                                                                                                     aGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGacaaac                                                 <  1:4890859/149‑1 (MQ=255)
                                                                                                      ggcggcGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGacaaaca                                                <  1:5128883/149‑1 (MQ=255)
                                                                                                         ggcggATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACaaacaaac                                             <  1:3791430/149‑1 (MQ=255)
                                                                                                            ggATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACAcc                                          >  1:4679453/1‑149 (MQ=255)
                                                                                                             gATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACAccc                                         >  1:4446293/1‑149 (MQ=255)
                                                                                                                gctgcCATCGCCACGACCATCGCGAATGGATAAGGCTACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAg                                      <  2:336026/149‑1 (MQ=255)
                                                                                                                  tgcCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGca                                    <  2:3852657/149‑1 (MQ=255)
                                                                                                                      aTCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatcat                                <  2:42785/149‑1 (MQ=255)
                                                                                                                      aTCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGcatcat                                <  1:686618/149‑1 (MQ=255)
                                                                                                                       tCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                               <  2:4726112/149‑1 (MQ=255)
                                                                                                                       tCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                               >  1:1271437/1‑149 (MQ=255)
                                                                                                                       tCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCata                               <  2:1493280/149‑1 (MQ=255)
                                                                                                                        cGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCatat                              >  1:1099642/1‑149 (MQ=255)
                                                                                                                        cGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCatat                              <  2:1671255/149‑1 (MQ=255)
                                                                                                                          ccACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATAtgt                            <  2:683476/149‑1 (MQ=255)
                                                                                                                          ccACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATAtgt                            <  1:1717549/149‑1 (MQ=255)
                                                                                                                           cACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACACACAAACACCCCAGCATCATAtgtg                           >  1:4522820/1‑149 (MQ=255)
                                                                                                                            aCGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTgg                          >  1:1152749/1‑149 (MQ=255)
                                                                                                                             cGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGc                         <  1:271853/149‑1 (MQ=255)
                                                                                                                              gACCATCGCGAATGGATAAGGCTACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGcc                        >  1:1939815/1‑149 (MQ=255)
                                                                                                                               aCCATCGCGAATGGATAAGGCTACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGccc                       >  2:4346919/1‑149 (MQ=255)
                                                                                                                               aCCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGccc                       >  1:2190099/1‑149 (MQ=255)
                                                                                                                               aCCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGccc                       <  1:574233/149‑1 (MQ=255)
                                                                                                                               aCCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGccc                       <  1:4841031/149‑1 (MQ=255)
                                                                                                                                   tCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTaa                   >  2:1343508/1‑149 (MQ=255)
                                                                                                                                   tCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTaa                   >  2:4576060/1‑149 (MQ=255)
                                                                                                                                     gcgAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAacc                 <  1:3023249/149‑1 (MQ=255)
                                                                                                                                       gAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAaccac               <  1:3501960/149‑1 (MQ=255)
                                                                                                                                       gAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAaccac               <  2:3909974/149‑1 (MQ=255)
                                                                                                                                        aaTGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAaccacc              >  2:3246466/1‑149 (MQ=255)
                                                                                                                                                aGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGt      >  1:203567/1‑149 (MQ=255)
                                                                                                                                                aGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGt      <  2:3623470/149‑1 (MQ=255)
                                                                                                                                                 ggCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTc     <  1:3903225/149‑1 (MQ=255)
                                                                                                                                                  gCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTcc    <  1:4835098/149‑1 (MQ=255)
                                                                                                                                                   cGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCa   <  2:2241798/149‑1 (MQ=255)
                                                                                                                                                    gACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCAt  >  1:565468/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GATAATGGTGAACGGTACCCCCAACTTCACAAAGTCACTGAAGCTGTAATTTCCCGGACCTAAAACCAGTGTGTTAACAGGTGAAGAAACCGGTGTCATAAAGGCGGCGGATGCTGCCATCGCCACGACCATCGCGAATGGATAAGGCGACACCCCCATCGTTTTGGCAGCAGCCAGCGCAATCGGAGCCATCAACACCGCCGTCGCGGTATTAGAGATAAATAGCCCAATAACCGCCGACAAGACAAACAAACACCCCAGCATCATATGTGGCCCGTAACCACCGCCAATGTCCAT  >  CP009273/2403047‑2403343

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: