Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,881,645 C→A 82.3% D202Y (GAT→TAT)  rnd ← ribonuclease D

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,881,6450CA82.3% 286.5 / 82.4 153D202Y (GAT→TAT) rndribonuclease D
Reads supporting (aligned to +/- strand):  ref base C (17/10);  new base A (60/66);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 2.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.04e-01

Only 100 of 154 total aligned reads displayed.

ACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCGC  >  CP009273/1881502‑1881789
                                                                                                                                               |                                                                                                                                                
aCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGc                                                                                                                                              <  2:420943/148‑1 (MQ=255)
aCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGc                                                                                                                                              <  1:14849/148‑1 (MQ=255)
 cgcgCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                             >  2:1579588/1‑148 (MQ=255)
   cgcTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAg                                                                                                                                           >  2:1954431/1‑148 (MQ=255)
     cTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGc                                                                                                                                         >  1:2900546/1‑148 (MQ=255)
      tACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCa                                                                                                                                        <  1:2782899/148‑1 (MQ=255)
      tACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCa                                                                                                                                        <  2:174518/148‑1 (MQ=255)
       aCCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCAt                                                                                                                                       <  2:1940167/148‑1 (MQ=255)
        ccGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATc                                                                                                                                      <  2:2332251/148‑1 (MQ=255)
        ccGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATc                                                                                                                                      >  2:350286/1‑148 (MQ=255)
        ccGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATc                                                                                                                                      >  2:350290/1‑148 (MQ=255)
         cGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCt                                                                                                                                     >  1:2568773/1‑148 (MQ=255)
         cGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCt                                                                                                                                     >  1:2090393/1‑148 (MQ=255)
          gACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCtt                                                                                                                                    <  2:2128203/148‑1 (MQ=255)
          gACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCtt                                                                                                                                    <  2:2315898/148‑1 (MQ=255)
           aCCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTc                                                                                                                                   >  2:2043120/1‑148 (MQ=255)
            ccACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTcc                                                                                                                                  >  1:2785036/1‑148 (MQ=255)
            ccACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTcc                                                                                                                                  <  1:2217418/148‑1 (MQ=255)
            ccACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTcc                                                                                                                                  >  2:1661969/1‑148 (MQ=255)
              acaAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCgg                                                                                                                                >  1:726909/1‑148 (MQ=255)
               caAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgc                                                                                                                               >  1:1627856/1‑148 (MQ=255)
                aaaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGAGATAACGCCAGGCATCTTCCGgcg                                                                                                                              >  2:1201201/1‑148 (MQ=255)
                  aTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCa                                                                                                                            >  2:2856965/1‑148 (MQ=255)
                  aTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCa                                                                                                                            >  2:2857813/1‑148 (MQ=255)
                    gcactTCACGCACGACAAAGTTCACCGCCCGATCGCGCTCTCGCGCCTTGCGCATTCGCCAGACGGCTAACAGTTGAAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGTCATAGGTGATATCACGCCAGGCATCTTCCGGCGCAAc                                                                                                                          >  1:2563425/4‑148 (MQ=255)
                    gCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAAc                                                                                                                          >  1:909825/1‑148 (MQ=255)
                    gCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAAc                                                                                                                          >  2:2570773/1‑148 (MQ=255)
                     ctctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACg                                                                                                                         >  1:2703053/1‑148 (MQ=255)
                       ctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGAc                                                                                                                       >  2:3046675/1‑148 (MQ=255)
                       ctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGAc                                                                                                                       <  2:656301/148‑1 (MQ=255)
                         tCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGNCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACtt                                                                                                                     <  1:1678749/148‑1 (MQ=255)
                          cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTc                                                                                                                    <  1:118743/148‑1 (MQ=255)
                          cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTc                                                                                                                    <  1:118741/148‑1 (MQ=255)
                           acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                   <  2:1908109/148‑1 (MQ=255)
                           acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                   <  2:51484/148‑1 (MQ=255)
                             gcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTg                                                                                                                 <  1:1063458/148‑1 (MQ=255)
                              cacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTga                                                                                                                <  1:209928/148‑1 (MQ=255)
                              cacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTga                                                                                                                <  1:2519917/148‑1 (MQ=255)
                                cgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacg                                                                                                              <  1:3297190/148‑1 (MQ=255)
                                cgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacg                                                                                                              <  1:1661510/148‑1 (MQ=255)
                                 gACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacga                                                                                                             <  1:1873065/148‑1 (MQ=255)
                                  aCAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacgac                                                                                                            >  2:2422833/1‑148 (MQ=255)
                                  aCAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTgacgac                                                                                                            >  2:2422212/1‑148 (MQ=255)
                                     aaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCa                                                                                                         >  2:1633481/1‑148 (MQ=255)
                                     aaGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCa                                                                                                         >  2:3092016/1‑148 (MQ=255)
                                      aGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAt                                                                                                        <  2:631432/148‑1 (MQ=255)
                                      aGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAt                                                                                                        <  2:631445/148‑1 (MQ=255)
                                        ttCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAttt                                                                                                      >  1:1179385/1‑148 (MQ=255)
                                            ccGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCAt                                                                                                  <  1:1391107/148‑1 (MQ=255)
                                                 aGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCCTTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggc                                                                                             >  2:1532494/1‑148 (MQ=255)
                                                 aGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggc                                                                                             >  2:2763065/1‑148 (MQ=255)
                                                  gATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcg                                                                                            >  1:3138787/1‑148 (MQ=255)
                                                  gATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcg                                                                                            >  1:3136213/1‑148 (MQ=255)
                                                  gATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcg                                                                                            >  2:1353118/1‑148 (MQ=255)
                                                         ctctCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCa                                                                                     <  1:1416570/148‑1 (MQ=255)
                                                          tctcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCAt                                                                                    <  1:237963/148‑1 (MQ=255)
                                                          tctcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCAt                                                                                    <  1:1056880/148‑1 (MQ=255)
                                                            tcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATcc                                                                                  <  2:3161200/148‑1 (MQ=255)
                                                                   ttGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                           <  2:2562974/148‑1 (MQ=255)
                                                                   ttGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                           <  2:726909/148‑1 (MQ=255)
                                                                   ttGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                           <  2:2840412/148‑1 (MQ=255)
                                                                       gcAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAgg                                                                       <  2:900304/148‑1 (MQ=255)
                                                                        cAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGt                                                                      <  1:3281039/148‑1 (MQ=255)
                                                                        cAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGt                                                                      <  2:2806588/148‑1 (MQ=255)
                                                                        cAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGt                                                                      <  2:2806166/148‑1 (MQ=255)
                                                                           tCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagc                                                                   <  1:2570157/148‑1 (MQ=255)
                                                                           tCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagc                                                                   >  1:456634/1‑148 (MQ=255)
                                                                                      cTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGc                                                        <  1:2326189/148‑1 (MQ=255)
                                                                                       tAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGcc                                                       >  1:2310579/1‑148 (MQ=255)
                                                                                        aaCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCt                                                      <  1:3092016/148‑1 (MQ=255)
                                                                                         aCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTc                                                     >  1:2513526/1‑148 (MQ=255)
                                                                                                       gCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATaa                                       <  2:2937010/148‑1 (MQ=255)
                                                                                                       gCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATaa                                       <  2:1288331/148‑1 (MQ=255)
                                                                                                       gCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATaa                                       <  1:3201463/148‑1 (MQ=255)
                                                                                                         cAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGc                                     <  1:2006300/148‑1 (MQ=255)
                                                                                                           gTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCtt                                   <  1:2837848/148‑1 (MQ=255)
                                                                                                           gTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCtt                                   <  1:1001860/148‑1 (MQ=255)
                                                                                                            ttGGCGTGTTGGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTg                                  <  2:2681310/148‑1 (MQ=255)
                                                                                                            ttGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTg                                  <  2:2309124/148‑1 (MQ=255)
                                                                                                            ttGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTg                                  <  2:2236878/148‑1 (MQ=255)
                                                                                                               gCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCg                               >  1:1246241/1‑148 (MQ=255)
                                                                                                               gCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCg                               >  1:486295/1‑148 (MQ=255)
                                                                                                                   gtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATCAGCTTGGCGGTGa                           >  2:1161441/1‑148 (MQ=255)
                                                                                                                   gtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGa                           <  1:992278/148‑1 (MQ=255)
                                                                                                                   gtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGa                           >  2:3096972/1‑148 (MQ=255)
                                                                                                                   gtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGACAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGa                           <  1:991204/148‑1 (MQ=255)
                                                                                                                     gCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATc                         >  1:2880351/1‑148 (MQ=255)
                                                                                                                        tAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGt                      <  1:3144621/148‑1 (MQ=255)
                                                                                                                         aaTTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTa                     >  1:1931739/1‑148 (MQ=255)
                                                                                                                          aTTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTaa                    >  1:2552367/1‑148 (MQ=255)
                                                                                                                          aTTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTaa                    >  1:2552408/1‑148 (MQ=255)
                                                                                                                           ttGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAAc                   >  1:81536/1‑148 (MQ=255)
                                                                                                                           ttGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAAc                   >  1:2946917/1‑148 (MQ=255)
                                                                                                                                aGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa              <  2:3138787/148‑1 (MQ=255)
                                                                                                                                aGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa              >  1:1885554/1‑148 (MQ=255)
                                                                                                                                 ggCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAc             <  2:3115447/148‑1 (MQ=255)
                                                                                                                                     ttGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGa         >  2:2034796/1‑148 (MQ=255)
                                                                                                                                     ttGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGa         >  1:142605/1‑148 (MQ=255)
                                                                                                                                         tGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATc     <  1:1150159/148‑1 (MQ=255)
                                                                                                                                          gATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATcc    >  1:2162868/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                
ACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCGC  >  CP009273/1881502‑1881789

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