Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 37,903 C→T D405N (GAC→AAC)  caiB ← crotonobetainyl CoA:carnitine CoA transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091337,9030CT100.0% 239.8 / NA 75D405N (GAC→AAC) caiBcrotonobetainyl CoA:carnitine CoA transferase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (36/39);  total (36/39)

TCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTAT  >  NC_000913/37772‑38031
                                                                                                                                   |                                                                                                                                
tCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           <  1:1344250/139‑1 (MQ=255)
           gTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGTCAGAcc                                                                                                                >  2:355578/1‑139 (MQ=255)
            tCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCt                                                                                                               >  1:45496/1‑139 (MQ=255)
                      gACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCa                                                                                                     >  2:775546/1‑139 (MQ=255)
                       aCGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           >  2:90285/1‑116 (MQ=255)
                       aCGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           <  1:90285/116‑1 (MQ=255)
                          tAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTc                                                                                                 <  2:1190266/139‑1 (MQ=255)
                            gATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCt                                                                                               >  2:977823/1‑139 (MQ=255)
                                  gtCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGt                                                                                         >  1:412292/1‑139 (MQ=255)
                                     cGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           >  1:761370/1‑102 (MQ=255)
                                     cGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           <  2:761370/102‑1 (MQ=255)
                                     cGCCAATGATATCCACTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGt                                                                                      >  2:954549/1‑139 (MQ=255)
                                             atatCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTg                                                                              >  1:576713/1‑139 (MQ=255)
                                                 ccATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGc                                                                          <  2:1280872/139‑1 (MQ=255)
                                                    ttGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTg                                                                                                         <  1:388757/105‑1 (MQ=255)
                                                    ttGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTg                                                                                                         >  2:388757/1‑105 (MQ=255)
                                                     tGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGat                                                                      >  2:771974/1‑139 (MQ=255)
                                                      gCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGata                                                                     <  1:771974/139‑1 (MQ=255)
                                                      gCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGata                                                                     <  1:775546/139‑1 (MQ=255)
                                                       cACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           >  2:667363/1‑84 (MQ=255)
                                                       cACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           <  1:667363/84‑1 (MQ=255)
                                                            ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                                      >  2:1241702/1‑84 (MQ=255)
                                                            ctATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                                      <  1:1241702/84‑1 (MQ=255)
                                                               tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCt                                                                                                                                >  2:638533/1‑71 (MQ=255)
                                                               tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           >  2:984472/1‑76 (MQ=255)
                                                               tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           <  1:984472/76‑1 (MQ=255)
                                                               tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           <  2:869920/76‑1 (MQ=255)
                                                               tCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACt                                                                                                                           >  1:869920/1‑76 (MQ=255)
                                                                   tttttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATgg                                                        <  2:45496/139‑1 (MQ=255)
                                                                   tttttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATgg                                                        <  1:1265529/139‑1 (MQ=255)
                                                                    ttttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCGACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGc                                                       >  2:1426916/1‑139 (MQ=255)
                                                                    ttttGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGc                                                       >  1:1113812/1‑139 (MQ=255)
                                                                       tGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                                    >  2:217058/1‑139 (MQ=255)
                                                                        gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGTCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   <  1:355578/139‑1 (MQ=255)
                                                                        gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   <  1:977823/139‑1 (MQ=255)
                                                                        gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   <  2:938190/139‑1 (MQ=255)
                                                                        gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   <  1:217058/139‑1 (MQ=255)
                                                                        gTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   <  2:1345720/139‑1 (MQ=255)
                                                                         ttCGTTTGTTATTGGGCGGGCGCTAGTCCGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCATCTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                                  >  2:65267/1‑139 (MQ=255)
                                                                         ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCa                                                                                                                    <  1:1053834/73‑1 (MQ=255)
                                                                         ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCa                                                                                                                    >  2:1053834/1‑73 (MQ=255)
                                                                         ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCt                                                                                               <  2:813372/94‑1 (MQ=255)
                                                                         ttCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCt                                                                                               >  1:813372/1‑94 (MQ=255)
                                                                             tttGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTc                                              >  2:634575/1‑139 (MQ=255)
                                                                             tttGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTc                                              <  1:954549/139‑1 (MQ=255)
                                                                               tGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTcca                                            <  1:65267/139‑1 (MQ=255)
                                                                                       ggcgggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCAtgt                                    <  2:801348/139‑1 (MQ=255)
                                                                                        gcgggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGAc                                                                                                       >  1:826008/1‑71 (MQ=255)
                                                                                        gcgggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTACTTCTCAACTTTGGCCAGACCTTTGCTGAc                                                                                                       <  2:826008/71‑1 (MQ=255)
                                                                                           ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGAc                                                                                                                 >  1:765330/1‑58 (MQ=255)
                                                                                           ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGAc                                                                                                                 <  2:765330/58‑1 (MQ=255)
                                                                                           ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   >  1:1362523/1‑120 (MQ=255)
                                                                                           ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   <  2:1362523/120‑1 (MQ=255)
                                                                                           ggcgCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAgg                                >  2:571121/1‑139 (MQ=255)
                                                                                                 aGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTc                                              <  2:178913/119‑1 (MQ=255)
                                                                                                 aGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTc                                              >  1:178913/1‑119 (MQ=255)
                                                                                                 aGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTc                          >  1:284257/1‑139 (MQ=255)
                                                                                                     aGGCAAGCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCcgcg                      <  1:1426916/139‑1 (MQ=255)
                                                                                                          aGCCGACTGACGCCACGCGTTTAGTTCTCAACTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGAt                <  2:1196193/139‑1 (MQ=255)
                                                                                                           gCCGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGAt                <  2:576713/139‑1 (MQ=255)
                                                                                                            ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                                      >  2:140425/1‑36 (MQ=255)
                                                                                                            ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGc                                                                                                                      <  1:140425/36‑1 (MQ=255)
                                                                                                            ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   <  1:390142/103‑1 (MQ=255)
                                                                                                            ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   >  2:390142/1‑103 (MQ=255)
                                                                                                            ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                                    >  2:1353327/1‑102 (MQ=255)
                                                                                                            ccGACTGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGc                                                    <  1:1353327/102‑1 (MQ=255)
                                                                                                                 tGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTc                          <  1:280978/123‑1 (MQ=255)
                                                                                                                 tGACGCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTc                          >  2:280978/1‑123 (MQ=255)
                                                                                                                     gCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGtttt                                                                                   >  2:1224828/1‑62 (MQ=255)
                                                                                                                     gCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGtttt                                                                                   <  1:1224828/62‑1 (MQ=255)
                                                                                                                     gCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                                  >  2:402490/1‑95 (MQ=255)
                                                                                                                     gCCACGCGTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCg                                                  <  1:402490/95‑1 (MQ=255)
                                                                                                                         cgcgTTTAGTTCTCAACTTTGGCCAGGCCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   >  2:943043/1‑90 (MQ=255)
                                                                                                                         cgcgTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGcc                                                   <  1:943043/90‑1 (MQ=255)
                                                                                                                         cgcgTTTAGTTCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTAt  <  1:634575/139‑1 (MQ=255)
                                                                                                                                   |                                                                                                                                
TCCGCAAGATCGTCCCACATTTGACGTAGATGTTGTCCGCCAATGATATCCATTGCACCTCTATCCATTTTTGTTCGTTTGTTATTGGGCGGGCGCTAGTCAGGCAAGCCGACTGACGCCACGCGTTTAGTCCTCAACTTTGGCCAGACCTTTGCTGACCAACTCCTGAATGTCGTTTTCGCTGTAGCCGATATTTTTCAAAATGGCAGCCGTGTCCATGCCATGTGAGGGCATTCCGCGCCAGATTTGTCCGGGGTTAT  >  NC_000913/37772‑38031

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