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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,881,6450CA97.8% 573.3 / 3.3 184D202Y (GAT→TAT) rndribonuclease D
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (89/91);  total (90/94)
Fisher's exact test for biased strand distribution p-value = 6.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.40e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 186 total aligned reads displayed.

TAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCGC  >  CP009273/1881500‑1881789
                                                                                                                                                 |                                                                                                                                                
taACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAAc                                                                                                                                                >  1:2598387/1‑148 (MQ=255)
  aCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGc                                                                                                                                              >  2:1580475/1‑148 (MQ=255)
  aCGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGc                                                                                                                                              <  2:2050846/148‑1 (MQ=255)
   cgcgCTACCGACCACAAATGCTCTTCACGTACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                             >  1:1893748/1‑148 (MQ=255)
   cgcgCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                             >  1:2884294/1‑148 (MQ=255)
   cgcgCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                             >  1:1875177/1‑148 (MQ=255)
   cgcgCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                             >  1:1893942/1‑148 (MQ=255)
   cgcgCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGcc                                                                                                                                             >  1:118439/1‑148 (MQ=255)
       cTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGc                                                                                                                                         <  1:1060023/148‑1 (MQ=255)
          ccGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATc                                                                                                                                      >  2:1904077/1‑148 (MQ=255)
          ccGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATc                                                                                                                                      >  2:615254/1‑148 (MQ=255)
           cGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCt                                                                                                                                     <  2:1726192/148‑1 (MQ=255)
           cGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCt                                                                                                                                     >  1:1476929/1‑148 (MQ=255)
             aCCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTc                                                                                                                                   >  2:1172413/1‑148 (MQ=255)
               cacaAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCg                                                                                                                                 >  2:326223/1‑148 (MQ=255)
                 caAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgc                                                                                                                               >  1:881404/1‑148 (MQ=255)
                  aaaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgcg                                                                                                                              <  2:2679377/148‑1 (MQ=255)
                  aaaTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGgcg                                                                                                                              >  2:1654365/1‑148 (MQ=255)
                    aTGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCa                                                                                                                            >  2:1703591/1‑148 (MQ=255)
                     tGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                           >  2:3409551/1‑148 (MQ=255)
                     tGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                           >  2:3408907/1‑148 (MQ=255)
                     tGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                           >  2:3342490/1‑148 (MQ=255)
                     tGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCaa                                                                                                                           >  2:2249894/1‑148 (MQ=255)
                      gCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAAc                                                                                                                          >  2:2095399/1‑148 (MQ=255)
                        tctTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGa                                                                                                                        >  1:669045/1‑148 (MQ=255)
                           tCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACtt                                                                                                                     >  1:95307/1‑148 (MQ=255)
                            cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTc                                                                                                                    <  2:1655425/148‑1 (MQ=255)
                            cacgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTc                                                                                                                    <  2:366569/148‑1 (MQ=255)
                             acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                   <  2:3061708/148‑1 (MQ=255)
                             acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                   >  2:2489202/1‑148 (MQ=255)
                             acgcacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTcc                                                                                                                   <  2:206159/148‑1 (MQ=255)
                                cacgACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTga                                                                                                                <  1:2553083/148‑1 (MQ=255)
                                     cAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGacgacg                                                                                                           >  1:1100753/1‑148 (MQ=255)
                                        aGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAt                                                                                                        >  2:3274630/1‑148 (MQ=255)
                                          ttCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCAttt                                                                                                      >  1:3324851/1‑148 (MQ=255)
                                             aCCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCa                                                                                                   <  1:3110405/148‑1 (MQ=255)
                                             aCCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCa                                                                                                   <  1:1363178/148‑1 (MQ=255)
                                                     aTCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAggcgg                                                                                           >  2:301440/1‑148 (MQ=255)
                                                       cgcgCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCa                                                                                         <  2:426557/148‑1 (MQ=255)
                                                       cgcgCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCa                                                                                         <  2:1996010/148‑1 (MQ=255)
                                                            tctcGCGCCTTTCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCAt                                                                                    <  2:920293/148‑1 (MQ=255)
                                                            tctcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCAt                                                                                    >  2:1767844/1‑148 (MQ=255)
                                                             ctcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGCCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATc                                                                                   <  2:1833666/148‑1 (MQ=255)
                                                              tcGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATcc                                                                                  >  2:2941607/1‑148 (MQ=255)
                                                               cgcgCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCa                                                                                 >  1:2169144/1‑148 (MQ=255)
                                                                   ccTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAgcgc                                                                             <  2:2884294/148‑1 (MQ=255)
                                                                     ttGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccg                                                                           >  1:402252/1‑148 (MQ=255)
                                                                      tGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGcgccgc                                                                          <  1:505069/148‑1 (MQ=255)
                                                                       gcgcAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCa                                                                         <  1:1435809/148‑1 (MQ=255)
                                                                          cAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGt                                                                      >  1:1856/1‑148 (MQ=255)
                                                                          cAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGt                                                                      <  2:1662998/148‑1 (MQ=255)
                                                                          cAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGt                                                                      >  2:115384/1‑148 (MQ=255)
                                                                             tCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTagc                                                                   <  1:1088376/148‑1 (MQ=255)
                                                                                ccAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCTGCATTCATCCAGCGCCGCAGGTagccag                                                                >  1:1280835/1‑148 (MQ=255)
                                                                                     cGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGa                                                           <  2:1512210/148‑1 (MQ=255)
                                                                                      ggCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAg                                                          <  1:2516198/148‑1 (MQ=255)
                                                                                           aCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTc                                                     >  2:1078428/1‑148 (MQ=255)
                                                                                           aCAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTc                                                     >  2:838331/1‑148 (MQ=255)
                                                                                             aGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCg                                                   <  1:1554035/148‑1 (MQ=255)
                                                                                               ttGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGtt                                                 >  1:1653835/1‑148 (MQ=255)
                                                                                                     aCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTAcc                                           <  1:1235664/148‑1 (MQ=255)
                                                                                                     aCAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTAcc                                           <  1:1235752/148‑1 (MQ=255)
                                                                                                       aGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCAt                                         >  2:252768/1‑148 (MQ=255)
                                                                                                        ggCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATa                                        >  1:3438104/1‑148 (MQ=255)
                                                                                                         gCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATaa                                       >  1:1426670/1‑148 (MQ=255)
                                                                                                             gTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCtt                                   <  1:2095399/148‑1 (MQ=255)
                                                                                                                 gCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCg                               <  1:2249894/148‑1 (MQ=255)
                                                                                                                 gCGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCg                               <  1:3342490/148‑1 (MQ=255)
                                                                                                                  cGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCgg                              <  2:1100753/148‑1 (MQ=255)
                                                                                                                  cGTGTGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCgg                              <  2:1738440/148‑1 (MQ=255)
                                                                                                                   gtgtGCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGt                             <  1:1703587/148‑1 (MQ=255)
                                                                                                                    tgtgCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTg                            <  1:2697188/148‑1 (MQ=255)
                                                                                                                    tgtgCGTAAGTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTg                            <  1:1958401/148‑1 (MQ=255)
                                                                                                                      tgCGTTATTGCAAGGCATTGGTGATATAACGCCAGTCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGAt                          <  2:1947147/148‑1 (MQ=255)
                                                                                                                       tcGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATc                         >  1:1035520/2‑148 (MQ=255)
                                                                                                                       gCGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATc                         <  2:3438104/148‑1 (MQ=255)
                                                                                                                        cGTAATTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCg                        <  1:269672/148‑1 (MQ=255)
                                                                                                                           aaTTGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTa                     <  2:3325334/148‑1 (MQ=255)
                                                                                                                             ttGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAAc                   >  2:1974153/1‑148 (MQ=255)
                                                                                                                              tGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACa                  >  2:1697991/1‑148 (MQ=255)
                                                                                                                              tGCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACa                  >  2:1082850/1‑148 (MQ=255)
                                                                                                                               gCCAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAg                 <  1:2489202/148‑1 (MQ=255)
                                                                                                                                ccAGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGa                >  2:2137613/1‑148 (MQ=255)
                                                                                                                                  aGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa              >  1:3381593/1‑148 (MQ=255)
                                                                                                                                  aGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa              <  2:2598387/148‑1 (MQ=255)
                                                                                                                                  aGGCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATa              <  1:1198011/148‑1 (MQ=255)
                                                                                                                                   ggCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAc             >  1:1582397/1‑148 (MQ=255)
                                                                                                                                   ggCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAc             <  2:1856/148‑1 (MQ=255)
                                                                                                                                    gCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAcc            <  1:379798/148‑1 (MQ=255)
                                                                                                                                    gCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAcc            <  1:2928260/148‑1 (MQ=255)
                                                                                                                                    gCATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATAcc            <  1:2008818/148‑1 (MQ=255)
                                                                                                                                     cATTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCa           >  2:2773465/1‑148 (MQ=255)
                                                                                                                                      aTTGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAg          >  1:2957854/1‑148 (MQ=255)
                                                                                                                                       ttGGTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGa         <  1:985769/148‑1 (MQ=255)
                                                                                                                                          gTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACAt      <  1:1767844/148‑1 (MQ=255)
                                                                                                                                          gTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACAt      >  1:1023054/1‑148 (MQ=255)
                                                                                                                                          gTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACAt      <  2:1396626/148‑1 (MQ=255)
                                                                                                                                          gTGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACAt      <  2:2292215/148‑1 (MQ=255)
                                                                                                                                           tGATATAACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATc     >  2:2906711/1‑148 (MQ=255)
                                                                                                                                              tataACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCTGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATccgc  <  2:2443454/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TAACGCGCTACCGACCACAAATGCTCTTCACGCACGACAAAGTTCACCGCCAGATCGCGCTCTCGCGCCTTGCGCAGTCGCCAGTCGGCTAACAGTTGCAGACAGGCCAGTTGGCGTGTGCGTAATTGCCAGGCATTGGTGATATCACGCCAGGCATCTTCCGGCGCAACGACTTCCTGACGACGCATTTGCATCAGGCGGCATTCATCCAGCGCCGCAGGTAGCCAGCCGGAGGCCTCCGTTTCTACCATAAGCTTGGCGGTGATCGGTAACAGATACCAGACATCCGC  >  CP009273/1881500‑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: