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,1880CA95.6% 655.0 / 20.8 227L107L (CTG→CTTnudEadenosine nucleotide hydrolase; substrates include Ap3A, Ap2A, ADP‑ribose, NADH
Reads supporting (aligned to +/- strand):  ref base C (7/3);  new base A (95/122);  total (102/125)
Fisher's exact test for biased strand distribution p-value = 1.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.08e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 227 total aligned reads displayed.

TCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCA  >  CP009273/3519043‑3519330
                                                                                                                                                 |                                                                                                                                              
 catATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGa                                                                                                                                             >  2:273222/1‑148 (MQ=255)
  atatGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGAt                                                                                                                                            <  2:1482654/148‑1 (MQ=255)
   tatGCGACAGCGTCCAGCGCACCTGTGGTAGAGGCTCTGGCTCATCGCCTACCAGTGATTCCGGGTAGATATCTAGCGCTACCACGATATTAATTTTTCTGGAAAAGTAAGACGGAGCCATGCTGAGCTTCTTCAAAAAAGTAAGATc                                                                                                                                           >  2:3175134/1‑148 (MQ=255)
    atGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCg                                                                                                                                          >  2:2162965/1‑148 (MQ=255)
    atGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCg                                                                                                                                          >  2:2036675/1‑148 (MQ=255)
    atGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCg                                                                                                                                          <  1:1820265/148‑1 (MQ=255)
    atGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCg                                                                                                                                          <  1:632226/148‑1 (MQ=255)
     tGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGt                                                                                                                                         <  1:3538386/148‑1 (MQ=255)
             cggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCa                                                                                                                                 >  1:2415806/1‑148 (MQ=255)
              ggccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCaa                                                                                                                                <  1:1481304/148‑1 (MQ=255)
               gccagcgCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCaaa                                                                                                                               >  1:868888/1‑148 (MQ=255)
                   gcgcACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCa                                                                                                                           <  1:2327200/148‑1 (MQ=255)
                      cACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAAcc                                                                                                                        <  2:191134/148‑1 (MQ=255)
                         cTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACctct                                                                                                                     <  1:1227763/148‑1 (MQ=255)
                          tgtgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCtctt                                                                                                                    >  2:3369036/1‑148 (MQ=255)
                           gtgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCtcttc                                                                                                                   <  2:1654402/148‑1 (MQ=255)
                            tgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCtcttct                                                                                                                  <  1:991577/148‑1 (MQ=255)
                            tgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCtcttct                                                                                                                  <  1:1822615/148‑1 (MQ=255)
                            tgGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCtcttct                                                                                                                  <  1:17188/148‑1 (MQ=255)
                             ggTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTcttctt                                                                                                                 <  1:1255939/148‑1 (MQ=255)
                               tAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCtttt                                                                                                               <  1:3125982/148‑1 (MQ=255)
                                 gCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAg                                                                                                             >  2:3001732/1‑148 (MQ=255)
                                   ggCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCt                                                                                                           <  1:1370612/148‑1 (MQ=255)
                                    gCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTc                                                                                                          >  1:1581545/1‑148 (MQ=255)
                                     ctctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcg                                                                                                         >  1:2518756/1‑148 (MQ=255)
                                      tctGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTcgc                                                                                                        <  2:2838431/148‑1 (MQ=255)
                                         ggCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGt                                                                                                     <  1:366745/148‑1 (MQ=255)
                                           cTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTa                                                                                                   >  2:1847495/1‑148 (MQ=255)
                                            tCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAg                                                                                                  <  1:2273621/148‑1 (MQ=255)
                                             cATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTAACACGATATTCATTTTCCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTCCTTGTAGATCGCGCTTAGc                                                                                                 >  1:2782492/1‑148 (MQ=255)
                                                cGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGc                                                                                              <  1:2832900/148‑1 (MQ=255)
                                                 gCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCt                                                                                             >  2:1288607/1‑148 (MQ=255)
                                                   cTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTc                                                                                           >  2:2569667/1‑148 (MQ=255)
                                                        aGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCCACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGa                                                                                      >  1:884346/1‑148 (MQ=255)
                                                          tGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACg                                                                                    <  2:1010304/148‑1 (MQ=255)
                                                          tGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACg                                                                                    <  2:768067/148‑1 (MQ=255)
                                                          tGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACg                                                                                    <  1:3492988/148‑1 (MQ=255)
                                                          tGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACg                                                                                    <  2:2018469/148‑1 (MQ=255)
                                                            aTTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCCACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCt                                                                                  >  1:1031374/1‑148 (MQ=255)
                                                            aTTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCt                                                                                  >  1:605045/1‑148 (MQ=255)
                                                              tCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCttt                                                                                <  2:2173987/148‑1 (MQ=255)
                                                                 gggTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCAc                                                                             <  2:2081509/148‑1 (MQ=255)
                                                                 gggTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCAc                                                                             <  2:238428/148‑1 (MQ=255)
                                                                      gagaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCCACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGa                                                                        >  1:2707418/1‑148 (MQ=255)
                                                                      gagaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGa                                                                        <  1:3025519/148‑1 (MQ=255)
                                                                      gagaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGa                                                                        >  2:2975342/1‑148 (MQ=255)
                                                                       agaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGAt                                                                       >  2:1549705/1‑148 (MQ=255)
                                                                       agaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGAt                                                                       >  2:1549583/1‑148 (MQ=255)
                                                                        gaTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATc                                                                      >  1:2182742/1‑148 (MQ=255)
                                                                         aTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCa                                                                     >  1:686132/1‑148 (MQ=255)
                                                                         aTCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCa                                                                     >  1:1308126/1‑148 (MQ=255)
                                                                          tCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCaa                                                                    <  2:1264836/148‑1 (MQ=255)
                                                                          tCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCaa                                                                    <  1:2116496/148‑1 (MQ=255)
                                                                           cTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAAt                                                                   >  1:2025217/1‑148 (MQ=255)
                                                                             tGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTa                                                                 <  1:2752597/148‑1 (MQ=255)
                                                                              gcgcTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTaa                                                                <  2:2797974/148‑1 (MQ=255)
                                                                               cgcTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAAt                                                               <  2:1713772/148‑1 (MQ=255)
                                                                                   aCCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTATTCTTTTAGCTCGCGGTTATCGGCTTCGTAGACGCTGTCACCCGGATCAATTAATCCtt                                                           >  2:2399213/1‑148 (MQ=255)
                                                                                    ccACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                          >  1:2506961/1‑148 (MQ=255)
                                                                                    ccACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                          >  1:1683698/1‑148 (MQ=255)
                                                                                    ccACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTATCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCttt                                                          >  1:2046041/1‑148 (MQ=255)
                                                                                      aCGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCg                                                        >  2:885334/1‑148 (MQ=255)
                                                                                          tatTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGaaaa                                                    <  1:486761/148‑1 (MQ=255)
                                                                                            ttCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAc                                                  <  1:982205/148‑1 (MQ=255)
                                                                                             tCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAAcc                                                 <  2:1113037/148‑1 (MQ=255)
                                                                                               aTTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTa                                               <  2:1427780/148‑1 (MQ=255)
                                                                                               aTTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTa                                               <  2:3447767/148‑1 (MQ=255)
                                                                                               aTTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTa                                               <  1:2649460/148‑1 (MQ=255)
                                                                                                   tGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTc                                           >  1:2771106/1‑148 (MQ=255)
                                                                                                     cTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGt                                         <  1:3558678/148‑1 (MQ=255)
                                                                                                      tGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTa                                        <  2:2518756/148‑1 (MQ=255)
                                                                                                      tGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTa                                        <  1:3369036/148‑1 (MQ=255)
                                                                                                       ggAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAg                                       >  1:3549433/1‑148 (MQ=255)
                                                                                                       ggAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAg                                       >  1:1968552/1‑148 (MQ=255)
                                                                                                        gAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAgg                                      >  1:2447848/1‑148 (MQ=255)
                                                                                                         aaaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGa                                     <  1:3297821/148‑1 (MQ=255)
                                                                                                         aaaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGa                                     <  1:2597722/148‑1 (MQ=255)
                                                                                                          aaaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGAt                                    <  1:2481679/148‑1 (MQ=255)
                                                                                                           aaGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGAtt                                   <  2:2267495/148‑1 (MQ=255)
                                                                                                              tAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAg                                >  2:709391/1‑148 (MQ=255)
                                                                                                               aaGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGt                               >  2:766257/1‑148 (MQ=255)
                                                                                                                aGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGtt                              <  1:1087121/148‑1 (MQ=255)
                                                                                                                 gACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTc                             >  2:1619343/1‑148 (MQ=255)
                                                                                                                 gACGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTc                             >  2:1477439/1‑148 (MQ=255)
                                                                                                                   cGGTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTccc                           <  2:955418/148‑1 (MQ=255)
                                                                                                                     gTGCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCAc                         >  1:264285/1‑148 (MQ=255)
                                                                                                                       gCCATGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTg                       <  1:296650/148‑1 (MQ=255)
                                                                                                                           tGCTGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTa                   <  2:2484799/148‑1 (MQ=255)
                                                                                                                              tGAGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTc                <  1:1129136/148‑1 (MQ=255)
                                                                                                                                aGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTcgc              >  1:1188030/1‑148 (MQ=255)
                                                                                                                                aGCTTCTTCAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGAATAcgc              >  1:3449758/1‑148 (MQ=255)
                                                                                                                                    tcttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAACCCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggat          <  2:3523795/148‑1 (MQ=255)
                                                                                                                                    tcttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggat          <  2:450766/148‑1 (MQ=255)
                                                                                                                                    tcttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggat          >  2:814538/1‑148 (MQ=255)
                                                                                                                                      ttcAAAAAAGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatca        <  2:620422/148‑1 (MQ=255)
                                                                                                                                         aaaaaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcagga     <  2:1031533/148‑1 (MQ=255)
                                                                                                                                          aaaaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcaggat    <  1:2327743/148‑1 (MQ=255)
                                                                                                                                            aaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTACTTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcaggatca  <  2:1808929/148‑1 (MQ=255)
                                                                                                                                            aaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcaggatca  >  1:2290224/1‑148 (MQ=255)
                                                                                                                                            aaaGTAAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCggatcaggatca  <  1:2569667/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
TCATATGCGCCAGCGGCCAGCGCACCTGTGGTAGCGGCTCTGGCTCATCGCCTTCCAGTGATTCCGGGTAGAGATCTTGCGCTACCACGATATTCATTTTGCTGGAAAAGTAAGACGGTGCCATGCTGAGCTTCTTCAAAAAAGTCAGATCGTTCGCTCCAAATCCAACCTCTTCTTTTAGCTCGCGGTTAGCGGCTTCGTAGACGCTTTCACCCGGATCAATTAATCCTTTCGAAAAACCTAATTCGTAGGATTCAGTTCCCACTGCGTATTCGCGGATCAGGATCA  >  CP009273/3519043‑3519330

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: