Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,284 T→G 15.0% T181P (ACC→CCC)  caiC ← putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091337,2840TG15.0% 122.3 / 2.1 68T181P (ACC→CCC) caiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (16/37);  new base G (13/2);  total (29/39)
Fisher's exact test for biased strand distribution p-value = 1.77e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCC  >  NC_000913/37051‑37498
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                         
tACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTAAACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACCCGGCCAGGCGGAGTAATATCCAGCGCAGCGCAGGATGTCATGGGTAAT‑CACCCCACCATCCGGTCGGGCGGTGGTGCCGGGGGGGAagg                                                                                                                                                                                                                      >  2:499463/1‑238 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTCTTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCAACGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCCGGTAGACGTCATCCTCCCGCAGTGCACACTGCCAGACCGAGTAATATCCATCGCACCGCAGGTTGTAATGGGTAAT‑CACCAAACCTTTCTGGCGGGAGGTGGGGCGGGGGGGGAAggga                                                                                                                                                                                                                    >  2:481078/1‑239 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACCCCTTTCGGTCGGCAGGGGGGGCCGGGGGTGAAgaga                                                                                                                                                                                                                    >  1:63220/1‑239 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGGGGTGGTGCCGGAGGGGAAgaga                                                                                                                                                                                                                    >  1:212264/1‑239 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAgaga                                                                                                                                                                                                                    >  1:358393/1‑239 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGGAAT‑CACACCACCTTTCGGTCGGGAGGGGGGGCGGGGGGGGaaaaaa                                                                                                                                                                                                                    >  2:254437/1‑235 (MQ=255)
     gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGGGCCGGCGGTGAAGAGAAtt                                                                                                                                                                                                                 >  1:55759/1‑239 (MQ=255)
     gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACAGCCAGGCGGAGTAATATCCAGCGAAGCGCCGGGTGTAATGGGAAAT‑CACCACACCATTCGGTCGGGAGGTGGGGCCGGAGGGGAAGAAAAtt                                                                                                                                                                                                                 >  2:185626/1‑239 (MQ=255)
     gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGGAGACGTCCTCTTCACGCAGTGCACACTGCCAGACGGCGGAATATCCACCGAAGCGCCGGTTGTAATGGGTAAT‑AACCACCACCTTCGGTCGGGAGGGGGGGCCGGAGGGGAAAGGAAtt                                                                                                                                                                                                                 >  2:115046/1‑239 (MQ=255)
      tCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGAACACTGCCAGGCGGAGTAATATCCCGCGAGGCGCAGGTTGTAATGGGTAAT‑CACCAAACCTTTCGGTCGGGAGGTGGGGCCGGGGGGGGagaggaat                                                                                                                                                                                                                 >  2:826409/1‑234 (MQ=255)
             aaGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCCGGTTGTAATGGGGAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGt                                                                                                                                                                                                         >  1:532869/1‑239 (MQ=255)
                gcgcGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GCTATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGGCACGCAGTGCACACTGCCAGGCGGGGGAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAG‑AACCACACCTTTCGGTCGGGCGGTGGGGCCGGAGGGGGAAAGAATTTcacccggctc                                                                                                                                                                                                      >  2:96052/1‑230 (MQ=255)
                gcgcGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCGTATc                                                                                                                                                                                                      >  1:314014/1‑239 (MQ=255)
                 cgcgGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCCGGGTGGAATGGGTAAT‑CACCCCACCCTTAGGTCGGGAGGGGGGGCCGGGGGGGAAGAGAATATCCCCCGTATCg                                                                                                                                                                                                     >  2:795337/1‑239 (MQ=255)
                    gggCGCTGTATTTCGCGACCAGCACAAAGGTGGCCCCGGCAGAAAACCCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGGCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCa                                                                                                                                                                                                  <  1:51563/239‑1 (MQ=255)
                            tttttCCCGCCCAGCAAGAAGGGGCCCCCGGCAAAAAACCCCC‑CCACCGC‑GGCAGTACCCTGGCAACC‑GATATGAACCGCAGGCATTCCCATCAGGTAGCCGTCATCGTCCCGCAGGGCACACGGCCAGGCGGAGTAATATCCAGCGAAGCGCAGTTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAg                                                                                                                                                                                          <  2:575293/237‑1 (MQ=255)
                             aTTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGAAATC‑GATATGAACCGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAg                                                                                                                                                                                          <  1:790995/238‑1 (MQ=255)
                                 ctcGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTccatgggaaagggcg                                                                                                                                                                                     >  2:790995/1‑225 (MQ=255)
                                  ccacgccgacaaaaGGGGGCCCCGGGAAAAAACCCCC‑CCCCTGCGCGGAGGACCATGGCAAGCGGATA‑GAAACCGAGGAATTCCCGGAGGGAAGACGTCCTCGTCACGCAGTGCCCCCTGCAAGGCGGAGTAATACCCAGCGAAGCGGAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGGCCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAgcggcg                                                                                                                                                                                     <  2:180724/228‑1 (MQ=255)
                                      ccaccaCAAAGGTACCCCGGGAAAAAAACCCCG‑CCATCCC‑CGCAGTACCCTGGCAATC‑GAAATGAAACGAAGGCATTACCGTAAGGTAGCCGTCAGCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                                                <  2:817500/239‑1 (MQ=255)
                                      ccAGCAAAAGGGTGGCCCCGGAAAAAACCGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                                                <  1:866592/239‑1 (MQ=255)
                                      caagaacaGAGGTGGCCCCGGCAAAAAACGGCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                                                <  1:795337/234‑1 (MQ=255)
                                         gCAAAATGGTGGCCCCGGCAGAAACCGCCG‑CCATCGC‑CGCAGTACACTGGCAACC‑GATATGAACCGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAg                                                                                                                                                                             <  1:390756/239‑1 (MQ=255)
                                         cccaacaGGTCGCCCCGG‑AAAAAACCCCGCCCAGCGC‑GCCGGTCCAGGGGAAATC‑GATATGAAACGCAGGCATTCCGGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAg                                                                                                                                                                             <  2:202946/233‑1 (MQ=255)
                                          cacaAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGGCAGTCGATAGCGGCGGTGGATAGc                                                                                                                                                                            >  1:56969/1‑239 (MQ=255)
                                          cacaAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTGGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTAACGTCAGTCGATACcggcgggccgaag                                                                                                                                                                             >  1:747298/1‑231 (MQ=255)
                                            caAAGGGTGCCCCGGCAGAAAACGCCC‑CCATGCC‑CGCAGTACCCTGGCAATC‑GATATGAAACGCAGGCATTCCCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGcac                                                                                                                                                                          <  1:760655/239‑1 (MQ=255)
                                                         ggCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGAAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGGGCCGGGGGGGAAGAGAATTTCCGCCGTTACGGCCGTCGAAAGgggggggggggaacccaacaggggggggg                                                                                                                                                             >  2:802276/1‑216 (MQ=255)
                                                           cAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGGGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATTGCCAAAGgggggggggggt                                                                                                                                                           >  1:257119/1‑231 (MQ=255)
                                                                                          tGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATCGCGGCGGTGCATCGCCCAAGGAGGCGGGATGTTGATTTTTCCCGTGA‑GGAAAcacaccaa                                                                                                                              >  1:910764/1‑231 (MQ=255)
                                                                                              aaTC‑GATATGAACCGCAGGCTTTCCCGGCAGGTAGACGCCATCGTCACGCAGTGCACACGGCCAGGGGGAGTAATATCCAGCGAGGCGCAGGTTGTAAGGGGTAAT‑CACCACCCCTTCCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGcc                                                                                                                          <  2:906441/238‑1 (MQ=255)
                                                                                               aTC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGCAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGAGAGGCGCGGTGCAGATGACAAAGGGGCCGCGTGTGTGTGTTTTAGTTGA‑GAAAACGaaaccacacccc                                                                                                                        >  2:367196/1‑229 (MQ=255)
                                                                                                 c‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GGAAAAGAACTcctgcccc                                                                                                                        >  1:464642/1‑230 (MQ=255)
                                                                                                     tatGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCCTAGCACAAGGTGGGCGGTTGTTGCTTTTTCAGTTGA‑GGTACAGAACTCACGCCCtcacca                                                                                                                   >  1:102518/1‑239 (MQ=255)
                                                                                                          aaCCCAGGCATTTCCGGCAGGTAGCCGTCACCGTCCCGCAGTGCACACTGCCAGGCGGAGTAATATCCAGGGAAGCGCAGGTTGTAATGGGTAAT‑CCCCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCggg                                                                                                              <  1:488526/239‑1 (MQ=255)
                                                                                                           ccGCAGGCATTACCGTCAGGTAGACGTCATCGTCCCGCAGTGCACACTGCCAGGCGGAGTAATACCCAGCGAAGGGCAGGTTGTAATTGGTAAT‑CCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGa                                                                                                             <  1:185626/238‑1 (MQ=255)
                                                                                                               aGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTGGTTGATTTTTCAGTTGG‑GTAAAAGAACTCACCCCATACTCACCGGGAAgg                                                                                                          >  1:669604/1‑237 (MQ=255)
                                                                                                                                     tcttcttcccggtgttcacctgcccAGGCGGTGTAATACCCGGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTCGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAAtgtg                                                                                  <  2:358393/216‑1 (MQ=255)
                                                                                                                                        gctgcccggaGTGCCAACTGCCGGGGGGAGTAATTTCCAGGGAGGGGGAGGTTGTAATGGGTAAC‑CACCCCTCCTTTCGGTCGGGAGGGGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCg                                                                               <  2:63220/230‑1 (MQ=255)
                                                                                                                                                   tGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCGGGTTTGTGTGTTTTCAAGTGATGTAACACAAC‑CCCCCCAACCTCACCGGGAAGTGGCCCCACTGTC‑AGGaaaactgtcccccattgaaggg                                                                    >  2:768410/1‑221 (MQ=255)
                                                                                                                                                        acTGCCAGGGGGATTAATTCCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCt                                                              <  1:510359/239‑1 (MQ=255)
                                                                                                                                                          tGCGAGGGGGAGTAATATCCGGGGGAGGGCGGGTTGTAAGGGGTAAC‑CCCCCCCCCTTTCGGGCGGGGGGGGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTc                                                            <  2:140164/239‑1 (MQ=255)
                                                                                                                                                                   gAGTATTATCCGGCGAGGGGCGGGTGGTAAGGGTTAAT‑CACCCCCCCGTTCGGGGGGGGGGGGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGGCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGaa                                                    <  2:669604/238‑1 (MQ=255)
                                                                                                                                                                    aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAAGTGGCCGGTGTTTGATTTTTTCGTTGA‑GTAAACGAACCCACGCCATCATCAGCGGGAAGTGCCCCACCTGTCAAGGCAAAT‑TGCCGCCAATTG‑GTGGGAGCTTccgcctggaaa                                                  >  2:195219/1‑230 (MQ=255)
                                                                                                                                                                     gtatgatCCAGCAAAGGCCAGGTTTGAATGGGTAGT‑CACCAAACCTTTGGGTCGGGAGGTGGTGCCGGAGTTGAAGAGAATTTACGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGTAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATc                                                  <  1:717923/233‑1 (MQ=255)
                                                                                                                                                                     gTAATATCCAGCGAGGCGCAGGTTGTAAGGGGTAAC‑CACCACCCCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATc                                                  <  1:826409/238‑1 (MQ=255)
                                                                                                                                                                     gTAATACCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATc                                                  <  1:254437/238‑1 (MQ=255)
                                                                                                                                                                       tatgccccggggggggcgggTTGT‑ATGGGGAATCCCCCCCACCTTTGGGGTGGGGGGTGTGGCGGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTTGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGt                                               <  2:108907/221‑1 (MQ=255)
                                                                                                                                                                       aattgccggggAAGGGCAGGTTTTAATGGGTAAT‑CACCACCCCTTTCGGGCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTg                                                <  1:662126/228‑1 (MQ=255)
                                                                                                                                                                         tacccggcgagggtcAGGTTTGAAGGGGTAAT‑GCCCCCCCCTGTCGGTGGGGGGGGGGTGCGGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGTGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTg                                             <  2:55759/225‑1 (MQ=255)
                                                                                                                                                                              tgCGAAGCTGAGTTTGTAATGGGTAAT‑CACCACACCTTTCGTTGGGGAGGTGGTGCCGGAGGTGAAGAGAATGTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:802699/238‑1 (MQ=255)
                                                                                                                                                                              ggggAAGCGGAGGTTGTAAGGGGAAAC‑CCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:876383/236‑1 (MQ=255)
                                                                                                                                                                              ggCGAAGCGCAGGTTGTAATGGGTAAG‑CCCCCCCCCTTCCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:499463/238‑1 (MQ=255)
                                                                                                                                                                              cgCGATGCGAAGGTTGTAATGGGTAAT‑CCCCTCTGCTTGCGGTCGGGAGGTGTTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:479318/238‑1 (MQ=255)
                                                                                                                                                                                 gaaggggggtttGTAAGGGGTAAT‑CACCACTCCTTTGGGTCGGGGGGTGGTGGCGGAGGTGAAGAAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  2:257119/226‑1 (MQ=255)
                                                                                                                                                                                 gAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:713925/236‑1 (MQ=255)
                                                                                                                                                                                 gAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTGTGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:738444/236‑1 (MQ=255)
                                                                                                                                                                                 agaGCCCAGGTTGTAAGGGGTAAT‑CACCACACCTTTCGGTCGGGAGGGGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:638/234‑1 (MQ=255)
                                                                                                                                                                                     cgcAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGCTAGCGGCGGGGCATAACCCAAGGGGGCGGGTTGTTGATTTTTTAGTTGA‑GTAAACGAACTCCACCCATCACCAACGGGAAGATGCACACCTTTC‑AGGCAAATATGCCC‑CCATTGAGGGGGCTATTCCTCCTCaaaatttttgtacataagaaa                                 >  2:738444/1‑220 (MQ=255)
                                                                                                                                                                                         ggTTGGAATGGGTAAC‑GCCCACACCTTTCGGTTGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGaa                              <  1:92925/238‑1 (MQ=255)
                                                                                                                                                                                                aaggggaaacccccccacCTTTGGGTCGGGAGGTGGTGCCGGAGGTGAAGGGAATTTCCGCCGTATGGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTgcgcg                       <  1:717039/224‑2 (MQ=255)
                                                                                                                                                                                                aTGGTTACT‑CCCCACCCTTTTGGGGCGGGGGGGGTTGCCGGGGGTGAAGAGAATTTCCGCCTTATCTGCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAGTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCAc                      <  2:102518/239‑1 (MQ=255)
                                                                                                                                                                                                      tattcccccctcttttGGGGCGGGGGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGt                  <  2:314014/225‑1 (MQ=255)
                                                                                                                                                                                                        t‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCC‑CAATTTAGTGGCATCTTCCTGCTGAAACTGTTTGTACATAGGGTAGAATTGCGCCCTGGGCAcc              >  1:245982/1‑239 (MQ=255)
                                                                                                                                                                                                               cgccTTTCGGGCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGGCGATAGCGGCGGTGCATAGCCCAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAgg        <  1:481078/237‑1 (MQ=255)
                                                                                                                                                                                                               cccttttgggggggggggggggGCGGGAGGTAGAGGAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCCCATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAgg        <  2:775960/227‑1 (MQ=255)
                                                                                                                                                                                                                     tcggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGGTTTTTCAGTTGA‑GTAAACGAACTCACGCCCTCATCAGCGGGAAGTGCCCCATCTGTC‑AGGCAAAATTGCCC‑CCATTGAGTGGGATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTTCGCACTGGTCAccctgtctctttta  >  1:644321/1‑227 (MQ=255)
                                                                                                                                                                                                                                     gCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTaaa                                                                                                                                      <  1:65947/93‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                         
TACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCG‑CCATCGC‑CGCAGTACACTGGCAATC‑GATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAT‑CACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGA‑GTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTC‑AGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCC  >  NC_000913/37051‑37498

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑