Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,519,188 C→A 100% 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,1880CA99.2% 805.8 / ‑2.0 254L107L (CTG→CTTnudEadenosine nucleotide hydrolase; substrates include Ap3A, Ap2A, ADP‑ribose, NADH
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (119/133);  total (120/134)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 254 total aligned reads displayed.

TCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCAGG  >  CP009273/3519043‑3519332
                                                                                                                                                 |                                                                                                                                                
tcatATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAg                                                                                                                                                <  1:3535015/148‑1 (MQ=255)
tcatATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAg                                                                                                                                                <  2:2090819/148‑1 (MQ=255)
tcatATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAg                                                                                                                                                >  1:423483/1‑148 (MQ=255)
 catATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGa                                                                                                                                               <  1:642039/148‑1 (MQ=255)
   tatGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATc                                                                                                                                             >  2:1612033/1‑148 (MQ=255)
    atGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCg                                                                                                                                            >  1:2021273/1‑148 (MQ=255)
     tGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGt                                                                                                                                           >  2:38865/1‑148 (MQ=255)
      gcgcCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGtt                                                                                                                                          <  1:1875676/148‑1 (MQ=255)
        gccagcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCg                                                                                                                                        <  2:449536/148‑1 (MQ=255)
         ccagcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGc                                                                                                                                       <  1:2312635/148‑1 (MQ=255)
         ccagcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGc                                                                                                                                       <  1:3130998/148‑1 (MQ=255)
          cagcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCt                                                                                                                                      <  2:684701/148‑1 (MQ=255)
           agcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTc                                                                                                                                     <  1:1728416/148‑1 (MQ=255)
           agcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTc                                                                                                                                     <  1:900365/148‑1 (MQ=255)
            gcggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTcc                                                                                                                                    >  2:3642616/1‑148 (MQ=255)
               gccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCaaa                                                                                                                                 <  2:351797/148‑1 (MQ=255)
               gccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCaaa                                                                                                                                 <  2:27460/148‑1 (MQ=255)
               gccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCaaa                                                                                                                                 <  2:2209679/148‑1 (MQ=255)
                ccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAAt                                                                                                                                <  1:1508066/148‑1 (MQ=255)
                  agcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATcc                                                                                                                              >  1:2580316/1‑148 (MQ=255)
                       aCCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACct                                                                                                                         <  1:1677156/148‑1 (MQ=255)
                         cTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACctct                                                                                                                       <  1:3050169/148‑1 (MQ=255)
                            tgGTTGCGGCTCTGGCTGATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCAGGATATTCTTTTTGCTGGAAAAGTAATAGGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAATTCCAACCtcttct                                                                                                                    <  2:3787959/148‑1 (MQ=255)
                            tgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCtcttct                                                                                                                    >  1:2473155/1‑148 (MQ=255)
                                 gCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAg                                                                                                               <  2:2403959/148‑1 (MQ=255)
                                        tGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCgg                                                                                                        <  1:547177/148‑1 (MQ=255)
                                         ggCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGt                                                                                                       <  1:2325082/148‑1 (MQ=255)
                                           cTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTa                                                                                                     >  1:3379714/1‑148 (MQ=255)
                                             cATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGc                                                                                                   >  2:2077697/1‑148 (MQ=255)
                                                cGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGc                                                                                                >  1:2178004/1‑148 (MQ=255)
                                                cGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGc                                                                                                >  2:1610158/1‑148 (MQ=255)
                                                cGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGc                                                                                                >  1:1239932/1‑148 (MQ=255)
                                                   cTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTc                                                                                             >  1:1663144/1‑148 (MQ=255)
                                                       cAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAg                                                                                         <  1:681116/148‑1 (MQ=255)
                                                       cAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAg                                                                                         <  1:3688759/148‑1 (MQ=255)
                                                         gTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGAc                                                                                       >  1:2653141/1‑148 (MQ=255)
                                                          tGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACg                                                                                      <  2:1939755/148‑1 (MQ=255)
                                                          tGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACg                                                                                      >  2:571662/1‑148 (MQ=255)
                                                           gATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGc                                                                                     <  1:2714204/148‑1 (MQ=255)
                                                           gATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGc                                                                                     <  1:2045178/148‑1 (MQ=255)
                                                                cGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCa                                                                                >  2:2928728/1‑148 (MQ=255)
                                                                  ggTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTCGACGCTTTCAcc                                                                              >  1:1390905/1‑148 (MQ=255)
                                                                  ggTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCAcc                                                                              >  1:3294252/1‑148 (MQ=255)
                                                                     agagATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCgg                                                                           >  2:3097270/1‑148 (MQ=255)
                                                                     agagATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCgg                                                                           >  1:557930/1‑148 (MQ=255)
                                                                       agaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGAt                                                                         >  1:1010091/1‑148 (MQ=255)
                                                                       agaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGAt                                                                         >  1:3338231/1‑148 (MQ=255)
                                                                        gaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATc                                                                        >  1:3503409/1‑148 (MQ=255)
                                                                         aTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCa                                                                       <  2:2659784/148‑1 (MQ=255)
                                                                         aTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCa                                                                       <  2:1028506/148‑1 (MQ=255)
                                                                           cTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTCAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAAt                                                                     >  2:1422967/1‑148 (MQ=255)
                                                                           cTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTCAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAAt                                                                     >  2:3326214/1‑148 (MQ=255)
                                                                            ttGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAAtt                                                                    <  1:993370/148‑1 (MQ=255)
                                                                            ttGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAAtt                                                                    <  1:2904461/148‑1 (MQ=255)
                                                                             tGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTa                                                                   <  1:1794921/148‑1 (MQ=255)
                                                                              gcgcTACCACGATATTCATTTTGCTGTAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTaa                                                                  <  2:3414493/148‑1 (MQ=255)
                                                                              gcgcTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTaa                                                                  <  1:667166/148‑1 (MQ=255)
                                                                                gcTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATc                                                                <  2:2360390/148‑1 (MQ=255)
                                                                                 cTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATcc                                                               >  1:2426999/1‑148 (MQ=255)
                                                                                 cTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATcc                                                               <  1:3623751/148‑1 (MQ=255)
                                                                                 cTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATcc                                                               <  2:545914/148‑1 (MQ=255)
                                                                                 cTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATcc                                                               <  2:154511/148‑1 (MQ=255)
                                                                                  tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCt                                                              <  2:313798/148‑1 (MQ=255)
                                                                                  tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCt                                                              >  1:2973603/1‑148 (MQ=255)
                                                                                  tACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCt                                                              <  2:2988152/148‑1 (MQ=255)
                                                                                   aCCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCtt                                                             <  1:73365/148‑1 (MQ=255)
                                                                                    ccACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                            >  1:1450271/1‑148 (MQ=255)
                                                                                    ccACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                            <  2:884513/148‑1 (MQ=255)
                                                                                         atatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaaa                                                       >  1:749490/1‑148 (MQ=255)
                                                                                          tatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaaaa                                                      >  1:3270161/1‑148 (MQ=255)
                                                                                          tatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaaaa                                                      >  2:180665/1‑148 (MQ=255)
                                                                                           atTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaaaaa                                                     <  1:474018/148‑1 (MQ=255)
                                                                                            ttCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAc                                                    >  1:394831/1‑148 (MQ=255)
                                                                                             tCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAcc                                                   <  2:2963371/148‑1 (MQ=255)
                                                                                               aTTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTa                                                 <  2:1745650/148‑1 (MQ=255)
                                                                                                ttttGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTaa                                                <  2:2609657/148‑1 (MQ=255)
                                                                                                       ggAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAg                                         >  2:3504154/1‑148 (MQ=255)
                                                                                                        gAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAgg                                        <  2:578251/148‑1 (MQ=255)
                                                                                                         aaaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGa                                       <  2:935648/148‑1 (MQ=255)
                                                                                                            aGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTc                                    >  2:1917310/1‑148 (MQ=255)
                                                                                                             gTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCa                                   <  2:2272738/148‑1 (MQ=255)
                                                                                                              tAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATACGTAGGATTCAg                                  >  1:2927680/1‑148 (MQ=255)
                                                                                                                    ggTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCa                            <  2:2221843/148‑1 (MQ=255)
                                                                                                                    ggTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCa                            <  2:1323199/148‑1 (MQ=255)
                                                                                                                          aTGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGt                      <  2:1349704/148‑1 (MQ=255)
                                                                                                                           tGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTa                     <  2:2021273/148‑1 (MQ=255)
                                                                                                                           tGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTa                     >  2:2877949/1‑148 (MQ=255)
                                                                                                                            gCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTAt                    <  2:557860/148‑1 (MQ=255)
                                                                                                                             cTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTAtt                   <  1:1320147/148‑1 (MQ=255)
                                                                                                                              tGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTc                  >  1:874740/1‑148 (MQ=255)
                                                                                                                                aGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTcgc                >  2:2806714/1‑148 (MQ=255)
                                                                                                                                    tcttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggat            >  1:3748013/1‑148 (MQ=255)
                                                                                                                                     cttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatc           >  1:3210444/1‑148 (MQ=255)
                                                                                                                                     cttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatc           >  1:2808431/1‑148 (MQ=255)
                                                                                                                                      ttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatca          >  1:2529543/1‑148 (MQ=255)
                                                                                                                                      ttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatca          <  1:106905/148‑1 (MQ=255)
                                                                                                                                       tcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcag         >  2:428057/1‑148 (MQ=255)
                                                                                                                                         aaaaaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcagga       <  1:3590444/148‑1 (MQ=255)
                                                                                                                                          aaaaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcaggat      >  1:2227070/1‑148 (MQ=255)
                                                                                                                                            aaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcaggatca    >  1:105678/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCAGG  >  CP009273/3519043‑3519332

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: