Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 39,079 C→T A13T (GCC→ACC)  caiB ← crotonobetainyl CoA:carnitine CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091339,0790CT98.5% 252.3 / NA 66A13T (GCC→ACC) caiBcrotonobetainyl CoA:carnitine CoA transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/40);  minor base G (0/1);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTC  >  NC_000913/38931‑39223
                                                                                                                                                    |                                                                                                                                                
gagaGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCa                                                                                                                                                <  2:452629/151‑1 (MQ=255)
     ttGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACg                                                                                                                                             >  2:169895/1‑149 (MQ=255)
     ttGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGc                                                                                                                                           >  2:251212/1‑151 (MQ=255)
      tGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGc                                                                                                                                           <  2:714716/150‑1 (MQ=255)
      tGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGc                                                                                                                                           <  1:1000538/150‑1 (MQ=255)
       gCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGcc                                                                                                                                          <  2:462337/150‑1 (MQ=255)
                   gTTGAACGCGAATGGGGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCa                                                                                                                              >  2:430835/1‑150 (MQ=255)
                        aCGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTa                                                                                                                        <  2:174591/151‑1 (MQ=255)
                          gcgAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGt                                                                                                                         >  1:1059884/1‑148 (MQ=255)
                           cgAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGGCAACGGCCCGAATTTCGGCATGGGTAgat                                                                                                                     <  2:988028/151‑1 (MQ=255)
                             aaTGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAgatga                                                                                                                   <  2:144123/151‑1 (MQ=255)
                                ggTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATcc                                                                                                                <  1:709621/151‑1 (MQ=255)
                                ggTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATcc                                                                                                                >  1:863542/1‑151 (MQ=255)
                                ggTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATcc                                                                                                                >  2:289483/1‑151 (MQ=255)
                                 gtgtCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCa                                                                                                               <  2:664909/151‑1 (MQ=255)
                                       gCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTAttt                                                                                                         <  2:429237/151‑1 (MQ=255)
                                        cccAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTg                                                                                                        <  1:118619/151‑1 (MQ=255)
                                             gCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCt                                                                                                   <  1:383854/151‑1 (MQ=255)
                                             gCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCt                                                                                                   <  2:828692/151‑1 (MQ=255)
                                             gCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCt                                                                                                   >  1:770924/1‑151 (MQ=255)
                                                aCGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGa                                                                                                 >  1:874603/1‑150 (MQ=255)
                                                 cGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGaaa                                                                                               <  2:791731/151‑1 (MQ=255)
                                                 cGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGaaa                                                                                               <  2:289780/151‑1 (MQ=255)
                                                  gTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGaaa                                                                                               >  1:68423/1‑150 (MQ=255)
                                                  gTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGaaa                                                                                               >  2:427432/1‑150 (MQ=255)
                                                       cGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCTGAGAAGACAACGCGCAATCCTGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAAt                                                                                            <  2:207926/148‑1 (MQ=255)
                                                       cGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAAtt                                                                                           <  2:127381/149‑1 (MQ=255)
                                                         aTCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATg                                                                                       <  1:424912/151‑1 (MQ=255)
                                                                tAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAgcgc                                                                                 >  1:315902/1‑150 (MQ=255)
                                                                 aaCTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCa                                                                                <  1:1048816/150‑1 (MQ=255)
                                                                       cgcgCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGc                                                                         >  2:471066/1‑151 (MQ=255)
                                                                        gcgcCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGcc                                                                        <  1:121881/151‑1 (MQ=255)
                                                                         cgcCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCg                                                                       >  2:933532/1‑151 (MQ=255)
                                                                                   tGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGt                                                             >  1:676563/1‑151 (MQ=255)
                                                                                    gCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATgcggc                                                               <  1:765793/148‑1 (MQ=255)
                                                                                    gCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGt                                                             <  2:455768/150‑1 (MQ=255)
                                                                                    gCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTa                                                            <  1:291323/151‑1 (MQ=255)
                                                                                       aaCATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCTAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAAc                                                         <  2:614623/151‑1 (MQ=255)
                                                                                                 cGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCg                                               >  1:1007438/1‑151 (MQ=255)
                                                                                                 cGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCg                                               <  1:544505/151‑1 (MQ=255)
                                                                                                   gCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGc                                             <  1:261610/151‑1 (MQ=255)
                                                                                                     aaaCGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCt                                           >  2:258256/1‑151 (MQ=255)
                                                                                                     aaaCGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCt                                           >  1:31776/1‑151 (MQ=255)
                                                                                                      aaCGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCt                                           >  2:710582/1‑150 (MQ=255)
                                                                                                           tCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTc                                     >  1:391323/1‑151 (MQ=255)
                                                                                                            ccGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTc                                     <  2:179640/150‑1 (MQ=255)
                                                                                                            ccGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTc                                     <  1:660377/150‑1 (MQ=255)
                                                                                                             cGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGt                                   <  2:622350/151‑1 (MQ=255)
                                                                                                             cGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGt                                   >  2:469989/1‑151 (MQ=255)
                                                                                                             cGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGt                                   <  2:511140/151‑1 (MQ=255)
                                                                                                                 gATTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTc                               <  1:1028143/151‑1 (MQ=255)
                                                                                                                  aTTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTcc                              <  1:12022/151‑1 (MQ=255)
                                                                                                                  aTTTCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTcc                              >  1:465869/1‑151 (MQ=255)
                                                                                                                    ttCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTccc                             <  1:1006010/150‑1 (MQ=255)
                                                                                                                    ttCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCg                            >  1:504748/1‑151 (MQ=255)
                                                                                                                     tCGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCg                            <  2:165718/150‑1 (MQ=255)
                                                                                                                      cGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGt                           <  1:1007151/150‑1 (MQ=255)
                                                                                                                      cGATACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGt                           <  1:251221/150‑1 (MQ=255)
                                                                                                                        aTACCGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAgg                        <  1:788825/151‑1 (MQ=255)
                                                                                                                            cGGAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTg                    >  2:415464/1‑151 (MQ=255)
                                                                                                                             ggAGAAGACAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGa                   <  2:589769/151‑1 (MQ=255)
                                                                                                                                 aagaCAACGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATa                 >  1:829781/1‑149 (MQ=255)
                                                                                                                                       aCGCGCAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCAt         >  2:121941/1‑151 (MQ=255)
                                                                                                                                         gcgcAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCa       <  1:652782/151‑1 (MQ=255)
                                                                                                                                           gcAATCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGc     <  1:15859/151‑1 (MQ=255)
                                                                                                                                              aTCCGGTCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTc  <  1:47186/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GAGAGTTGCGGGTAGTTCGGTTGAACGCGAATGGTGTCGGCCCAGGCGACGTTCTCGATCCAGATAACTTCCGCGCCCCATTCTGCGAACATTTGCCCGGCAAACGGTCCGGCGATTTCGATACCGGAGAAGACAACGCGCAATCCGGCCAACGGCCCGAATTTCGGCATGGGTAGATGATCCATTATTTGCTCCTGAAAAATTTATGTAGCGCATGACTGCCGGATGCGGCGTAAACGCTTTATCCGGCCTACATTCGTGCTCCCGTAGGCCTGATAAGACGCATCAGCGTC  >  NC_000913/38931‑39223

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