Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,250 A→C 21.8% V192G (GTG→GGG)  caiC ← putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091337,2500AC21.8% 113.4 / 3.8 65V192G (GTG→GGG) caiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base A (26/20);  new base C (3/15);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 5.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCG  >  NC_000913/37015‑37475
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
cgtattACGTTCGGTACGGGGGGCCCGTGACTTTTGGTCCTGTCCCCCGGAGGCGGGGGCTCTTTTTTTCCCGACCAGCACAAAGGTGGCCTCGGCAGAAAACGTCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATcacca                                                                                                                                                                                                                                   <  2:259245/234‑1 (MQ=255)
cGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTAATTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATcacca                                                                                                                                                                                                                                   <  1:510981/239‑1 (MQ=255)
          tCGGTAACGGTGGCGCGGTTCTTCCGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTtcgg                                                                                                                                                                                                                         <  1:494595/239‑1 (MQ=255)
             gTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGGGCAGGTTGTAATGGGTAATCACCCAACCTTtcggtcg                                                                                                                                                                                                                      >  1:139061/1‑239 (MQ=255)
               aaCGGTGGGGCTGTACTTCTGTACCGTTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCCAAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGGGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCggg                                                                                                                                                                                                                    <  2:562892/239‑1 (MQ=255)
                 cGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAAAATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGCAg                                                                                                                                                                                                                  >  1:39172/1‑239 (MQ=255)
                      gcgcGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTACTATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCCCACCTTTCGGTCGGGCggtggg                                                                                                                                                                                                             >  1:634360/1‑238 (MQ=255)
                        gcgGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAAAGCAGGCCTTACCGCCAGGTACACG‑TCATCGTCCCGCAGTGCACACTCCCAGGCGGAGTAATATCAAGCGAACCGCAGGTGGTAATGGGGAATACCAACACCTTTCGGTCGGGAgggggggg                                                                                                                                                                                                           >  2:436227/1‑236 (MQ=255)
                               ttCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGt                                                                                                                                                                                                    <  1:370842/239‑1 (MQ=255)
                                         gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGGGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGGGGTGGTGCCGGAGGTGAAAGGAAtt                                                                                                                                                                                          >  1:687246/1‑239 (MQ=255)
                                         gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAAtt                                                                                                                                                                                          >  1:501685/1‑239 (MQ=255)
                                         gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAGGAGAAtt                                                                                                                                                                                          >  1:510286/1‑239 (MQ=255)
                                         gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATACCCACACCTTTCGGTCGGGAGGGGGGGCCGgggggggaaggaaat                                                                                                                                                                                          >  2:52007/1‑228 (MQ=255)
                                           ccccAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCAAGCACTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTAGTAATGGGTAATCACCACACCTTTCGGACGGGGGGTGGGGCGGGGGGGGAAGAGAAGTTc                                                                                                                                                                                        >  2:122110/1‑239 (MQ=255)
                                            cccAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTc                                                                                                                                                                                        >  1:150734/1‑238 (MQ=255)
                                              cAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGGAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTccg                                                                                                                                                                                      >  1:484423/1‑238 (MQ=255)
                                                 aaGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGt                                                                                                                                                                                  >  1:788140/1‑239 (MQ=255)
                                                    gcgcGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGGGCCGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGGGGGGGGGCCGGAGGGGGAGAGAAGTTccgccggagg                                                                                                                                                                               >  2:452554/1‑235 (MQ=255)
                                                        gggCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGGGCACACTGCCAGGCGGAGTAAAATCCAGCGAAGACCAGGTTGAAATGGGAAATCCCCACACCTTTCGGTCGGGGGGCGGGGCCGGAGGTGAAGAGAAATTCCACCATATCGtaa                                                                                                                                                                           >  2:260018/1‑237 (MQ=255)
                                                                     ctcGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGAAATCACCCAACCTTTCGGTCGGGAGGTGGTCCCGGAGGTGACGAGAATATCCGCAGTATCGTCAGTCGATAgcgtcg                                                                                                                                                              >  1:847176/1‑239 (MQ=255)
                                                                     ctcGACCAGCAAAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAgcggcg                                                                                                                                                              <  1:715794/239‑1 (MQ=255)
                                                                          ccaccaCAAGGGGGGCCCCGGCAGAAACCGCCGCCATCGCCGCAGTCCACGGGGAACCGATTGAAAACCCGGGATTTACCGGAAGGT‑GACGTTCATCGGCACGCAGTGCACCCTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCCCACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                         <  2:24123/239‑1 (MQ=255)
                                                                          ccaccaCAAAGGGGGCCCCGGAAAAAACCGCCCCCATCGCCGCAGTACACTGGCAATCGATAGGAAACGCAGGCATTACCGTCAGGTAGGCG‑TCATCGTCACGCAGTGCACACTGCCGGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAAGGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                         <  2:139061/239‑1 (MQ=255)
                                                                          ccAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAAGGGAAGAGAATTTCCCCCGGATCGTCAGTCGAAACcggcgggggc                                                                                                                                                         >  1:625328/1‑235 (MQ=255)
                                                                                         ccccGGCAGAAAACGCCGCCAGCGCGGCAGTACACTGGAAATCGTTATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACGGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGc                                                                                                                                           <  1:298931/238‑1 (MQ=255)
                                                                                                         cgccATCGCCGCAGTACACTGGCAAGCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCAAACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCa                                                                                                                         <  1:50758/239‑1 (MQ=255)
                                                                                                             aTCCCCGCAGTACCCGGGCAATCGATAGGAAACGCCGGCTTATCCGTCAGGTAGCGG‑TCCCCATCACGCAGTGCACACTGCCAGGCGAAGTAATATCCAGGGAAGCGCGGGTTGTAAGGGGTAATCACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTg                                                                                                                     <  2:510286/239‑1 (MQ=255)
                                                                                                                  cgcAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCGTAGGACAAGGGGGCAGGTTTGTGA‑TTTTTTAGTTGAggga                                                                                                                >  1:823867/1‑236 (MQ=255)
                                                                                                                                 aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCGGGGTGTTGG‑GTTTTTAGTTGGGTCAACAGACTccccccc                                                                                                  >  1:162099/1‑232 (MQ=255)
                                                                                                                                 aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCCTAGCCCAAGGGGGCCGGTGTTTGAGTTTTTCTTTTG‑GTAAACGAACTCACGccc                                                                                                  >  1:188641/1‑237 (MQ=255)
                                                                                                                                                             aGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGGACAGCCCAAAGGGGCAGGGTGTTGA‑TTTTTTAGTTGAGTAAACGAACTCACGCCCTC‑ATCAGCGGGAAATGCCCAATCTGTccg                                                                    >  1:99090/1‑237 (MQ=255)
                                                                                                                                                                      ggcATCGTCACGCAGTGCACACGGCCAGGCGGAGTAATATCCAGCGAAGGGCAGGTTGTAATGGGTAATCCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAAtg                                                             <  1:100717/236‑1 (MQ=255)
                                                                                                                                                                          tCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATAGCACAAGGGGGCAGGGTGTTGA‑TTTTTCTGTTTAGTAAAAGAAATCACCCCACC‑ATCAACGGGGAGTGGCACAACTGTCAGGaaaaaggggcg                                                        >  2:646764/1‑232 (MQ=255)
                                                                                                                                                                             tcccGCAGTGCACACTGCAAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑Caa                                                    <  1:680417/236‑1 (MQ=255)
                                                                                                                                                                               aCGGAGGGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAAGGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAAtt                                                  <  1:421/239‑1 (MQ=255)
                                                                                                                                                                                  gcgtgccactcgCCAGGGCGAGTAATTCCCCGCGGAGGGCAGGTTGAAATGGGTAACCCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCAAAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAg                                               <  1:242298/228‑1 (MQ=255)
                                                                                                                                                                                  catgtcaccCTGCCGGGGGGAGTAATTTCCAGCGAAGCTCAGGTTGTAATGGGTAAACCCCAACCCTTTCGGGCGGAAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAg                                               <  1:377379/234‑1 (MQ=255)
                                                                                                                                                                                            tGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATAGAACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATAAGCGGGGAAGGCCACATCTGTTAGGGAAATGTGCCG‑CAATTGGGTGGGATAtt                                      >  1:873661/1‑238 (MQ=255)
                                                                                                                                                                                             gCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGCAGAGAATTTCCGCCGTAGCGTCCGTCGATAGCGGCGGTGCAGAGCACAAAGGGGCAGGTGTTTGT‑TTTTTTAGTTTGGGAAAGCAACCACCGCCAAC‑ATTACCGGGGAGGGCTAAATCTTTTCGGCAAAAGTTCCGCCAATTG‑GGGGCATCtcc                                     >  2:836281/1‑236 (MQ=255)
                                                                                                                                                                                                  gcgggtttttttCCAGGTAAGCTCAGGTTGTAAGGGGTAATCCCCCCCCCTTTCGGGCGGGGGGTGGTGCGGGAGGTAAAGGGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCctgct                                <  2:99090/227‑1 (MQ=255)
                                                                                                                                                                                                      aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATTTGCCG‑CAATTGAGTGGGACTTTCCTGCTGaa                             >  1:410456/1‑237 (MQ=255)
                                                                                                                                                                                                      aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGGTGTTGA‑TTTTTCAGTTTAGTAAACGAACTCACGCCACC‑ATCAGCGGGAAGGGCCACATCTGTAAGGCAAAAGTGCCG‑CAATTGAGGGGGATCTTCCTGCTGGAta                           >  1:289945/1‑238 (MQ=255)
                                                                                                                                                                                                      aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGACGATAGCGGCGGGGCCCAGACAAAGGGGGCAGGTTGGTGA‑TTTTTCCATTGAGTTAACGAAATCAAGCCACC‑CACCACAGGAGGTGCCCCACCTGTCAGGCAAAAGTGCCG‑CAATTGAGGGGCGACTTACTGCTGaaaa                           >  2:44401/1‑237 (MQ=255)
                                                                                                                                                                                                       ggaaattCCAGCGAGGGGCAGTTTGTATTGGGTAATCCCCACACCTTGCGGTCGGGAGGTGGTGCCGGAGGTAAAGGGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACAGCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATc                           <  1:415831/232‑1 (MQ=255)
                                                                                                                                                                                                       gTAATTGCCAGCGAAGCGAAGGTTGGAATGGGAAATCACCCCCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATc                           <  1:260018/238‑1 (MQ=255)
                                                                                                                                                                                                          atttCCAGGGAAGCGCAGGTTGAAATGGGTAATCACCCCCCCTTTCGGTCGGGAGATGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGt                        <  1:436227/235‑1 (MQ=255)
                                                                                                                                                                                                             cccAGCGAAGCTCAGGTGGTAAGTGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGc                                                          <  1:190951/201‑1 (MQ=255)
                                                                                                                                                                                                                ggggAAGCGGAGGTTGTAAGGGGTAATCCCCCCCCCTTGCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:339987/236‑1 (MQ=255)
                                                                                                                                                                                                                ggggAAGCGCAGGTTGTAATGGGTAATCACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:452554/236‑1 (MQ=255)
                                                                                                                                                                                                                gcgggggggcAGGTTGTAATGGGTAATGCCCCCCCCTTTCGGTGGGGAGGGGGTGGCGGGGGTGAAGAAAATTTCCGCCGTTTCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  2:873661/231‑1 (MQ=255)
                                                                                                                                                                                                                cgggaggcggAGTTTGTAATGGGTAATCCCCCGCCCTTGCGGGCGGGGGGTGGGTCCGGGGGTGAAGAAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  2:410456/233‑1 (MQ=255)
                                                                                                                                                                                                                cgcggaggggaGGTTGAAAGGGGTAATCCCCACCCTTTTGGGTCGGGGGGTGGGGCGGGGGGTAAAGGGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  2:447256/229‑1 (MQ=255)
                                                                                                                                                                                                                aGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGGAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCCCATCTTTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTTGTACa                 >  1:400776/1‑239 (MQ=255)
                                                                                                                                                                                                                aGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGGATAGCCCAAGGGGGGCGGGTGTTGA‑GTTTTTAGTTGAGTAAACGAACTTCCGCCAACAACCACCGGGAGGTGCC‑CCTCTGTCCGGGCAAAGTGGCG‑CCAATTTGTGGCCTCTCCCTTCCGAATCTGTTgtataa                 >  2:292120/1‑234 (MQ=255)
                                                                                                                                                                                                                   gaaggggggtttGTAAGGGGTAATCACCCCCCCTTTCGGGCGGGAGGGGGTGCCGGGGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  2:150734/226‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTTTAATGGGTAACCCCCCCCCCTTTCGGTCGGGAGGGTGTGCCGGAGGTTAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:368109/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAAGTGAAGAGAATTTGCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:1267/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTGTAATGGGTAAGCACCCCACCTTTGGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:572510/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTGTAATGGGTAACCTCCACCCCGTTCGGGCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:351381/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTGAAATGGGTAATCCCCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:52007/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCCCAGGTTGAAATGGGTAATCACCACCCCTTTGGGGCGGGAGGGGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:113281/236‑1 (MQ=255)
                                                                                                                                                                                                                      gccgagttgttAAGGGGTAATCCCCCCTCCTGGCGGGCGGGAGGTGGGGCCGGGGGTGAGGAGAATTGCCGCCGTAACGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGa            <  2:188641/228‑1 (MQ=255)
                                                                                                                                                                                                                       cgcAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGGTGTTTA‑TTTTTCAGTTGGGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGGGCCCCATTTTTCAGGCAAATGTGCCC‑CAATTGAGTGGCATCTTCCTGCTTAATCTTTTGATAACTAGGATa          >  1:447256/1‑239 (MQ=255)
                                                                                                                                                                                                                        ggaGGTTGAAATGGGTAAACCCCACCCCTTTCGGTCGGGAGGGGGGCCCGGAGGTGAAGAGAATTGCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATa          <  1:451267/236‑1 (MQ=255)
                                                                                                                                                                                                                          ggtttttAATGGGTAATCCCCACACCTTGCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGaa       <  2:289945/233‑1 (MQ=255)
                                                                                                                                                                                                                             gtttaagtggTAAGCCCCACCCCTTTCGGGCGGGAGGTGGGGCCGGGGGTGAAGGGAATTTCCGCCGGATCGGCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAAtt     <  1:321965/230‑1 (MQ=255)
                                                                                                                                                                                                                               gTAATGGGAAATCACCACACCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTgcg  <  1:63613/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
CGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACG‑TCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGA‑TTTTTCAGTTGAGTAAACGAACTCACGCCATC‑ATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCG‑CAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCG  >  NC_000913/37015‑37475

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