Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,601,871 C→A 31.5% P90Q (CCG→CAG)  tam → trans‑aconitate methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,601,8710CA31.5% 577.2 / 247.9 323P90Q (CCG→CAG) tamtrans‑aconitate methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (89/132);  new base A (45/57);  total (134/189)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 327 total aligned reads displayed.

TGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCAT  >  CP009273/1601722‑1602019
                                                                                                                                                     |                                                                                                                                                    
    ccAGGTAACAGCACCGCCCTTCTACAACACCGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTa                                                                                                                                                  <  1:1016868/150‑1 (MQ=255)
     cAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTAc                                                                                                                                                 <  2:4550082/150‑1 (MQ=255)
     cAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTAc                                                                                                                                                 >  2:4764311/1‑150 (MQ=255)
     cAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTAc                                                                                                                                                 <  2:5048301/150‑1 (MQ=255)
      aGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACa                                                                                                                                                >  2:4572405/1‑150 (MQ=255)
      aGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACa                                                                                                                                                >  2:1785460/1‑150 (MQ=255)
        gTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAgg                                                                                                                                              <  2:336929/150‑1 (MQ=255)
          aaCAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCa                                                                                                                                            >  1:1836208/1‑150 (MQ=255)
           aCAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCAc                                                                                                                                           <  2:1296871/150‑1 (MQ=255)
            cAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACt                                                                                                                                          >  2:4072334/1‑150 (MQ=255)
            cAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGAGAACAGGCATCGACTCGTCTCCGGCGATGAGTGCTGAAGCGCGCAGTGCTTTGCCAGACAGCCAGTTCGTGGAAGCAGACAGCCGCAACAGGCAACCGGTACAGGCACt                                                                                                                                          >  2:134948/1‑150 (MQ=255)
            cAGCACCGCCCTTCTACAAAAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGACAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACt                                                                                                                                          >  1:4839876/1‑150 (MQ=255)
             aGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTc                                                                                                                                         >  1:5314710/1‑150 (MQ=255)
             aGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTc                                                                                                                                         >  1:71872/1‑150 (MQ=255)
                aCCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGAt                                                                                                                                      >  2:4991641/1‑150 (MQ=255)
                   gCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTg                                                                                                                                   >  1:1230924/1‑150 (MQ=255)
                   gCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCAATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTg                                                                                                                                   >  1:1869724/1‑150 (MQ=255)
                    cccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGa                                                                                                                                  <  1:971226/150‑1 (MQ=255)
                     ccTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGAt                                                                                                                                 >  1:5115743/1‑150 (MQ=255)
                      cTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAtt                                                                                                                                <  2:1886392/150‑1 (MQ=255)
                       ttCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAttt                                                                                                                               <  1:2316329/150‑1 (MQ=255)
                       ttCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGAttt                                                                                                                               <  1:1490637/150‑1 (MQ=255)
                          tACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTg                                                                                                                            <  2:2741143/150‑1 (MQ=255)
                          tACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTg                                                                                                                            <  1:4340933/150‑1 (MQ=255)
                             aacaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTa                                                                                                                         >  2:4218980/1‑150 (MQ=255)
                             aacaacGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTa                                                                                                                         >  1:1902168/1‑150 (MQ=255)
                                      ggCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCAc                                                                                                                <  2:2058850/150‑1 (MQ=255)
                                      ggCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCAc                                                                                                                <  2:3710713/150‑1 (MQ=255)
                                       gCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACt                                                                                                               <  2:3549179/150‑1 (MQ=255)
                                          tGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCa                                                                                                            >  1:2153122/1‑150 (MQ=255)
                                            cGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAAt                                                                                                          <  2:248378/150‑1 (MQ=255)
                                            cGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAAt                                                                                                          >  2:4342997/1‑150 (MQ=255)
                                             ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATg                                                                                                         <  2:4943521/150‑1 (MQ=255)
                                             ggCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATg                                                                                                         <  2:21234/150‑1 (MQ=255)
                                              gCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATgg                                                                                                        <  2:1157341/150‑1 (MQ=255)
                                                cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                                                                                                      <  2:2602895/150‑1 (MQ=255)
                                                cAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCt                                                                                                      <  2:4534557/150‑1 (MQ=255)
                                                 aGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTg                                                                                                     >  1:312167/1‑150 (MQ=255)
                                                  ggATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGc                                                                                                    >  2:39188/1‑150 (MQ=255)
                                                   gATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGcc                                                                                                   <  1:4685824/150‑1 (MQ=255)
                                                    aTAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGccc                                                                                                  <  2:169577/150‑1 (MQ=255)
                                                    aTAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGccc                                                                                                  >  1:278344/1‑150 (MQ=255)
                                                    aTAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGccc                                                                                                  >  1:5291366/1‑150 (MQ=255)
                                                         aGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCa                                                                                             >  2:647359/1‑150 (MQ=255)
                                                         aGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCa                                                                                             <  1:3394788/150‑1 (MQ=255)
                                                          ggCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCAc                                                                                            <  2:729546/150‑1 (MQ=255)
                                                            cATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTa                                                                                          >  2:2756707/1‑150 (MQ=255)
                                                            cATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTa                                                                                          >  2:3024465/1‑150 (MQ=255)
                                                             aTCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTAc                                                                                         >  1:3818499/1‑150 (MQ=255)
                                                                gACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGaa                                                                                      >  1:2934190/1‑150 (MQ=255)
                                                                 aCTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAt                                                                                     <  1:5238864/150‑1 (MQ=255)
                                                                  cTCGTCTCCGGCTATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAtt                                                                                    >  1:1666875/1‑150 (MQ=255)
                                                                  cTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAAtt                                                                                    <  2:3639028/150‑1 (MQ=255)
                                                                   tCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTg                                                                                   <  2:1411776/150‑1 (MQ=255)
                                                                     gTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGtt                                                                                 <  1:2710819/150‑1 (MQ=255)
                                                                      tctcCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGttt                                                                                <  2:3868763/150‑1 (MQ=255)
                                                                         ccGGTGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCt                                                                             <  2:2822075/150‑1 (MQ=255)
                                                                         ccGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCt                                                                             <  1:3528242/150‑1 (MQ=255)
                                                                             cGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAATCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCTATGGCTGCCCGACCACTACGAATTGTTTCCTCATc                                                                         <  2:2213106/150‑1 (MQ=255)
                                                                             cGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATc                                                                         >  1:1812853/1‑150 (MQ=255)
                                                                                tgatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTgg                                                                      >  1:2390914/1‑150 (MQ=255)
                                                                                 gatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGt                                                                     >  1:4943526/1‑150 (MQ=255)
                                                                                 gatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGt                                                                     <  2:4470584/150‑1 (MQ=255)
                                                                                 gatTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGt                                                                     >  1:1974307/1‑150 (MQ=255)
                                                                                   ttGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGttt                                                                   <  1:3652921/150‑1 (MQ=255)
                                                                                    tGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTc                                                                  <  2:1947578/150‑1 (MQ=255)
                                                                                        gAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTa                                                              <  1:4880443/150‑1 (MQ=255)
                                                                                        gAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTa                                                              >  1:2495967/1‑150 (MQ=255)
                                                                                          aGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACt                                                            >  2:3111193/1‑150 (MQ=255)
                                                                                              cgcAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAAt                                                        <  2:388529/150‑1 (MQ=255)
                                                                                                cAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATcc                                                      <  2:236598/150‑1 (MQ=255)
                                                                                                cAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATcc                                                      <  1:2211470/150‑1 (MQ=255)
                                                                                                    gCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAg                                                  <  2:802143/150‑1 (MQ=255)
                                                                                                    gCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAg                                                  <  2:3244681/150‑1 (MQ=255)
                                                                                                      tttGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAggg                                                <  1:3560232/150‑1 (MQ=255)
                                                                                                           cAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGc                                           >  2:4423759/1‑150 (MQ=255)
                                                                                                            aGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCt                                          <  1:4072334/150‑1 (MQ=255)
                                                                                                                  ccAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGt                                    <  2:242115/150‑1 (MQ=255)
                                                                                                                  ccAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGt                                    <  1:4986938/150‑1 (MQ=255)
                                                                                                                   cAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTa                                   <  1:2627962/150‑1 (MQ=255)
                                                                                                                     gTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACa                                 <  2:4853430/150‑1 (MQ=255)
                                                                                                                        cGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGAt                              <  2:3796949/150‑1 (MQ=255)
                                                                                                                          tGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGc                            >  1:5184110/1‑150 (MQ=255)
                                                                                                                          tGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGc                            <  1:4609614/150‑1 (MQ=255)
                                                                                                                           ggAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGcc                           <  2:1055127/150‑1 (MQ=255)
                                                                                                                           ggAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGcc                           <  2:5288637/150‑1 (MQ=255)
                                                                                                                              aGCAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGa                        >  1:3224983/1‑150 (MQ=255)
                                                                                                                               gCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGAt                       >  2:112295/1‑150 (MQ=255)
                                                                                                                                cAGACATCCGCAACTGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATa                      <  1:4791101/150‑1 (MQ=255)
                                                                                                                                 aGACATCCGCAACTGGCAACCGGTACGGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATaa                     <  1:675762/150‑1 (MQ=255)
                                                                                                                                      tCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATCActggc                >  1:4143865/1‑150 (MQ=255)
                                                                                                                                      tCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCATGCTGGCAGTACAGATGCCAGATAActggc                <  2:2194652/150‑1 (MQ=255)
                                                                                                                                       ccGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAActggct               <  1:618548/150‑1 (MQ=255)
                                                                                                                                              tGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCg        <  1:1345158/150‑1 (MQ=255)
                                                                                                                                              tGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCg        <  1:1038605/150‑1 (MQ=255)
                                                                                                                                              tGGCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCg        <  1:4553157/150‑1 (MQ=255)
                                                                                                                                               ggCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGa       <  2:2495967/150‑1 (MQ=255)
                                                                                                                                               ggCAACAGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGa       <  1:3941244/150‑1 (MQ=255)
                                                                                                                                                  aaCCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGAccc    >  2:4632104/1‑150 (MQ=255)
                                                                                                                                                  aaCCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGAccc    <  2:1582019/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGCCCAGGTAACAGCACCGCCCTTCTACAACAACGTTGGCCTGCGGCCAGGATAACAGGCATCGACTCGTCTCCGGCGATGATTGCTGAAGCGCGCAGTGCTTTGCCAGACTGCCAGTTCGTGGAAGCAGACATCCGCAACTGGCAACCGGTACAGGCACTCGATCTGATTTTTGCTAATGCCTCACTGCAATGGCTGCCCGACCACTACGAATTGTTTCCTCATCTGGTTTCTTTACTTAATCCGCAGGGCGTGCTGGCAGTACAGATGCCAGATAACTGGCTGGAGCCGACCCAT  >  CP009273/1601722‑1602019

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: