Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,676,886 C→A 5.8% intergenic (‑79/+104) dctA ← / ← yhjK C4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,676,8860CA 5.8% 1070.1 / 25.4 276intergenic (‑79/+104)dctA/yhjKC4‑dicarboxylic acid, orotate and citrate transporter/cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (127/133);  new base A (7/9);  total (134/142)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

Only 100 of 286 total aligned reads displayed.

CAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAG  >  CP009273/3676738‑3677034
                                                                                                                                                    |                                                                                                                                                    
cAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATc                                                                                                                                                      <  1:4790503/149‑1 (MQ=255)
cAAGGAGAAGACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATc                                                                                                                                                      <  2:4296680/149‑1 (MQ=255)
  aGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcac                                                                                                                                                    <  1:1178730/149‑1 (MQ=255)
  aGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcac                                                                                                                                                    <  2:3402220/149‑1 (MQ=255)
   ggAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATcaca                                                                                                                                                   >  2:1017095/1‑149 (MQ=255)
    gagaATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATAACAg                                                                                                                                                  >  1:3757512/1‑149 (MQ=255)
      gaATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaa                                                                                                                                                <  1:316889/149‑1 (MQ=255)
      gaATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaa                                                                                                                                                <  1:2081812/149‑1 (MQ=255)
       aaTACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGaaa                                                                                                                                               <  1:437267/149‑1 (MQ=255)
         tACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATg                                                                                                                                             <  1:769331/149‑1 (MQ=255)
          aCCAATGGCTATCGCTGTCAGGACCTGAAAGTACAGTCTATTAAACAGAGAGGTTTTCATAGGGTGTCCTATAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGt                                                                                                                                            <  2:3634354/149‑1 (MQ=255)
          aCCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGt                                                                                                                                            <  2:1245027/149‑1 (MQ=255)
          aCCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGt                                                                                                                                            <  2:2527587/149‑1 (MQ=255)
          aCCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGt                                                                                                                                            >  1:1571467/1‑149 (MQ=255)
          aCCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGt                                                                                                                                            <  2:3524312/149‑1 (MQ=255)
            cAATGGCTATCTCTGTCAGGACCTTAAATTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTAc                                                                                                                                          <  2:3265197/149‑1 (MQ=255)
             aaTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTAcc                                                                                                                                         <  1:4048402/149‑1 (MQ=255)
              aTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCt                                                                                                                                        >  2:506657/1‑149 (MQ=255)
              aTGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCt                                                                                                                                        <  1:3065019/149‑1 (MQ=255)
               tGGCTATCGCTGTCAGGACCTGAAAGTAAAGGATTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTg                                                                                                                                       <  1:1238911/149‑1 (MQ=255)
                  cTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGAt                                                                                                                                    <  2:1451785/149‑1 (MQ=255)
                         tGTCAGGACCTTAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGTGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAAtt                                                                                                                             <  2:2214210/149‑1 (MQ=255)
                           tCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGa                                                                                                                           >  1:4681882/1‑149 (MQ=255)
                             aGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGaaa                                                                                                                         <  2:1982589/149‑1 (MQ=255)
                               gACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACAACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACa                                                                                                                       >  2:4739024/1‑149 (MQ=255)
                               gACCTGAAAGTAAAGGCTTTTAAACAGAGAGGATTTCATAGGGTGTCCTTTAGTAAAACCACAGGTAATGTAAGGTGATGGGGTACCAGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACa                                                                                                                       >  1:886140/1‑149 (MQ=255)
                                aCCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaa                                                                                                                      <  1:1806197/149‑1 (MQ=255)
                                aCCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACaa                                                                                                                      <  2:3106567/149‑1 (MQ=255)
                                    gAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAg                                                                                                                  <  2:1289180/149‑1 (MQ=255)
                                    gAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAg                                                                                                                  <  1:394228/149‑1 (MQ=255)
                                     aaaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAgg                                                                                                                 <  1:3309717/149‑1 (MQ=255)
                                      aaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGt                                                                                                                >  1:1891827/1‑149 (MQ=255)
                                      aaGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGt                                                                                                                <  1:4575540/149‑1 (MQ=255)
                                            aGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTaaaaa                                                                                                          <  2:4668387/149‑1 (MQ=255)
                                            aGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTATAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTaaaaa                                                                                                          <  2:1177005/149‑1 (MQ=255)
                                             ggCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAg                                                                                                         >  1:1541428/1‑149 (MQ=255)
                                             ggCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAg                                                                                                         >  1:574310/1‑149 (MQ=255)
                                                ttttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGttt                                                                                                      <  2:4301239/149‑1 (MQ=255)
                                                 tttACACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTg                                                                                                     <  2:2905603/149‑1 (MQ=255)
                                                 tttAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTg                                                                                                     >  2:548073/1‑149 (MQ=255)
                                                   tAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAg                                                                                                   <  2:2592450/149‑1 (MQ=255)
                                                    aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGatc                                                                                                  >  1:168135/1‑147 (MQ=255)
                                                    aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGc                                                                                                  >  1:2168786/1‑149 (MQ=255)
                                                    aaaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGc                                                                                                  <  2:1723017/149‑1 (MQ=255)
                                                     aaCAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCt                                                                                                 <  1:2734794/149‑1 (MQ=255)
                                                        agagagGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGc                                                                                              <  2:3742128/149‑1 (MQ=255)
                                                            agGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGc                                                                                          >  2:3115303/1‑149 (MQ=255)
                                                            agGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGc                                                                                          <  1:3516097/149‑1 (MQ=255)
                                                                  tCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAc                                                                                    <  2:3477703/149‑1 (MQ=255)
                                                                   cATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAAcc                                                                                   >  1:1712340/1‑149 (MQ=255)
                                                                        ggTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATAACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCAc                                                                              <  2:3340681/149‑1 (MQ=255)
                                                                        ggTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCAGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCAc                                                                              >  1:959033/1‑149 (MQ=255)
                                                                          tgtCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtt                                                                            >  1:588219/1‑149 (MQ=255)
                                                                          tgtCCTCTAGTAAAACCACAGGTCGTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACtt                                                                            <  2:4439584/149‑1 (MQ=255)
                                                                           gtCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACttt                                                                           >  2:725707/1‑149 (MQ=255)
                                                                              cTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                                                                        >  1:1001959/1‑149 (MQ=255)
                                                                              cTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                                                                        >  1:2326454/1‑149 (MQ=255)
                                                                              cTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTa                                                                        >  1:3242059/1‑149 (MQ=255)
                                                                                 tAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAg                                                                     <  1:3564168/149‑1 (MQ=255)
                                                                                   gTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtt                                                                   <  1:1364040/149‑1 (MQ=255)
                                                                                   gTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtt                                                                   <  1:3214617/149‑1 (MQ=255)
                                                                                   gTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATAACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtt                                                                   <  2:2506000/149‑1 (MQ=255)
                                                                                     aaaaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGtttt                                                                 <  1:4544233/149‑1 (MQ=255)
                                                                                      aaaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg                                                                >  2:1099016/1‑149 (MQ=255)
                                                                                      aaaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATAACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTg                                                                >  2:4567954/1‑149 (MQ=255)
                                                                                       aaCCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGt                                                               >  2:1490147/1‑149 (MQ=255)
                                                                                           acaGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATc                                                           >  1:1598638/1‑149 (MQ=255)
                                                                                            caGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCa                                                          <  1:1421270/149‑1 (MQ=255)
                                                                                              ggTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGt                                                        <  2:2069794/149‑1 (MQ=255)
                                                                                              ggTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATAACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGt                                                        >  1:2829461/1‑149 (MQ=255)
                                                                                                tCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGttt                                                      <  2:2168670/149‑1 (MQ=255)
                                                                                                 cTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTg                                                     <  2:4544937/149‑1 (MQ=255)
                                                                                                   tGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTggg                                                   <  1:3031066/149‑1 (MQ=255)
                                                                                                       aGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAg                                               >  1:4121495/1‑149 (MQ=255)
                                                                                                          ttATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTa                                            >  1:1932913/1‑149 (MQ=255)
                                                                                                            aTGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACt                                          <  1:3390440/149‑1 (MQ=255)
                                                                                                               gggTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACtttt                                       >  1:1588282/1‑149 (MQ=255)
                                                                                                               gggTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACtttt                                       <  2:1543818/149‑1 (MQ=255)
                                                                                                               gggTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGGAGCTACtttt                                       >  1:3642760/1‑149 (MQ=255)
                                                                                                                   aCCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTttcttc                                   >  1:2994517/1‑149 (MQ=255)
                                                                                                                    ccTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTcc                                  >  1:3131162/1‑149 (MQ=255)
                                                                                                                    ccTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTcc                                  >  2:1859145/1‑149 (MQ=255)
                                                                                                                     cTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCa                                 >  2:4257545/1‑149 (MQ=255)
                                                                                                                        cGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGt                              <  2:3023585/149‑1 (MQ=255)
                                                                                                                          gCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTaa                            <  2:910197/149‑1 (MQ=255)
                                                                                                                             ttAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAActc                         >  2:4787937/1‑149 (MQ=255)
                                                                                                                              tAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAActct                        >  1:587525/1‑149 (MQ=255)
                                                                                                                               aaaaTAACACCCAGACAACATCACAGAAATGTACCTGGATCAGAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAActctc                       <  2:2399537/149‑1 (MQ=255)
                                                                                                                                 aaTAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCtt                     >  2:56920/1‑149 (MQ=255)
                                                                                                                                  aTAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTc                    >  2:3809148/1‑149 (MQ=255)
                                                                                                                                       acCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGa               <  1:1033551/149‑1 (MQ=255)
                                                                                                                                        cccAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGAt              <  1:4315175/149‑1 (MQ=255)
                                                                                                                                         ccAGACAACATAACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGAtt             <  1:2447262/149‑1 (MQ=255)
                                                                                                                                          cAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGAttt            <  2:4195233/149‑1 (MQ=255)
                                                                                                                                           aGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTc           >  2:648014/1‑149 (MQ=255)
                                                                                                                                             acaacaTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCaa         <  2:1915367/149‑1 (MQ=255)
                                                                                                                                                  aTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAggg    >  1:658964/1‑149 (MQ=255)
                                                                                                                                                  aTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAggg    >  1:4071552/1‑149 (MQ=255)
                                                                                                                                                  aTCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAggg    >  1:3646117/1‑149 (MQ=255)
                                                                                                                                                   tCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGa   >  2:77386/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CAAGGAGAATACCAATGGCTATCGCTGTCAGGACCTGAAAGTAAAGGCTTTTAAACAGAGAGGTTTTCATAGGGTGTCCTTTAGTAAAACCACAGGTCTTGTAAGGTTATGGGGTACCTGCGGCCTTAAAATAACACCCAGACAACATCACAGAAATGTACCTGGATCATAATTGAAACAAAAAGGTTAAAAAGTTTGAGCTGGCTCGCACAAACCAGCACTTTTTAAAGTTTTGTAATCAGTTTGGGGTAGCTACTTTTCTTCCAGGTAACTCTCTTCGAAGATTTCAATAGGGAG  >  CP009273/3676738‑3677034

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: