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

CGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCG  >  NC_000913/37015‑37475
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
cgtattACGTTCGGTACGGGGGGCCCGTGACTTTTGGTCCTGTCCCCCGGAGGCGGGGGCTCTTTTTTTCCCGACCAGCACAAAGGTGGCCTCGGCAGAAAACGTCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATcacca                                                                                                                                                                                                                                   <  2:259245/234‑1 (MQ=255)
cGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTAATTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATcacca                                                                                                                                                                                                                                   <  1:510981/239‑1 (MQ=255)
          tCGGTAACGGTGGCGCGGTTCTTCCGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTtcgg                                                                                                                                                                                                                         <  1:494595/239‑1 (MQ=255)
             gTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGGGCAGGTTGTAATGGGTAATCACCCAACCTTtcggtcg                                                                                                                                                                                                                      >  1:139061/1‑239 (MQ=255)
               aaCGGTGGGGCTGTACTTCTGTACCGTTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCCAAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGGGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCggg                                                                                                                                                                                                                    <  2:562892/239‑1 (MQ=255)
                 cGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAAAATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGCAg                                                                                                                                                                                                                  >  1:39172/1‑239 (MQ=255)
                      gcgcGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTACTATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCCCACCTTTCGGTCGGGCggtggg                                                                                                                                                                                                             >  1:634360/1‑238 (MQ=255)
                        gcgGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAAAGCAGGCCTTACCGCCAGGTACACGTCATCGTCCCGCAGTGCACACTCCCAGGCGGAGTAATATCAAGCGAACCGCAGGTGGTAATGGGGAATACCAACACCTTTCGGTCGGGAgggggggg                                                                                                                                                                                                           >  2:436227/1‑236 (MQ=255)
                               ttCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGt                                                                                                                                                                                                    <  1:370842/239‑1 (MQ=255)
                                         gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGGGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGGGGTGGTGCCGGAGGTGAAAGGAAtt                                                                                                                                                                                          >  1:687246/1‑239 (MQ=255)
                                         gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAAtt                                                                                                                                                                                          >  1:501685/1‑239 (MQ=255)
                                         gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAGGAGAAtt                                                                                                                                                                                          >  1:510286/1‑239 (MQ=255)
                                         gTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATACCCACACCTTTCGGTCGGGAGGGGGGGCCGgggggggaaggaaat                                                                                                                                                                                          >  2:52007/1‑228 (MQ=255)
                                           ccccAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCAAGCACTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTAGTAATGGGTAATCACCACACCTTTCGGACGGGGGGTGGGGCGGGGGGGGAAGAGAAGTTc                                                                                                                                                                                        >  2:122110/1‑239 (MQ=255)
                                            cccAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTc                                                                                                                                                                                        >  1:150734/1‑238 (MQ=255)
                                              cAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGGAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTccg                                                                                                                                                                                      >  1:484423/1‑238 (MQ=255)
                                                 aaGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGt                                                                                                                                                                                  >  1:788140/1‑239 (MQ=255)
                                                    gcgcGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGGGCCGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGGGGGGGGGCCGGAGGGGGAGAGAAGTTccgccggagg                                                                                                                                                                               >  2:452554/1‑235 (MQ=255)
                                                        gggCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGGGCACACTGCCAGGCGGAGTAAAATCCAGCGAAGACCAGGTTGAAATGGGAAATCCCCACACCTTTCGGTCGGGGGGCGGGGCCGGAGGTGAAGAGAAATTCCACCATATCGtaa                                                                                                                                                                           >  2:260018/1‑237 (MQ=255)
                                                                     ctcGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGAAATCACCCAACCTTTCGGTCGGGAGGTGGTCCCGGAGGTGACGAGAATATCCGCAGTATCGTCAGTCGATAgcgtcg                                                                                                                                                              >  1:847176/1‑239 (MQ=255)
                                                                     ctcGACCAGCAAAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAgcggcg                                                                                                                                                              <  1:715794/239‑1 (MQ=255)
                                                                          ccaccaCAAGGGGGGCCCCGGCAGAAACCGCCGCCATCGCCGCAGTCCACGGGGAACCGATTGAAAACCCGGGATTTACCGGAAGGTGACGTTCATCGGCACGCAGTGCACCCTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCCCACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                         <  2:24123/239‑1 (MQ=255)
                                                                          ccaccaCAAAGGGGGCCCCGGAAAAAACCGCCCCCATCGCCGCAGTACACTGGCAATCGATAGGAAACGCAGGCATTACCGTCAGGTAGGCGTCATCGTCACGCAGTGCACACTGCCGGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAAGGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                                         <  2:139061/239‑1 (MQ=255)
                                                                          ccAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAAGGGAAGAGAATTTCCCCCGGATCGTCAGTCGAAACcggcgggggc                                                                                                                                                         >  1:625328/1‑235 (MQ=255)
                                                                                         ccccGGCAGAAAACGCCGCCAGCGCGGCAGTACACTGGAAATCGTTATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACGGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGc                                                                                                                                           <  1:298931/238‑1 (MQ=255)
                                                                                                         cgccATCGCCGCAGTACACTGGCAAGCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCAAACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCa                                                                                                                         <  1:50758/239‑1 (MQ=255)
                                                                                                             aTCCCCGCAGTACCCGGGCAATCGATAGGAAACGCCGGCTTATCCGTCAGGTAGCGGTCCCCATCACGCAGTGCACACTGCCAGGCGAAGTAATATCCAGGGAAGCGCGGGTTGTAAGGGGTAATCACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTg                                                                                                                     <  2:510286/239‑1 (MQ=255)
                                                                                                                  cgcAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCGTAGGACAAGGGGGCAGGTTTGTGATTTTTTAGTTGAggga                                                                                                                >  1:823867/1‑236 (MQ=255)
                                                                                                                                 aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGGGGCGGGGTGTTGGGTTTTTAGTTGGGTCAACAGACTccccccc                                                                                                  >  1:162099/1‑232 (MQ=255)
                                                                                                                                 aTCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCCTAGCCCAAGGGGGCCGGTGTTTGAGTTTTTCTTTTGGTAAACGAACTCACGccc                                                                                                  >  1:188641/1‑237 (MQ=255)
                                                                                                                                                             aGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGGACAGCCCAAAGGGGCAGGGTGTTGATTTTTTAGTTGAGTAAACGAACTCACGCCCTCATCAGCGGGAAATGCCCAATCTGTccg                                                                    >  1:99090/1‑237 (MQ=255)
                                                                                                                                                                      ggcATCGTCACGCAGTGCACACGGCCAGGCGGAGTAATATCCAGCGAAGGGCAGGTTGTAATGGGTAATCCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAtg                                                             <  1:100717/236‑1 (MQ=255)
                                                                                                                                                                          tCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATAGCACAAGGGGGCAGGGTGTTGATTTTTCTGTTTAGTAAAAGAAATCACCCCACCATCAACGGGGAGTGGCACAACTGTCAGGaaaaaggggcg                                                        >  2:646764/1‑232 (MQ=255)
                                                                                                                                                                             tcccGCAGTGCACACTGCAAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCaa                                                    <  1:680417/236‑1 (MQ=255)
                                                                                                                                                                               aCGGAGGGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAAGGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAAtt                                                  <  1:421/239‑1 (MQ=255)
                                                                                                                                                                                  gcgtgccactcgCCAGGGCGAGTAATTCCCCGCGGAGGGCAGGTTGAAATGGGTAACCCCCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCAAAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAg                                               <  1:242298/228‑1 (MQ=255)
                                                                                                                                                                                  catgtcaccCTGCCGGGGGGAGTAATTTCCAGCGAAGCTCAGGTTGTAATGGGTAAACCCCAACCCTTTCGGGCGGAAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAg                                               <  1:377379/234‑1 (MQ=255)
                                                                                                                                                                                            tGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGCATAGAACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATAAGCGGGGAAGGCCACATCTGTTAGGGAAATGTGCCGCAATTGGGTGGGATAtt                                      >  1:873661/1‑238 (MQ=255)
                                                                                                                                                                                             gCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGCAGAGAATTTCCGCCGTAGCGTCCGTCGATAGCGGCGGTGCAGAGCACAAAGGGGCAGGTGTTTGTTTTTTTAGTTTGGGAAAGCAACCACCGCCAACATTACCGGGGAGGGCTAAATCTTTTCGGCAAAAGTTCCGCCAATTGGGGGCATCtcc                                     >  2:836281/1‑236 (MQ=255)
                                                                                                                                                                                                  gcgggtttttttCCAGGTAAGCTCAGGTTGTAAGGGGTAATCCCCCCCCCTTTCGGGCGGGGGGTGGTGCGGGAGGTAAAGGGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgct                                <  2:99090/227‑1 (MQ=255)
                                                                                                                                                                                                      aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATTTGCCGCAATTGAGTGGGACTTTCCTGCTGaa                             >  1:410456/1‑237 (MQ=255)
                                                                                                                                                                                                      aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGGTGTTGATTTTTCAGTTTAGTAAACGAACTCACGCCACCATCAGCGGGAAGGGCCACATCTGTAAGGCAAAAGTGCCGCAATTGAGGGGGATCTTCCTGCTGGAta                           >  1:289945/1‑238 (MQ=255)
                                                                                                                                                                                                      aGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGACGATAGCGGCGGGGCCCAGACAAAGGGGGCAGGTTGGTGATTTTTCCATTGAGTTAACGAAATCAAGCCACCCACCACAGGAGGTGCCCCACCTGTCAGGCAAAAGTGCCGCAATTGAGGGGCGACTTACTGCTGaaaa                           >  2:44401/1‑237 (MQ=255)
                                                                                                                                                                                                       ggaaattCCAGCGAGGGGCAGTTTGTATTGGGTAATCCCCACACCTTGCGGTCGGGAGGTGGTGCCGGAGGTAAAGGGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACAGCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc                           <  1:415831/232‑1 (MQ=255)
                                                                                                                                                                                                       gTAATTGCCAGCGAAGCGAAGGTTGGAATGGGAAATCACCCCCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc                           <  1:260018/238‑1 (MQ=255)
                                                                                                                                                                                                          atttCCAGGGAAGCGCAGGTTGAAATGGGTAATCACCCCCCCTTTCGGTCGGGAGATGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGt                        <  1:436227/235‑1 (MQ=255)
                                                                                                                                                                                                             cccAGCGAAGCTCAGGTGGTAAGTGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGc                                                          <  1:190951/201‑1 (MQ=255)
                                                                                                                                                                                                                ggggAAGCGGAGGTTGTAAGGGGTAATCCCCCCCCCTTGCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:339987/236‑1 (MQ=255)
                                                                                                                                                                                                                ggggAAGCGCAGGTTGTAATGGGTAATCACCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:452554/236‑1 (MQ=255)
                                                                                                                                                                                                                gcgggggggcAGGTTGTAATGGGTAATGCCCCCCCCTTTCGGTGGGGAGGGGGTGGCGGGGGTGAAGAAAATTTCCGCCGTTTCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  2:873661/231‑1 (MQ=255)
                                                                                                                                                                                                                cgggaggcggAGTTTGTAATGGGTAATCCCCCGCCCTTGCGGGCGGGGGGTGGGTCCGGGGGTGAAGAAAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  2:410456/233‑1 (MQ=255)
                                                                                                                                                                                                                cgcggaggggaGGTTGAAAGGGGTAATCCCCACCCTTTTGGGTCGGGGGGTGGGGCGGGGGGTAAAGGGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  2:447256/229‑1 (MQ=255)
                                                                                                                                                                                                                aGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGGAAACGAACTCACGCCATCATCAGCGGGAAGTGCCCCATCTTTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTTGTACa                 >  1:400776/1‑239 (MQ=255)
                                                                                                                                                                                                                aGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGGGGATAGCCCAAGGGGGGCGGGTGTTGAGTTTTTAGTTGAGTAAACGAACTTCCGCCAACAACCACCGGGAGGTGCCCCTCTGTCCGGGCAAAGTGGCGCCAATTTGTGGCCTCTCCCTTCCGAATCTGTTgtataa                 >  2:292120/1‑234 (MQ=255)
                                                                                                                                                                                                                   gaaggggggtttGTAAGGGGTAATCACCCCCCCTTTCGGGCGGGAGGGGGTGCCGGGGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  2:150734/226‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTTTAATGGGTAACCCCCCCCCCTTTCGGTCGGGAGGGTGTGCCGGAGGTTAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:368109/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAAGTGAAGAGAATTTGCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:1267/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTGTAATGGGTAAGCACCCCACCTTTGGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:572510/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTGTAATGGGTAACCTCCACCCCGTTCGGGCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:351381/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCGCAGGTTGAAATGGGTAATCCCCACCCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:52007/236‑1 (MQ=255)
                                                                                                                                                                                                                   gAAGCCCAGGTTGAAATGGGTAATCACCACCCCTTTGGGGCGGGAGGGGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACa                 <  1:113281/236‑1 (MQ=255)
                                                                                                                                                                                                                      gccgagttgttAAGGGGTAATCCCCCCTCCTGGCGGGCGGGAGGTGGGGCCGGGGGTGAGGAGAATTGCCGCCGTAACGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGa            <  2:188641/228‑1 (MQ=255)
                                                                                                                                                                                                                       cgcAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGGTGTTTATTTTTCAGTTGGGTAAACGAACTCACGCCATCATCAGCGGGAAGGGCCCCATTTTTCAGGCAAATGTGCCCCAATTGAGTGGCATCTTCCTGCTTAATCTTTTGATAACTAGGATa          >  1:447256/1‑239 (MQ=255)
                                                                                                                                                                                                                        ggaGGTTGAAATGGGTAAACCCCACCCCTTTCGGTCGGGAGGGGGGCCCGGAGGTGAAGAGAATTGCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATa          <  1:451267/236‑1 (MQ=255)
                                                                                                                                                                                                                          ggtttttAATGGGTAATCCCCACACCTTGCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGaa       <  2:289945/233‑1 (MQ=255)
                                                                                                                                                                                                                             gtttaagtggTAAGCCCCACCCCTTTCGGGCGGGAGGTGGGGCCGGGGGTGAAGGGAATTTCCGCCGGATCGGCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAAtt     <  1:321965/230‑1 (MQ=255)
                                                                                                                                                                                                                               gTAATGGGAAATCACCACACCTTTCGGTCGGGAGGTGGGGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTgcg  <  1:63613/239‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                    
CGGAATACATTCGGTAACGGTGGCGCGGTACTTCTGTACCTGTCCCCAGAAGGCGCGGGCGCTGTATTTCTCGACCAGCACAAAGGTGGCCCCGGCAGAAAACGCCGCCATCGCCGCAGTACACTGGCAATCGATATGAAACGCAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGTATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCG  >  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: