Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,454 T→A 7.5% Q124L (CAA→CTA)  caiC ← putative crotonobetaine/carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091337,4540TA 7.5% 196.8 / 4.1 64Q124L (CAA→CTA) caiCputative crotonobetaine/carnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (24/35);  new base A (3/2);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGA  >  NC_000913/37306‑37598
                                                                                                                                                    |                                                                                                                                                
gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt                                                                                                                                                  >  1:285802/1‑149 (MQ=255)
gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt                                                                                                                                                  <  3:91962/149‑1 (MQ=255)
 tCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt                                                                                                                                                  >  8:84932/1‑148 (MQ=255)
 tCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt                                                                                                                                                  >  4:421317/1‑148 (MQ=255)
   aGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGa                                                                                                                                                <  8:363623/148‑1 (MQ=255)
   aGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGa                                                                                                                                                >  7:153263/1‑148 (MQ=255)
   aGTCAATAGCGGCGGTGCATAGCACAAGGTGGCAGGGTGTTGATCTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGAt                                                                                                                                               <  2:82605/149‑1 (MQ=255)
    gTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATa                                                                                                                                              <  2:285802/149‑1 (MQ=255)
       gATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACAt                                                                                                                                           <  6:142522/149‑1 (MQ=255)
        gatacggcggTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATa                                                                                                                                          <  2:353264/145‑1 (MQ=255)
         tAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAg                                                                                                                                         <  5:438058/149‑1 (MQ=255)
            cggcggTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGAt                                                                                                                                      <  1:24914/149‑1 (MQ=255)
             ggcggTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATa                                                                                                                                     <  1:410992/149‑1 (MQ=255)
             ggcggTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGATGAGTAACCGCACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAAATGAGTGCCATCTTCCTGCTGAATCAGTTGATACATAGGa                                                                                                                                       >  1:414561/1‑147 (MQ=255)
                 gTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAAt                                                                                                                                 <  8:420921/149‑1 (MQ=255)
                            aaGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGt                                                                                                                      >  2:55097/1‑149 (MQ=255)
                                    aGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACcagca                                                                                                              <  7:460332/149‑1 (MQ=255)
                                      gttgttGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAgg                                                                                                            >  7:249864/1‑149 (MQ=255)
                                      gttgttGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAgg                                                                                                            <  7:84932/149‑1 (MQ=255)
                                          ttgATTTTTCAGTTGGGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACg                                                                                                        >  1:455581/1‑149 (MQ=255)
                                          ttgATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACg                                                                                                        >  1:451641/1‑149 (MQ=255)
                                          ttgATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCAGTCAGGCAAATGTGCCGCAATTGAGTGGCATATTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACg                                                                                                        >  8:207586/1‑149 (MQ=255)
                                            gATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGATCACCAGCAGGCACGc                                                                                                       <  7:204578/148‑1 (MQ=255)
                                                  tCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCt                                                                                                 <  7:207586/148‑1 (MQ=255)
                                                         ggTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCa                                                                                         <  2:455581/148‑1 (MQ=255)
                                                         aGTAAACGAACTCACGCCATCATCAGCTGGAAGTGCCACATATGTCAGGCAAATGTGACAAAATTGATTTGCATCTTCCTGCGGAATATTGAGATAAATAGGATAGAATTGAGCACTTGTCACCAGAAGGAACGACTGGCTATTTTg                                                                                           >  7:323940/1‑147 (MQ=255)
                                                         aGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAAtt                                                                                                                                <  7:65612/110‑1 (MQ=255)
                                                         aGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAAtt                                                                                                                                >  8:65612/1‑110 (MQ=255)
                                                         aGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCa                                                                                         <  8:153263/149‑1 (MQ=255)
                                                             aaCGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGa                                                                                      >  6:294193/1‑148 (MQ=255)
                                                             aaCGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGAt                                                                                     <  2:907/149‑1 (MQ=255)
                                                             aaCGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGAt                                                                                     <  3:421317/149‑1 (MQ=255)
                                                                gAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCa                                                                                  >  2:278754/1‑149 (MQ=255)
                                                                    tCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCAcgcg                                                                              >  6:36290/1‑149 (MQ=255)
                                                                       cGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTAGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCtt                                                                           >  5:290666/1‑149 (MQ=255)
                                                                       cGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTAGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCtt                                                                           >  5:297068/1‑149 (MQ=255)
                                                                          catcatCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTcc                                                                        <  4:73635/149‑1 (MQ=255)
                                                                           agcatcaGCTTGTAGTGTCTTATCTGTAAGGGAAATTTTTCGCTATTGAGTGGCATCTTACTGCTGAATCTGTTGATACATATTATAGAATTGCGCACTGTTCTCCAGAAGGCTCGCGTTGGTATTTTGCAGGATCCAAGCGGTTTCCt                                                                       <  8:323940/147‑1 (MQ=255)
                                                                            tcatcaGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTc                                                                      <  5:36290/149‑1 (MQ=255)
                                                                            tcatcaGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTc                                                                      <  6:445145/149‑1 (MQ=255)
                                                                              atcaGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCg                                                                     <  4:375361/148‑1 (MQ=255)
                                                                                        aGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcg                                                          >  5:15076/1‑149 (MQ=255)
                                                                                        aGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCGGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcg                                                          <  2:402782/149‑1 (MQ=255)
                                                                                                   tGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTc                                                                      <  2:297044/126‑1 (MQ=255)
                                                                                                   tGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTc                                                                      >  1:297044/1‑126 (MQ=255)
                                                                                                    gTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCgg                                              >  4:82387/1‑149 (MQ=255)
                                                                                                    gTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCgg                                              >  1:11555/1‑149 (MQ=255)
                                                                                                     tCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATc                                                <  8:249864/146‑1 (MQ=255)
                                                                                                     tCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATc                                                <  2:451641/146‑1 (MQ=255)
                                                                                                     tCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGc                                             <  1:55097/149‑1 (MQ=255)
                                                                                                        ggCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCAc                                           >  6:162211/1‑148 (MQ=255)
                                                                                                         gCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCCCGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCa                                         >  2:208926/1‑149 (MQ=255)
                                                                                                          cAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCa                                         >  2:30911/1‑148 (MQ=255)
                                                                                                                  ccGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCgcgc                                <  6:15076/149‑1 (MQ=255)
                                                                                                                      aaTTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcgggc                                                       <  2:304239/122‑1 (MQ=255)
                                                                                                                      aaTTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAAt                            >  7:436609/1‑149 (MQ=255)
                                                                                                                           aggggCATCTTACTGCTGAATCTGTTTATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACCCCTGGCTATTTTGGAGGATCCACGCGCTTTCCTCGCACCCCAGGCGGGCGTTAATAGGCACCATAATCGCGCCAAttttt                        <  1:306187/145‑1 (MQ=255)
                                                                                                                           aGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAAttttt                        <  1:278754/148‑1 (MQ=255)
                                                                                                                              ggCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGc                      >  5:180952/1‑147 (MQ=255)
                                                                                                                                 aTCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcgg                                                         <  8:338511/109‑1 (MQ=255)
                                                                                                                                 aTCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAggcgg                                                         >  7:338511/1‑109 (MQ=255)
                                                                                                                                 aTCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGACACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGcc                 >  7:13531/1‑149 (MQ=255)
                                                                                                                                        tgctgAATCTGTTGATACATAGGATAGAAATGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTACTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAAc            >  5:41195/1‑147 (MQ=255)
                                                                                                                                          ctgAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCa          <  5:294193/147‑1 (MQ=255)
                                                                                                                                           tgAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCa                                                                                         >  8:104621/1‑67 (MQ=255)
                                                                                                                                           tgAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCa                                                                                         <  7:104621/67‑1 (MQ=255)
                                                                                                                                              aTCTGTAGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCaaaa    <  6:290666/149‑1 (MQ=255)
                                                                                                                                              aTCTGTAGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCaaaa    <  6:297068/149‑1 (MQ=255)
                                                                                                                                                   ttGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGa  <  5:162211/146‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATACATAGGATAGAATTGCGCACTGGTCACCAGCAGGCACGCCTGGCTATTTTGCAGGATCCACGCGCTTTCCTCGCACAACAGGCGGGCGTTAATCGGCACCATAATCGCGCCAATTTTTGCCAGCCCGAACCAGCAAAAGA  >  NC_000913/37306‑37598

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: