Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,519,188 C→A 40.8% 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,1880CA40.8% 355.0 / inf 311L107L (CTG→CTTnudEadenosine nucleotide hydrolase; substrates include Ap3A, Ap2A, ADP‑ribose, NADH
Reads supporting (aligned to +/- strand):  ref base C (76/108);  new base A (51/76);  total (127/184)
Fisher's exact test for biased strand distribution p-value = 9.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

Only 100 of 314 total aligned reads displayed.

CCATCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCAGGTGAT  >  CP009273/3519040‑3519336
                                                                                                                                                    |                                                                                                                                                    
ccATCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTc                                                                                                                                                      <  1:2867611/149‑1 (MQ=255)
 catcatATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCa                                                                                                                                                     <  2:4686037/149‑1 (MQ=255)
 catcatATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCa                                                                                                                                                     <  2:4706056/149‑1 (MQ=255)
  atcatATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAg                                                                                                                                                    <  1:4321096/149‑1 (MQ=255)
   tcatATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGa                                                                                                                                                   <  2:77765/149‑1 (MQ=255)
    catATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGAt                                                                                                                                                  <  1:1956291/149‑1 (MQ=255)
      tatGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCg                                                                                                                                                <  1:3798303/149‑1 (MQ=255)
      tatGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCg                                                                                                                                                >  2:2829589/1‑149 (MQ=255)
        tGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGtt                                                                                                                                              <  1:242505/149‑1 (MQ=255)
        tGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGtt                                                                                                                                              >  1:447731/1‑149 (MQ=255)
        tGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGtt                                                                                                                                              <  2:3601317/149‑1 (MQ=255)
        tGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGtt                                                                                                                                              >  1:441363/1‑149 (MQ=255)
         gcgcCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTc                                                                                                                                             <  1:4708098/149‑1 (MQ=255)
         gcgcCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTc                                                                                                                                             <  2:1580011/149‑1 (MQ=255)
         gcgcCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTc                                                                                                                                             <  2:1442600/149‑1 (MQ=255)
            ccagcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCt                                                                                                                                          <  1:2331920/149‑1 (MQ=255)
               gcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCa                                                                                                                                       >  2:2036470/1‑149 (MQ=255)
               gcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCa                                                                                                                                       >  2:1104462/1‑149 (MQ=255)
                 ggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCaaa                                                                                                                                     <  2:604410/149‑1 (MQ=255)
                 ggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCaaa                                                                                                                                     >  2:1071152/1‑149 (MQ=255)
                  gccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAAt                                                                                                                                    <  2:4783268/149‑1 (MQ=255)
                     agcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCa                                                                                                                                 <  2:556265/149‑1 (MQ=255)
                     agcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCa                                                                                                                                 >  2:2203747/1‑149 (MQ=255)
                      gcgcACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATccca                                                                                                                                >  1:1256812/1‑147 (MQ=255)
                         cACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACct                                                                                                                             >  1:139516/1‑149 (MQ=255)
                         cACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACct                                                                                                                             <  2:912330/149‑1 (MQ=255)
                          aCCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACctc                                                                                                                            <  2:2750003/149‑1 (MQ=255)
                           ccTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACctct                                                                                                                           <  1:2798414/149‑1 (MQ=255)
                            cTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCtctt                                                                                                                          >  2:754204/1‑149 (MQ=255)
                              gtgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCtcttct                                                                                                                        <  2:1385386/149‑1 (MQ=255)
                                ggTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCttt                                                                                                                      <  1:4438776/149‑1 (MQ=255)
                                      ggCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTc                                                                                                                <  1:2982303/149‑1 (MQ=255)
                                      ggCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTc                                                                                                                >  2:3808982/1‑149 (MQ=255)
                                       gCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcg                                                                                                               >  1:1345615/1‑149 (MQ=255)
                                                  tCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGc                                                                                                    >  2:1235459/1‑149 (MQ=255)
                                                  tCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGc                                                                                                    <  1:2069969/149‑1 (MQ=255)
                                                   cGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCt                                                                                                   <  2:3481237/149‑1 (MQ=255)
                                                    gCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCtt                                                                                                  >  2:327784/1‑149 (MQ=255)
                                                     ccTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTc                                                                                                 <  1:973314/149‑1 (MQ=255)
                                                        tCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTa                                                                                              <  1:1011507/149‑1 (MQ=255)
                                                         ccAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAg                                                                                             >  2:1505458/1‑149 (MQ=255)
                                                         ccAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAg                                                                                             >  2:3310085/1‑149 (MQ=255)
                                                            gTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACg                                                                                          >  2:1797789/1‑149 (MQ=255)
                                                            gTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACg                                                                                          >  1:2042895/1‑149 (MQ=255)
                                                               aTTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTCGCTCGCGGTTAGCGGCTTCGTAGACGCtt                                                                                       >  2:2477980/1‑149 (MQ=255)
                                                                ttCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCttt                                                                                      <  1:1902659/149‑1 (MQ=255)
                                                                ttCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCttt                                                                                      <  2:1622695/149‑1 (MQ=255)
                                                                  ccGGGTAGAGATCGTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCa                                                                                    <  2:4303948/149‑1 (MQ=255)
                                                                     ggTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCAccc                                                                                 >  2:1540942/1‑149 (MQ=255)
                                                                         gagaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGAt                                                                             <  2:1970843/149‑1 (MQ=255)
                                                                          agaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATc                                                                            >  1:2640008/1‑149 (MQ=255)
                                                                             tCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAAt                                                                         >  2:2993099/1‑149 (MQ=255)
                                                                               ttGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTa                                                                       >  1:3758312/1‑149 (MQ=255)
                                                                                tGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTaa                                                                      <  1:327784/149‑1 (MQ=255)
                                                                                tGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTaa                                                                      <  2:1741781/149‑1 (MQ=255)
                                                                                 gcgcTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAAt                                                                     >  1:1499055/1‑149 (MQ=255)
                                                                                    cTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCt                                                                  <  2:3786100/149‑1 (MQ=255)
                                                                                     tACCACGATATTCATTTTGCTTGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTGCAAATCCAACCTCTTCTATTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCtt                                                                 <  2:1074394/149‑1 (MQ=255)
                                                                                     tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCtt                                                                 >  1:346934/1‑149 (MQ=255)
                                                                                     tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCtt                                                                 <  2:3372186/149‑1 (MQ=255)
                                                                                     tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCtt                                                                 <  2:4058335/149‑1 (MQ=255)
                                                                                     tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCtt                                                                 <  2:3160322/149‑1 (MQ=255)
                                                                                     tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCtt                                                                 <  1:3080375/149‑1 (MQ=255)
                                                                                      aCCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                                >  2:858801/1‑149 (MQ=255)
                                                                                      aCCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                                >  2:1048816/1‑149 (MQ=255)
                                                                                      aCCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                                >  1:3521/1‑149 (MQ=255)
                                                                                      aCCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                                <  1:3767393/149‑1 (MQ=255)
                                                                                        cACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCg                                                              >  2:2696130/1‑149 (MQ=255)
                                                                                        cACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCg                                                              >  1:4774382/1‑149 (MQ=255)
                                                                                          cGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaa                                                            <  2:3538836/149‑1 (MQ=255)
                                                                                             tatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaaaaa                                                         >  1:271682/1‑149 (MQ=255)
                                                                                               ttCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAcc                                                       >  1:2103336/1‑149 (MQ=255)
                                                                                                tCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCATCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCt                                                      <  1:3961288/149‑1 (MQ=255)
                                                                                                   ttttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAAt                                                   >  1:2548773/1‑149 (MQ=255)
                                                                                                   ttttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAAt                                                   <  1:3670390/149‑1 (MQ=255)
                                                                                                    tgtGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAAtt                                                  <  1:159847/147‑1 (MQ=255)
                                                                                                       gCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGt                                               <  2:1308642/149‑1 (MQ=255)
                                                                                                       gCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGt                                               <  2:4109091/149‑1 (MQ=255)
                                                                                                        cTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAATTCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTa                                              <  2:2518993/149‑1 (MQ=255)
                                                                                                           taaaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGa                                           <  2:3695682/148‑1 (MQ=255)
                                                                                                            aaaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGAt                                          >  1:648084/1‑149 (MQ=255)
                                                                                                            aaaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGAt                                          <  2:2110360/149‑1 (MQ=255)
                                                                                                              aaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTc                                        >  1:4111986/1‑149 (MQ=255)
                                                                                                              aaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTc                                        <  2:2304333/149‑1 (MQ=255)
                                                                                                               aGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCa                                       >  2:2765548/1‑149 (MQ=255)
                                                                                                               aGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCa                                       <  2:441363/149‑1 (MQ=255)
                                                                                                                   aGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTc                                   >  1:904208/1‑149 (MQ=255)
                                                                                                                   aGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTc                                   >  1:1868747/1‑149 (MQ=255)
                                                                                                                    gACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTcc                                  >  2:1702481/1‑149 (MQ=255)
                                                                                                                     aCGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTccc                                 <  1:2563473/149‑1 (MQ=255)
                                                                                                                      cGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCa                                >  1:1534799/1‑149 (MQ=255)
                                                                                                                       ggTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCAc                               <  1:741317/149‑1 (MQ=255)
                                                                                                                             aTGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTa                         <  1:3480113/149‑1 (MQ=255)
                                                                                                                               gCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTAtt                       <  1:1556325/149‑1 (MQ=255)
                                                                                                                                   aGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTcgcg                   >  2:4750064/1‑149 (MQ=255)
                                                                                                                                      ttcttcAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggat                >  2:2402406/1‑149 (MQ=255)
                                                                                                                                      ttcttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggat                >  2:4531911/1‑149 (MQ=255)
                                                                                                                                           cAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcagga           <  2:3919878/149‑1 (MQ=255)
                                                                                                                                               aaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGgatcaggatcag       >  1:1741295/1‑149 (MQ=255)
                                                                                                                                                aaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGatcaggatcagg      >  1:1177044/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCATCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCAGGTGAT  >  CP009273/3519040‑3519336

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: