Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,519,188 C→A 87.2% L107L (CTG→CTT nudE ← adenosine nucleotide hydrolase; substrates include Ap3A, Ap2A, ADP‑ribose, NADH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,519,1880CA87.2% 497.4 / 78.7 220L107L (CTG→CTTnudEadenosine nucleotide hydrolase; substrates include Ap3A, Ap2A, ADP‑ribose, NADH
Reads supporting (aligned to +/- strand):  ref base C (14/14);  new base A (87/105);  total (101/119)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.77e-01

Only 100 of 220 total aligned reads displayed.

TCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCAGGTG  >  CP009273/3519043‑3519334
                                                                                                                                                 |                                                                                                                                                  
tcatATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAg                                                                                                                                                  <  1:165083/148‑1 (MQ=255)
 catATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGa                                                                                                                                                 >  1:900796/1‑148 (MQ=255)
 catATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGa                                                                                                                                                 <  1:196837/148‑1 (MQ=255)
    atGCGCCCGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCg                                                                                                                                              <  2:856039/148‑1 (MQ=255)
     tGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGt                                                                                                                                             <  2:1088292/148‑1 (MQ=255)
     tGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGt                                                                                                                                             <  2:223222/148‑1 (MQ=255)
      gcgcCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTCGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGtt                                                                                                                                            >  1:3168716/1‑148 (MQ=255)
         ccagcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGc                                                                                                                                         >  1:1624603/1‑148 (MQ=255)
           agcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTGCAAAAAAGTAAGATCGTTCGCTc                                                                                                                                       >  2:3292694/1‑148 (MQ=255)
            gcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTcc                                                                                                                                      >  2:689233/1‑148 (MQ=255)
             cggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCa                                                                                                                                     >  2:482601/1‑148 (MQ=255)
               gccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCaaa                                                                                                                                   >  2:2624281/1‑148 (MQ=255)
                 cagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATc                                                                                                                                 >  1:1821320/1‑148 (MQ=255)
                   gcgcACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCa                                                                                                                               >  1:2110024/1‑148 (MQ=255)
                    cgcACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGGTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCaa                                                                                                                              <  2:1273627/148‑1 (MQ=255)
                       aCCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACct                                                                                                                           <  1:456156/148‑1 (MQ=255)
                       aCCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACct                                                                                                                           >  1:3106114/1‑148 (MQ=255)
                       aCCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACct                                                                                                                           >  2:1756643/1‑148 (MQ=255)
                        ccTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACctc                                                                                                                          >  2:2172572/1‑148 (MQ=255)
                         cTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACctct                                                                                                                         >  1:795128/1‑148 (MQ=255)
                         cTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACctct                                                                                                                         >  2:1853895/1‑148 (MQ=255)
                           gtgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCtcttc                                                                                                                       <  1:3245927/148‑1 (MQ=255)
                             ggTAGCGGCTCTGGCTCATCGCCTTCTAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTcttctt                                                                                                                     <  1:374050/148‑1 (MQ=255)
                             ggTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTcttctt                                                                                                                     >  1:2294028/1‑148 (MQ=255)
                              gTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGTAATGTAAGACGGTGCCATGCTGAGGTTCTTCAATATAGTAAGATCGTTCGCTGCAAATCCAACCTGTTCttt                                                                                                                    >  1:656948/1‑148 (MQ=255)
                                aGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTa                                                                                                                  <  1:3050462/148‑1 (MQ=255)
                                aGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTa                                                                                                                  <  1:2160501/148‑1 (MQ=255)
                                aGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTa                                                                                                                  <  2:2330553/148‑1 (MQ=255)
                                 gCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAg                                                                                                                 <  1:20309/148‑1 (MQ=255)
                                    gCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTc                                                                                                              >  2:1547966/1‑148 (MQ=255)
                                     ctctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcg                                                                                                             >  1:2150135/1‑148 (MQ=255)
                                     ctctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcg                                                                                                             >  2:2521134/1‑148 (MQ=255)
                                     ctctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcg                                                                                                             >  1:487594/1‑148 (MQ=255)
                                     ctctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcg                                                                                                             <  2:2113386/148‑1 (MQ=255)
                                       ctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcgcg                                                                                                           >  2:899396/1‑148 (MQ=255)
                                          gCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGtt                                                                                                        >  2:607521/1‑148 (MQ=255)
                                          gCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGtt                                                                                                        <  1:1291847/148‑1 (MQ=255)
                                          gCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGtt                                                                                                        <  1:746074/148‑1 (MQ=255)
                                               tCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCgg                                                                                                   >  2:511929/1‑148 (MQ=255)
                                                cGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGc                                                                                                  <  1:2197209/148‑1 (MQ=255)
                                                cGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGc                                                                                                  >  2:620706/1‑148 (MQ=255)
                                                 gCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCt                                                                                                 <  1:2520262/148‑1 (MQ=255)
                                                     tCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGt                                                                                             <  2:2600055/148‑1 (MQ=255)
                                                      ccAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTa                                                                                            >  1:1382284/1‑148 (MQ=255)
                                                              tCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCttt                                                                                    >  2:3048652/1‑148 (MQ=255)
                                                                  ggTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCAcc                                                                                <  2:862634/148‑1 (MQ=255)
                                                                        gaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATc                                                                          >  1:1288814/1‑148 (MQ=255)
                                                                         aTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCa                                                                         <  1:979252/148‑1 (MQ=255)
                                                                           cTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAAt                                                                       <  1:1853895/148‑1 (MQ=255)
                                                                                gcTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATc                                                                  <  1:636675/148‑1 (MQ=255)
                                                                                 cTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATcc                                                                 <  2:669771/148‑1 (MQ=255)
                                                                                 cTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATcc                                                                 >  1:2749135/1‑148 (MQ=255)
                                                                                 cTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATcc                                                                 <  2:56081/148‑1 (MQ=255)
                                                                                  tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCt                                                                >  2:1491332/1‑148 (MQ=255)
                                                                                  tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCt                                                                <  2:47088/148‑1 (MQ=255)
                                                                                   aCCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCtt                                                               <  2:2190469/148‑1 (MQ=255)
                                                                                    ccACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                              >  2:253845/1‑148 (MQ=255)
                                                                                     cACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTc                                                             >  1:812727/1‑148 (MQ=255)
                                                                                      aCGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCg                                                            >  2:2929431/1‑148 (MQ=255)
                                                                                       cGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGa                                                           >  2:150447/1‑148 (MQ=255)
                                                                                         atatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaaa                                                         <  1:687966/148‑1 (MQ=255)
                                                                                          tatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaaaa                                                        >  2:1439764/1‑148 (MQ=255)
                                                                                            ttCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAc                                                      <  2:3168716/148‑1 (MQ=255)
                                                                                            ttCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAc                                                      <  1:1074118/148‑1 (MQ=255)
                                                                                             tCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAcc                                                     >  2:2681098/1‑148 (MQ=255)
                                                                                             tCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAcc                                                     <  1:1120984/148‑1 (MQ=255)
                                                                                                ttttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTaa                                                  >  1:2925237/1‑148 (MQ=255)
                                                                                                ttttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTaa                                                  >  1:1323639/1‑148 (MQ=255)
                                                                                                ttttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTaa                                                  >  1:95933/1‑148 (MQ=255)
                                                                                                 tttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAAt                                                 <  2:20999/148‑1 (MQ=255)
                                                                                                 tttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCTAAAAACCTAAt                                                 >  2:2164942/1‑148 (MQ=255)
                                                                                                  ttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAAtt                                                <  2:2514018/148‑1 (MQ=255)
                                                                                                  ttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAAtt                                                <  2:1947833/148‑1 (MQ=255)
                                                                                                     cTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGt                                             <  1:633970/148‑1 (MQ=255)
                                                                                                     cTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGt                                             <  2:2890432/148‑1 (MQ=255)
                                                                                                      tGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTa                                            >  1:1448545/1‑148 (MQ=255)
                                                                                                         aaaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGa                                         <  2:2826971/148‑1 (MQ=255)
                                                                                                          aaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGAt                                        >  2:1578358/1‑148 (MQ=255)
                                                                                                           aaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGAtt                                       >  2:1998191/1‑148 (MQ=255)
                                                                                                           aaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGAtt                                       <  1:2172572/148‑1 (MQ=255)
                                                                                                            aGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCAAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTc                                      >  2:2273279/1‑148 (MQ=255)
                                                                                                              tAAGACGGTGCCATGCTGAGCCTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAg                                    <  1:119314/148‑1 (MQ=255)
                                                                                                                aGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGtt                                  >  2:3286417/1‑148 (MQ=255)
                                                                                                                 gACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTc                                 <  2:1884150/148‑1 (MQ=255)
                                                                                                                 gACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTc                                 <  2:1692626/148‑1 (MQ=255)
                                                                                                                    ggTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCa                              <  1:564398/148‑1 (MQ=255)
                                                                                                                         cATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCg                         <  2:910873/148‑1 (MQ=255)
                                                                                                                         cATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAAGTCGTAGGATTCAGTTCCCACTGCg                         <  2:2565619/148‑1 (MQ=255)
                                                                                                                          aTGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGt                        <  2:294804/148‑1 (MQ=255)
                                                                                                                          aTGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGt                        >  2:1819800/1‑148 (MQ=255)
                                                                                                                            gCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTAt                      <  2:2702714/148‑1 (MQ=255)
                                                                                                                            gCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTAt                      <  1:2658687/148‑1 (MQ=255)
                                                                                                                                  cttcttCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCgg                >  2:2138443/1‑148 (MQ=255)
                                                                                                                                  cttcttCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCgg                >  1:583556/1‑148 (MQ=255)
                                                                                                                                   ttcttcAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCgga               <  1:511929/148‑1 (MQ=255)
                                                                                                                                   ttcttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCgga               >  2:1311206/1‑148 (MQ=255)
                                                                                                                                    tcttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggat              <  1:2709409/148‑1 (MQ=255)
                                                                                                                                    tcttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATACAGTTCCCACTGCGTATTCGCggat              >  2:274517/1‑148 (MQ=255)
                                                                                                                                        cAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcagg          <  1:2531428/148‑1 (MQ=255)
                                                                                                                                           aaaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcaggatc       <  2:812727/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCAGGTG  >  CP009273/3519043‑3519334

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: