Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 37,692 G→T 9.6% Q45K (CAG→AAG)  caiC ← putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927337,6920GT 9.6% 761.0 / 13.1 241Q45K (CAG→AAG) caiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base G (117/89);  new base T (2/33);  total (119/122)
Fisher's exact test for biased strand distribution p-value = 4.19e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 242 total aligned reads displayed.

GGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCA  >  CP009273/37543‑37838
                                                                                                                                                     |                                                                                                                                                  
  cgTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGAt                                                                                                                                                  <  2:4609495/150‑1 (MQ=255)
  cgTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGAt                                                                                                                                                  >  1:1818636/1‑150 (MQ=255)
   gTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGAtt                                                                                                                                                 <  1:1225233/150‑1 (MQ=255)
    ttAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGAttt                                                                                                                                                >  1:911696/1‑150 (MQ=255)
    ttAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGAttt                                                                                                                                                <  1:2228672/150‑1 (MQ=255)
     tAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTa                                                                                                                                               <  2:2186117/150‑1 (MQ=255)
      aaTCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTaa                                                                                                                                              <  2:1334532/150‑1 (MQ=255)
       aTCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAAGTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAAc                                                                                                                                             <  2:4070000/150‑1 (MQ=255)
        tCGGCACCATACTCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATACATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACt                                                                                                                                            <  2:1326374/150‑1 (MQ=255)
          ggCACCATAATCGCGCCAATTTGTGGCAGCCCGAACCAGCAAAAGATTAAGTCCGGGCAGTTGTCGAGATGTAGTGCTAGCTTGTCGCCTTTGCGAATGCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTTATCTCCTGATTTAACTCa                                                                                                                                          <  1:3656390/150‑1 (MQ=255)
           gCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCaa                                                                                                                                         >  1:4317608/1‑150 (MQ=255)
            cACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAg                                                                                                                                        <  2:120310/150‑1 (MQ=255)
             aCCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGa                                                                                                                                       >  1:3054956/1‑150 (MQ=255)
             aCCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGa                                                                                                                                       <  1:4845109/150‑1 (MQ=255)
                aTAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATaa                                                                                                                                    <  1:2790921/150‑1 (MQ=255)
                 tAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAAc                                                                                                                                   >  1:162142/1‑150 (MQ=255)
                 tAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAAc                                                                                                                                   >  2:2565618/1‑150 (MQ=255)
                   aTCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACta                                                                                                                                 <  2:845416/150‑1 (MQ=255)
                      gcgcCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATAc                                                                                                                              <  2:1291252/150‑1 (MQ=255)
                        gcCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCg                                                                                                                            <  2:2879114/150‑1 (MQ=255)
                           aaTTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGtt                                                                                                                         <  1:1061228/150‑1 (MQ=255)
                           aaTTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGtt                                                                                                                         >  1:2550989/1‑150 (MQ=255)
                            aTTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTa                                                                                                                        <  1:248884/150‑1 (MQ=255)
                               tttGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACg                                                                                                                     >  2:4702262/1‑150 (MQ=255)
                               tttGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACg                                                                                                                     <  1:3246336/150‑1 (MQ=255)
                                 tGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGAc                                                                                                                   <  1:2583472/150‑1 (MQ=255)
                                    cAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTcc                                                                                                                >  1:4122763/1‑150 (MQ=255)
                                        ccGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAATATAACTATACCGGTTAACGACTccgccg                                                                                                            <  2:3226192/150‑1 (MQ=255)
                                        ccGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTccgccg                                                                                                            >  1:4389384/1‑150 (MQ=255)
                                         cGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCcgccgc                                                                                                           >  1:547677/1‑150 (MQ=255)
                                          gAACCAGCAAAAGATAATTTCCGGGCTGTTGTGGAGTTGTAGTGCAACGTGGTCGCGTTTGCGAAGCCCCCGCGTATAAAACCGGTTTGCCGTGGGTTTAATCTCCTTATTTAACTCAAGATAACTATACGGGTTAACGACTCCGCCGCt                                                                                                          <  2:2444782/150‑1 (MQ=255)
                                          gAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAATATAACTATACCGTTTAACGACTCCGCCGCt                                                                                                          <  2:1090064/150‑1 (MQ=255)
                                             ccAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGa                                                                                                       >  2:1977884/1‑150 (MQ=255)
                                             ccAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGa                                                                                                       >  1:2928045/1‑150 (MQ=255)
                                                  aaaaGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTcac                                                                                                  >  2:3974887/1‑150 (MQ=255)
                                                  aaaaGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTcac                                                                                                  >  1:4231158/1‑150 (MQ=255)
                                                  aaaaGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTcac                                                                                                  >  2:4338515/1‑150 (MQ=255)
                                                  aaaaGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTcac                                                                                                  >  2:1137140/1‑150 (MQ=255)
                                                    aaGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACaa                                                                                                >  1:1458635/1‑150 (MQ=255)
                                                    aaGATAAAGTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACaa                                                                                                <  2:4837858/150‑1 (MQ=255)
                                                        tAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCTAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTTATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCa                                                                                            <  2:4009776/150‑1 (MQ=255)
                                                        tAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAATATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCa                                                                                            <  2:2010072/150‑1 (MQ=255)
                                                         aaaTTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAg                                                                                           >  1:1489984/1‑150 (MQ=255)
                                                         aaaTGCCGGGCAGTTGTCGAGTTGTAGTGCACCCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACTGGTTTGCCGTGCGTTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAg                                                                                           <  2:3403085/150‑1 (MQ=255)
                                                          aaTTCCGGGCAGTTGTCGAGATGTAGTGCAACCTGGTCGCCTTTGCGAATCCCCAGCTTATAAAACAGGTTGGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAgc                                                                                          <  2:4630933/150‑1 (MQ=255)
                                                           aTTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAgcg                                                                                         >  1:2366917/1‑150 (MQ=255)
                                                                gggCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGtt                                                                                    >  2:2898109/1‑150 (MQ=255)
                                                                 ggCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTGGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGttt                                                                                   <  2:4545393/150‑1 (MQ=255)
                                                                    aGTTGTCGAGATTTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTGGCCGTGCGGTTAATCTCCTGATTTAACTCAATATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTAt                                                                                <  2:1632896/150‑1 (MQ=255)
                                                                      ttGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGa                                                                              >  2:2621718/1‑150 (MQ=255)
                                                                       tGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGAc                                                                             >  1:387574/1‑150 (MQ=255)
                                                                       tGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGAGTCACAAATCAGCGCCGTTTTATGAc                                                                             >  1:2967556/1‑150 (MQ=255)
                                                                        gTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGAcc                                                                            >  2:2465338/1‑150 (MQ=255)
                                                                           gagaTGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCCGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTa                                                                         <  2:2838147/150‑1 (MQ=255)
                                                                            agaTGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTaa                                                                        >  2:4514562/1‑150 (MQ=255)
                                                                             gaTGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAAGCTCCTGATTTAACTCAATATAAATATCCCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTaaa                                                                       <  1:3103528/150‑1 (MQ=255)
                                                                               tGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCCGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACg                                                                     <  1:3181136/150‑1 (MQ=255)
                                                                               tGTAGTGCAACCTTGTCGCCTTTGCGAAGCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACg                                                                     <  2:216823/150‑1 (MQ=255)
                                                                               tGTAGTGCAACCTGGTCGCCTTTGCGAATCCCCCGCGTATAAAACCGGTTGGCGGTGCGGTTAATCGCCTGATTTAACTCAAGATAACTATACCGGTTCACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACg                                                                     <  1:524709/150‑1 (MQ=255)
                                                                                 tAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTc                                                                   >  2:4583376/1‑150 (MQ=255)
                                                                                 tAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTc                                                                   >  2:3563735/1‑150 (MQ=255)
                                                                                 tAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTc                                                                   >  2:3861559/1‑150 (MQ=255)
                                                                                 tAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGATTGCCGTGCGGTAAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAAAGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGAACCTAAAAGTc                                                                   >  1:1352980/1‑150 (MQ=255)
                                                                                  aGTGAAACCTTGTCGCCTTGGCTAAACCCCAGCGTATAAAACTGGTTTGCCGTGCGTTTAATCTCCTTTTTTAACTCAAGATAACTATCCCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTcc                                                                  <  2:1795814/150‑1 (MQ=255)
                                                                                   gTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCg                                                                 >  2:3585675/1‑150 (MQ=255)
                                                                                         ccTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGAt                                                           >  1:3930759/1‑150 (MQ=255)
                                                                                           ttGTCGCATTTGCTACTCCCCCGCGTATAAAACCGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGAAAACTATACCGGTTAACGACTCCGCCGCGGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCg                                                         <  2:937382/150‑1 (MQ=255)
                                                                                              tCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTcc                                                      >  2:302858/1‑150 (MQ=255)
                                                                                              tCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTcc                                                      >  2:1779772/1‑150 (MQ=255)
                                                                                              tCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTcc                                                      >  2:2734161/1‑150 (MQ=255)
                                                                                               cGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTccc                                                     >  1:1509134/1‑150 (MQ=255)
                                                                                                 ccTTTGCGACCCCCCCGCGTATAAAACAGGTTTGCCGTGCGTTAAATCACCTTATTTAACAAAATATAACTAAACCGGTAAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCcac                                                   <  2:2029570/150‑1 (MQ=255)
                                                                                                  cTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCcaca                                                  >  1:113436/1‑150 (MQ=255)
                                                                                                     tGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTTATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACAttt                                               >  2:373462/1‑150 (MQ=255)
                                                                                                       cGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTTATTTAACTCAATATAACTATACCGTTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGa                                             <  2:1863525/150‑1 (MQ=255)
                                                                                                        gAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGAc                                            >  1:871256/1‑150 (MQ=255)
                                                                                                         aaTCCCCAGCGTATAAAACTGGTTGGCCGTGCGTTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACg                                           <  2:1465595/150‑1 (MQ=255)
                                                                                                              cctgCGTATAAAACCGTTTTGCGGTGCGTTTAATCTCGTTTTTAAACTAAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGAt                                      <  2:4123789/147‑1 (MQ=255)
                                                                                                               cAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGAtg                                     >  1:1989853/1‑150 (MQ=255)
                                                                                                               cAGCGTAAAAAACTGTTTGGCGGTGCGTTTAATCTCATTATTTAACTAAATAAAACTAAACGGGTCAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGAtg                                     <  1:1809477/150‑1 (MQ=255)
                                                                                                                aGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTGTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGAtgt                                    >  1:2245953/1‑150 (MQ=255)
                                                                                                                     ataAAACCGGTTTGCCGTGCGGTTAATCTCCTGTTTTAACTCAAGATAACTATCCCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTcc                               <  2:2523974/150‑1 (MQ=255)
                                                                                                                     ataAAACAGGTTTGACGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATTAAAGTAAAAGTAAGCAAGATAGTACCACAATTGACGAAGAtgttgtga                               >  2:1608867/1‑148 (MQ=255)
                                                                                                                      taAAACAGGTTTGCCGTGCGTTTAATCTCCTGATTTAACACAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCg                              <  2:4288436/150‑1 (MQ=255)
                                                                                                                            aGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATg                        >  2:1764241/1‑150 (MQ=255)
                                                                                                                            aGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATg                        >  2:4220508/1‑150 (MQ=255)
                                                                                                                             ggTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCCATGa                       >  1:1303787/1‑150 (MQ=255)
                                                                                                                                ttGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATTACCGTAAACGTCCGCAAGATCGTCCCACAATTTACGTAGATGTTGTCCGCCAATGatat                    >  1:4070680/1‑150 (MQ=255)
                                                                                                                                ttGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGatat                    >  2:1917196/1‑150 (MQ=255)
                                                                                                                                 tGCGGTGGGTTTAATCTCCTTATTTAACTAATGATAATTATACGGGTAAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATc                   <  1:3379843/150‑1 (MQ=255)
                                                                                                                                 tGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTTTCCGCCAATGATATc                   >  1:4793440/1‑150 (MQ=255)
                                                                                                                                 tGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATc                   >  1:3844818/1‑150 (MQ=255)
                                                                                                                                     gTGGGGTTAATATCCTTATTTAACAAAATATAACTACACGGGTAAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCAtt               <  1:1708362/150‑1 (MQ=255)
                                                                                                                                      tGCGGTTAATCTCCTTATTTAACTCAAGATAACTATCCCGTTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTg              <  2:2356483/150‑1 (MQ=255)
                                                                                                                                        cgtttAATCTCCTTATTTAAATCAATATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCa            <  2:911696/147‑1 (MQ=255)
                                                                                                                                         ggTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCAc           >  2:1420558/1‑150 (MQ=255)
                                                                                                                                         ggTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCAc           >  2:3756767/1‑150 (MQ=255)
                                                                                                                                             aaTCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACctct       >  2:653058/1‑150 (MQ=255)
                                                                                                                                                ctcCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATc    <  2:1946938/150‑1 (MQ=255)
                                                                                                                                                  ccTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATTACCGTAAAAGTCCGCCAGATCGTCCCACAATTTACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCa  >  1:473242/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
GGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGATAAATTCCGGGCAGTTGTCGAGATGTAGTGCAACCTTGTCGCCTTTGCGAATCCCCAGCGTATAAAACAGGTTTGCCGTGCGGTTAATCTCCTGATTTAACTCAAGATAACTATACCGGTTAACGACTCCGCCGCTGGATTCACAAATCAGCGCCGTTTTATGACCGTAAACGTCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCA  >  CP009273/37543‑37838

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: