Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,059,577 G→T 19.4% D473Y (GAT→TAT)  scpC → propionyl‑CoA:succinate CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,059,5770GT19.4% 703.0 / 113.7 247D473Y (GAT→TAT) scpCpropionyl‑CoA:succinate CoA transferase
Reads supporting (aligned to +/- strand):  ref base G (103/96);  new base T (23/25);  total (126/121)
Fisher's exact test for biased strand distribution p-value = 7.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.32e-01

Only 100 of 252 total aligned reads displayed.

TGATCATCACCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGC  >  CP009273/3059428‑3059726
                                                                                                                                                     |                                                                                                                                                     
   tcatcaCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATc                                                                                                                                                    <  1:932824/150‑1 (MQ=255)
   tcatcaCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATc                                                                                                                                                    >  2:3189785/1‑150 (MQ=255)
   tcatcaCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATc                                                                                                                                                    >  2:2863790/1‑150 (MQ=255)
   tcatcaCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATc                                                                                                                                                    >  2:189569/1‑150 (MQ=255)
    catcaCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCt                                                                                                                                                   >  2:4303321/1‑150 (MQ=255)
     atcaCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCtt                                                                                                                                                  <  1:2307212/150‑1 (MQ=255)
        aCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGc                                                                                                                                               >  1:4039883/1‑150 (MQ=255)
        aCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGc                                                                                                                                               >  1:4021796/1‑150 (MQ=255)
        aCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGc                                                                                                                                               >  2:1013898/1‑150 (MQ=255)
        aCCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGc                                                                                                                                               >  1:803828/1‑150 (MQ=255)
          cGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCa                                                                                                                                             <  1:992252/150‑1 (MQ=255)
            aaCAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACg                                                                                                                                           >  2:3794086/1‑150 (MQ=255)
              cAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTc                                                                                                                                         <  2:3096076/150‑1 (MQ=255)
               aaGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCt                                                                                                                                        >  1:3646002/1‑150 (MQ=255)
               aaGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCt                                                                                                                                        <  2:1341932/150‑1 (MQ=255)
                  ggATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCg                                                                                                                                     <  2:4111527/150‑1 (MQ=255)
                    aTCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGAc                                                                                                                                   >  2:1824639/1‑150 (MQ=255)
                     tCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGACAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACt                                                                                                                                  >  1:4324066/1‑150 (MQ=255)
                       gCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTa                                                                                                                                >  1:4432292/1‑150 (MQ=255)
                        ccGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTac                                                                                                                               <  1:3036498/150‑1 (MQ=255)
                          gATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTacac                                                                                                                             >  2:1356781/1‑150 (MQ=255)
                          gATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTacac                                                                                                                             >  1:3804521/1‑150 (MQ=255)
                           aTCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACAcc                                                                                                                            >  2:1168874/1‑150 (MQ=255)
                                  cgGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGtaattt                                                                                                                     >  1:1530737/1‑150 (MQ=255)
                                    gTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGtaatttaa                                                                                                                   <  2:1344312/150‑1 (MQ=255)
                                      cTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGtaatttaatt                                                                                                                 <  2:1412774/150‑1 (MQ=255)
                                      cTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGtaatttaatt                                                                                                                 <  2:2351503/150‑1 (MQ=255)
                                       tttCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTg                                                                                                                >  1:1866350/1‑150 (MQ=255)
                                       tttCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTg                                                                                                                >  1:3971645/1‑150 (MQ=255)
                                       tttCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTg                                                                                                                >  2:440298/1‑150 (MQ=255)
                                        ttCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATCTTCACCACTAGCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGc                                                                                                               >  2:2002927/1‑150 (MQ=255)
                                         tCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCa                                                                                                              <  2:4510126/150‑1 (MQ=255)
                                           cccGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAAc                                                                                                            >  2:4284434/1‑150 (MQ=255)
                                            ccGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAAcc                                                                                                           <  2:91546/150‑1 (MQ=255)
                                             cGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCg                                                                                                          <  1:3199508/150‑1 (MQ=255)
                                                 tCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTc                                                                                                      <  2:2028821/150‑1 (MQ=255)
                                                   aaCGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGa                                                                                                    >  2:243230/1‑150 (MQ=255)
                                                       cgcCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCt                                                                                                >  1:2657641/1‑150 (MQ=255)
                                                          ccGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTggg                                                                                             >  2:1491341/1‑150 (MQ=255)
                                                           cGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGt                                                                                            >  2:1069443/1‑150 (MQ=255)
                                                            gCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGtt                                                                                           <  2:47468/150‑1 (MQ=255)
                                                              aCTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTaa                                                                                         >  1:2903344/1‑150 (MQ=255)
                                                              aCTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTaa                                                                                         <  2:2780357/150‑1 (MQ=255)
                                                                tATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATg                                                                                       <  2:420522/150‑1 (MQ=255)
                                                                tATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATg                                                                                       <  2:414467/150‑1 (MQ=255)
                                                                    aTTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTc                                                                                   >  2:1774051/1‑150 (MQ=255)
                                                                     ttGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCg                                                                                  <  2:1939357/150‑1 (MQ=255)
                                                                        aTAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTcc                                                                               >  1:1032156/1‑150 (MQ=255)
                                                                        aTAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTACTTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTcc                                                                               >  2:1217435/1‑150 (MQ=255)
                                                                         tAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTccc                                                                              <  2:2240929/150‑1 (MQ=255)
                                                                           aTTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGt                                                                            >  2:1507305/1‑150 (MQ=255)
                                                                           aTTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGAAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGt                                                                            >  1:3087635/1‑150 (MQ=255)
                                                                            ttGTGCACATCCTCTGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTc                                                                           <  2:2170084/150‑1 (MQ=255)
                                                                            ttGTGCACATCCTCTGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTc                                                                           <  2:2170096/150‑1 (MQ=255)
                                                                              gtgCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCAc                                                                         <  2:2015460/150‑1 (MQ=255)
                                                                                   caTCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCt                                                                    >  1:1666925/1‑150 (MQ=255)
                                                                                      ccTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCa                                                                 >  2:1387445/1‑150 (MQ=255)
                                                                                       cTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGACCCGTCACGATCTTCAg                                                                >  2:3970498/1‑150 (MQ=255)
                                                                                              tCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATg                                                         <  1:3361974/150‑1 (MQ=255)
                                                                                               cGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGc                                                        >  2:3523214/1‑150 (MQ=255)
                                                                                                 ggATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTg                                                      <  2:1147032/150‑1 (MQ=255)
                                                                                                   aTTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGaa                                                    <  2:344641/150‑1 (MQ=255)
                                                                                                     tATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGc                                                  >  2:2697033/1‑150 (MQ=255)
                                                                                                      aTCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCa                                                 >  1:877717/1‑150 (MQ=255)
                                                                                                      aTCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCa                                                 <  2:2554058/150‑1 (MQ=255)
                                                                                                      aTCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCa                                                 <  1:2532437/150‑1 (MQ=255)
                                                                                                            aTCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTgg                                           >  2:3073549/1‑150 (MQ=255)
                                                                                                            aTCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCGTAGCCACGACTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGACCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTgg                                           >  2:738346/1‑150 (MQ=255)
                                                                                                                 tATCTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGta                                      >  1:4259402/1‑150 (MQ=255)
                                                                                                                 tATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGTTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGta                                      >  2:677936/1‑150 (MQ=255)
                                                                                                                    cTGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATaa                                   >  1:2426541/1‑150 (MQ=255)
                                                                                                                    cTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATaa                                   <  2:4611494/150‑1 (MQ=255)
                                                                                                                    cTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATaa                                   <  1:438519/150‑1 (MQ=255)
                                                                                                                     tGGAAAATGCGCCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAAc                                  <  1:3010318/150‑1 (MQ=255)
                                                                                                                        aaaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgat                               >  1:3390013/1‑150 (MQ=255)
                                                                                                                         aaaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatg                              >  1:724146/1‑150 (MQ=255)
                                                                                                                          aaTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatga                             >  1:427515/1‑150 (MQ=255)
                                                                                                                           aTGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACgatgat                            <  1:1136070/150‑1 (MQ=255)
                                                                                                                              cgcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTc                         <  1:3758798/150‑1 (MQ=255)
                                                                                                                               gcCTGGCGGCCATATTCACCACGCTCTTAGCCACGTCTTCTACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGTTTAATGACTCGTCCAGTCTCGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                        <  2:2376401/150‑1 (MQ=255)
                                                                                                                               gcCTGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                        <  1:1818229/150‑1 (MQ=255)
                                                                                                                               gcCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCt                        >  1:2221257/1‑150 (MQ=255)
                                                                                                                                 cTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGa                      >  2:3380381/1‑150 (MQ=255)
                                                                                                                                  tGGCGGACATATTCACCACTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                     >  1:1756144/1‑150 (MQ=255)
                                                                                                                                  tGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGAt                     <  1:1272252/150‑1 (MQ=255)
                                                                                                                                    gCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTg                   <  2:2000692/150‑1 (MQ=255)
                                                                                                                                      ggACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGaa                 >  2:2745288/1‑150 (MQ=255)
                                                                                                                                       gACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAAt                <  1:2200697/150‑1 (MQ=255)
                                                                                                                                        aCATATTCACCACGTTCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATa               <  2:1408497/150‑1 (MQ=255)
                                                                                                                                        aCATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATa               >  1:2113997/1‑150 (MQ=255)
                                                                                                                                         cATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaa              >  1:676011/1‑150 (MQ=255)
                                                                                                                                          atatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa             <  2:1005588/150‑1 (MQ=255)
                                                                                                                                          atatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa             >  1:1290502/1‑150 (MQ=255)
                                                                                                                                          atatTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATaaa             <  1:4334887/150‑1 (MQ=255)
                                                                                                                                             ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGat          >  1:3013117/1‑150 (MQ=255)
                                                                                                                                             ttCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGat          <  1:2000918/150‑1 (MQ=255)
                                                                                                                                                accacTATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATAtt       <  2:391089/150‑1 (MQ=255)
                                                                                                                                                  cacGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTcc     <  1:1884081/150‑1 (MQ=255)
                                                                                                                                                   acGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCt    <  1:2446652/150‑1 (MQ=255)
                                                                                                                                                     gATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGc  >  1:4599976/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGATCATCACCGAACAAGGGATCGCCGATCTGCGCGGTCTTTCCCCGCTTCAACGCGCCCGCACTATCATTGATAATTGTGCACATCCTATGTATCGGGATTATCTGCATCGCTATCTGGAAAATGCGCCTGGCGGACATATTCACCACGATCTTAGCCACGTCTTCGACTTACACCGTAATTTAATTGCAACCGGCTCGATGCTGGGTTAATGACTCGTCCCGTCACGATCTTCAGCAATATGCTGAAGCATCGTGGCGGTATAACGATGATTTCTGATTGAATAAAGATATTCCTGC  >  CP009273/3059428‑3059726

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: ‑