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

TACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCC  >  NC_000913/37051‑37498
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                         
tACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTAAACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACCCGGCCAGGCGGAGTAATATCCAGCGCAGCGCAGGATGTCATGGGTAATCACCCCACCATCCGGTCGGGCGGTGGTGCCGGGGGGGAagg                                                                                                                                                                                                                      >  2:499463/1‑238 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTCTTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCAACGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCCGGTAGACGTCATCCTCCCGCAGTGCACACTGCCAGACCGAGTAATATCCATCGCACCGCAGGTTGTAATGGGTAATCACCAAACCTTTCTGGCGGGAGGTGGGGCGGGGGGGGAAggga                                                                                                                                                                                                                    >  2:481078/1‑239 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACCCCTTTCGGTCGGCAGGGGGGGCCGGGGGTGAAgaga                                                                                                                                                                                                                    >  1:63220/1‑239 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGGGGTGGTGCCGGAGGGGAAgaga                                                                                                                                                                                                                    >  1:212264/1‑239 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAgaga                                                                                                                                                                                                                    >  1:358393/1‑239 (MQ=255)
  ccTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGGAATCACACCACCTTTCGGTCGGGAGGGGGGGCGGGGGGGGaaaaaa                                                                                                                                                                                                                    >  2:254437/1‑235 (MQ=255)
     gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGGGCCGGCGGTGAAGAGAAtt                                                                                                                                                                                                                 >  1:55759/1‑239 (MQ=255)
     gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACAGCCAGGCGGAGTAATATCCAGCGAAGCGCCGGGTGTAATGGGAAATCACCACACCATTCGGTCGGGAGGTGGGGCCGGAGGGGAAGAAAAtt                                                                                                                                                                                                                 >  2:185626/1‑239 (MQ=255)
     gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGGAGACGTCCTCTTCACGCAGTGCACACTGCCAGACGGCGGAATATCCACCGAAGCGCCGGTTGTAATGGGTAATAACCACCACCTTCGGTCGGGAGGGGGGGCCGGAGGGGAAAGGAAtt                                                                                                                                                                                                                 >  2:115046/1‑239 (MQ=255)
      tCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGAACACTGCCAGGCGGAGTAATATCCCGCGAGGCGCAGGTTGTAATGGGTAATCACCAAACCTTTCGGTCGGGAGGTGGGGCCGGGGGGGGagaggaat                                                                                                                                                                                                                 >  2:826409/1‑234 (MQ=255)
             aaGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCCGGTTGTAATGGGGAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGt                                                                                                                                                                                                         >  1:532869/1‑239 (MQ=255)
                gcgcGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGCTATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGGCACGCAGTGCACACTGCCAGGCGGGGGAATATCCAGCGAAGCGCAGGTTGTAATGGGTAAGAACCACACCTTTCGGTCGGGCGGTGGGGCCGGAGGGGGAAAGAATTTcacccggctc                                                                                                                                                                                                      >  2:96052/1‑230 (MQ=255)
                gcgcGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCGTATc                                                                                                                                                                                                      >  1:314014/1‑239 (MQ=255)
                 cgcgGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCCGGGTGGAATGGGTAATCACCCCACCCTTAGGTCGGGAGGGGGGGCCGGGGGGGAAGAGAATATCCCCCGTATCg                                                                                                                                                                                                     >  2:795337/1‑239 (MQ=255)
                    gggCGCTGTATTTCGCGACCAGCACAAAGGTGGCCCCGGCAGAAAACCCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGGCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCa                                                                                                                                                                                                  <  1:51563/239‑1 (MQ=255)
                            tttttCCCGCCCAGCAAGAAGGGGCCCCCGGCAAAAAACCCCCCCACCGCGGCAGTACCCTGGCAACCGATATGAACCGCAGGCATTCCCATCAGGTAGCCGTCATCGTCCCGCAGGGCACACGGCCAGGCGGAGTAATATCCAGCGAAGCGCAGTTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAg                                                                                                                                                                                          <  2:575293/237‑1 (MQ=255)
                             aTTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGAAATCGATATGAACCGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAg                                                                                                                                                                                          <  1:790995/238‑1 (MQ=255)
                                 ctcGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTccatgggaaagggcg                                                                                                                                                                                     >  2:790995/1‑225 (MQ=255)
                                  ccacgccgacaaaaGGGGGCCCCGGGAAAAAACCCCCCCCCTGCGCGGAGGACCATGGCAAGCGGATAGAAACCGAGGAATTCCCGGAGGGAAGACGTCCTCGTCACGCAGTGCCCCCTGCAAGGCGGAGTAATACCCAGCGAAGCGGAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGGCCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAgcggcg                                                                                                                                                                                     <  2:180724/228‑1 (MQ=255)
                                      ccaccaCAAAGGTACCCCGGGAAAAAAACCCCGCCATCCCCGCAGTACCCTGGCAATCGAAATGAAACGAAGGCATTACCGTAAGGTAGCCGTCAGCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                                                <  2:817500/239‑1 (MQ=255)
                                      ccAGCAAAAGGGTGGCCCCGGAAAAAACCGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                                                <  1:866592/239‑1 (MQ=255)
                                      caagaacaGAGGTGGCCCCGGCAAAAAACGGCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                                                <  1:795337/234‑1 (MQ=255)
                                         gCAAAATGGTGGCCCCGGCAGAAACCGCCGCCATCGCCGCAGTACACTGGCAACCGATATGAACCGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAg                                                                                                                                                                             <  1:390756/239‑1 (MQ=255)
                                         cccaacaGGTCGCCCCGGAAAAAACCCCGCCCAGCGCGCCGGTCCAGGGGAAATCGATATGAAACGCAGGCATTCCGGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAg                                                                                                                                                                             <  2:202946/233‑1 (MQ=255)
                                          cacaAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGGCAGTCGATAGCGGCGGTGGATAGc                                                                                                                                                                            >  1:56969/1‑239 (MQ=255)
                                          cacaAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTGGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTAACGTCAGTCGATACcggcgggccgaag                                                                                                                                                                             >  1:747298/1‑231 (MQ=255)
                                            caAAGGGTGCCCCGGCAGAAAACGCCCCCATGCCCGCAGTACCCTGGCAATCGATATGAAACGCAGGCATTCCCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGcac                                                                                                                                                                          <  1:760655/239‑1 (MQ=255)
                                                         ggCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGAAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGGGCCGGGGGGGAAGAGAATTTCCGCCGTTACGGCCGTCGAAAGgggggggggggaacccaacaggggggggg                                                                                                                                                             >  2:802276/1‑216 (MQ=255)
                                                           cAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGGGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATTGCCAAAGgggggggggggt                                                                                                                                                           >  1:257119/1‑231 (MQ=255)
                                                                                          tGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATCGCGGCGGTGCATCGCCCAAGGAGGCGGGATGTTGATTTTTCCCGTGAGGAAAcacaccaa                                                                                                                              >  1:910764/1‑231 (MQ=255)
                                                                                              aaTCGATATGAACCGCAGGCTTTCCCGGCAGGTAGACGCCATCGTCACGCAGTGCACACGGCCAGGGGGAGTAATATCCAGCGAGGCGCAGGTTGTAAGGGGTAATCACCACCCCTTCCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGcc                                                                                                                          <  2:906441/238‑1 (MQ=255)
                                                                                               aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGCAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGAGAGGCGCGGTGCAGATGACAAAGGGGCCGCGTGTGTGTGTTTTAGTTGAGAAAACGaaaccacacccc                                                                                                                        >  2:367196/1‑229 (MQ=255)
                                                                                                 cGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGGAAAAGAACTcctgcccc                                                                                                                        >  1:464642/1‑230 (MQ=255)
                                                                                                     tatGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCCTAGCACAAGGTGGGCGGTTGTTGCTTTTTCAGTTGAGGTACAGAACTCACGCCCtcacca                                                                                                                   >  1:102518/1‑239 (MQ=255)
                                                                                                          aaCCCAGGCATTTCCGGCAGGTAGCCGTCACCGTCCCGCAGTGCACACTGCCAGGCGGAGTAATATCCAGGGAAGCGCAGGTTGTAATGGGTAATCCCCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCggg                                                                                                              <  1:488526/239‑1 (MQ=255)
                                                                                                           ccGCAGGCATTACCGTCAGGTAGACGTCATCGTCCCGCAGTGCACACTGCCAGGCGGAGTAATACCCAGCGAAGGGCAGGTTGTAATTGGTAATCCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGa                                                                                                             <  1:185626/238‑1 (MQ=255)
                                                                                                               aGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTGGTTGATTTTTCAGTTGGGTAAAAGAACTCACCCCATACTCACCGGGAAgg                                                                                                          >  1:669604/1‑237 (MQ=255)
                                                                                                                                     tcttcttcccggtgttcacctgcccAGGCGGTGTAATACCCGGCGAAGCGCAGGTTGTAATGGGTAATCACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTCGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAtgtg                                                                                  <  2:358393/216‑1 (MQ=255)
                                                                                                                                        gctgcccggaGTGCCAACTGCCGGGGGGAGTAATTTCCAGGGAGGGGGAGGTTGTAATGGGTAACCACCCCTCCTTTCGGTCGGGAGGGGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                                                               <  2:63220/230‑1 (MQ=255)
                                                                                                                                                   tGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCGGGTTTGTGTGTTTTCAAGTGATGTAACACAACCCCCCCAACCTCACCGGGAAGTGGCCCCACTGTCAGGaaaactgtcccccattgaaggg                                                                    >  2:768410/1‑221 (MQ=255)
                                                                                                                                                        acTGCCAGGGGGATTAATTCCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCt                                                              <  1:510359/239‑1 (MQ=255)
                                                                                                                                                          tGCGAGGGGGAGTAATATCCGGGGGAGGGCGGGTTGTAAGGGGTAACCCCCCCCCCTTTCGGGCGGGGGGGGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTc                                                            <  2:140164/239‑1 (MQ=255)
                                                                                                                                                                   gAGTATTATCCGGCGAGGGGCGGGTGGTAAGGGTTAATCACCCCCCCGTTCGGGGGGGGGGGGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGGCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGaa                                                    <  2:669604/238‑1 (MQ=255)
                                                                                                                                                                    aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAAGTGGCCGGTGTTTGATTTTTTCGTTGAGTAAACGAACCCACGCCATCATCAGCGGGAAGTGCCCCACCTGTCAAGGCAAATTGCCGCCAATTGGTGGGAGCTTccgcctggaaa                                                  >  2:195219/1‑230 (MQ=255)
                                                                                                                                                                     gtatgatCCAGCAAAGGCCAGGTTTGAATGGGTAGTCACCAAACCTTTGGGTCGGGAGGTGGTGCCGGAGTTGAAGAGAATTTACGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGTAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc                                                  <  1:717923/233‑1 (MQ=255)
                                                                                                                                                                     gTAATATCCAGCGAGGCGCAGGTTGTAAGGGGTAACCACCACCCCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc                                                  <  1:826409/238‑1 (MQ=255)
                                                                                                                                                                     gTAATACCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc                                                  <  1:254437/238‑1 (MQ=255)
                                                                                                                                                                       tatgccccggggggggcgggTTGTATGGGGAATCCCCCCCACCTTTGGGGTGGGGGGTGTGGCGGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTTGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGt                                               <  2:108907/221‑1 (MQ=255)
                                                                                                                                                                       aattgccggggAAGGGCAGGTTTTAATGGGTAATCACCACCCCTTTCGGGCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg                                                <  1:662126/228‑1 (MQ=255)
                                                                                                                                                                         tacccggcgagggtcAGGTTTGAAGGGGTAATGCCCCCCCCTGTCGGTGGGGGGGGGGTGCGGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGTGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg                                             <  2:55759/225‑1 (MQ=255)
                                                                                                                                                                              tgCGAAGCTGAGTTTGTAATGGGTAATCACCACACCTTTCGTTGGGGAGGTGGTGCCGGAGGTGAAGAGAATGTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:802699/238‑1 (MQ=255)
                                                                                                                                                                              ggggAAGCGGAGGTTGTAAGGGGAAACCCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:876383/236‑1 (MQ=255)
                                                                                                                                                                              ggCGAAGCGCAGGTTGTAATGGGTAAGCCCCCCCCCTTCCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:499463/238‑1 (MQ=255)
                                                                                                                                                                              cgCGATGCGAAGGTTGTAATGGGTAATCCCCTCTGCTTGCGGTCGGGAGGTGTTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:479318/238‑1 (MQ=255)
                                                                                                                                                                                 gaaggggggtttGTAAGGGGTAATCACCACTCCTTTGGGTCGGGGGGTGGTGGCGGAGGTGAAGAAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  2:257119/226‑1 (MQ=255)
                                                                                                                                                                                 gAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:713925/236‑1 (MQ=255)
                                                                                                                                                                                 gAAGCGCAGGTTGTAATGGGTAATCACCACACCTTGTGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:738444/236‑1 (MQ=255)
                                                                                                                                                                                 agaGCCCAGGTTGTAAGGGGTAATCACCACACCTTTCGGTCGGGAGGGGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                                        <  1:638/234‑1 (MQ=255)
                                                                                                                                                                                     cgcAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGCTAGCGGCGGGGCATAACCCAAGGGGGCGGGTTGTTGATTTTTTAGTTGAGTAAACGAACTCCACCCATCACCAACGGGAAGATGCACACCTTTCAGGCAAATATGCCCCCATTGAGGGGGCTATTCCTCCTCaaaatttttgtacataagaaa                                 >  2:738444/1‑220 (MQ=255)
                                                                                                                                                                                         ggTTGGAATGGGTAACGCCCACACCTTTCGGTTGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGaa                              <  1:92925/238‑1 (MQ=255)
                                                                                                                                                                                                aaggggaaacccccccacCTTTGGGTCGGGAGGTGGTGCCGGAGGTGAAGGGAATTTCCGCCGTATGGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGCGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTgcgcg                       <  1:717039/224‑2 (MQ=255)
                                                                                                                                                                                                aTGGTTACTCCCCACCCTTTTGGGGCGGGGGGGGTTGCCGGGGGTGAAGAGAATTTCCGCCTTATCTGCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAGTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCAc                      <  2:102518/239‑1 (MQ=255)
                                                                                                                                                                                                      tattcccccctcttttGGGGCGGGGGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGt                  <  2:314014/225‑1 (MQ=255)
                                                                                                                                                                                                        tCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCCCAATTTAGTGGCATCTTCCTGCTGAAACTGTTTGTACATAGGGTAGAATTGCGCCCTGGGCAcc              >  1:245982/1‑239 (MQ=255)
                                                                                                                                                                                                               cgccTTTCGGGCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGGCGATAGCGGCGGTGCATAGCCCAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAgg        <  1:481078/237‑1 (MQ=255)
                                                                                                                                                                                                               cccttttgggggggggggggggGCGGGAGGTAGAGGAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCCCATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAgg        <  2:775960/227‑1 (MQ=255)
                                                                                                                                                                                                                     tcggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGGTTTTTCAGTTGAGTAAACGAACTCACGCCCTCATCAGCGGGAAGTGCCCCATCTGTCAGGCAAAATTGCCCCCATTGAGTGGGATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTTCGCACTGGTCAccctgtctctttta  >  1:644321/1‑227 (MQ=255)
                                                                                                                                                                                                                                     gCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaaa                                                                                                                                      <  1:65947/93‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                         
TACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCC  >  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: