Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,474,882 C→T T331T (ACG→ACA argI ← ornithine carbamoyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,474,8820CT98.0% 685.8 / ‑3.2 202T331T (ACG→ACAargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (82/113);  total (85/117)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 205 total aligned reads displayed.

CTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGT  >  NC_000913/4474792‑4474971
                                                                                          |                                                                                         
agaatagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAgtgt                                                                                         >  2:1445179‑M2/78‑93 (MQ=255)
agaatagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAgtgt                                                                                         >  1:3517189‑M2/78‑93 (MQ=255)
 gaatagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTc                                                                                        <  1:1377050‑M2/17‑1 (MQ=255)
  aatagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTTAGTGTcg                                                                                       <  2:744033‑M2/18‑1 (MQ=255)
   atagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTcgc                                                                                      >  2:863189‑M2/75‑93 (MQ=255)
   atagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTcgc                                                                                      >  1:4513504‑M2/75‑93 (MQ=255)
   atagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTcgc                                                                                      <  2:3800639‑M2/19‑1 (MQ=255)
     agacctatcactgcataaaataatgtgatgccgggattgtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGa                                                                                    <  1:1621536‑M2/21‑1 (MQ=255)
        cctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCa                                                                                 <  2:2052597‑M2/24‑1 (MQ=255)
        cctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCa                                                                                 <  1:3239409‑M2/24‑1 (MQ=255)
        cctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCa                                                                                 >  1:1140704‑M2/70‑93 (MQ=255)
           atcactgcataaaataatgtgatgccgggatggtttgtatttctcggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCa                                                                              >  1:3453603‑M2/67‑93 (MQ=255)
            tcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCAc                                                                             >  1:1175883‑M2/66‑93 (MQ=255)
             cactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCAcc                                                                            >  1:1011401‑M2/65‑93 (MQ=255)
              actgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCg                                                                           >  2:408556‑M2/64‑93 (MQ=255)
               ctgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGc                                                                          <  1:755513‑M2/31‑1 (MQ=255)
                tgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCt                                                                         >  1:1119332‑M2/62‑93 (MQ=255)
                tgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCt                                                                         >  1:885960‑M2/62‑93 (MQ=255)
                tgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCt                                                                         <  1:2395478‑M2/32‑1 (MQ=255)
                tgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCt                                                                         >  2:4344897‑M2/62‑93 (MQ=255)
                  cataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCttt                                                                       >  1:4747234‑M2/60‑93 (MQ=255)
                    taaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGa                                                                     <  1:2669800‑M2/36‑1 (MQ=255)
                    taaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGa                                                                     <  1:2318563‑M2/36‑1 (MQ=255)
                     aaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGAt                                                                    >  1:2048780‑M2/57‑93 (MQ=255)
                      aaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATa                                                                   <  2:3220596‑M2/38‑1 (MQ=255)
                       aataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAg                                                                  <  1:3637671‑M2/39‑1 (MQ=255)
                         taatgtgatgccgggctagtttgtctttcccggcatctttatagcgatagcaATTATTTACTGCGTGTCGCGACCATCACCGCTTTGCTAGTa                                                                <  1:4220928‑M2/41‑1 (MQ=255)
                         taatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTa                                                                <  2:3614466‑M2/41‑1 (MQ=255)
                          aatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTAt                                                               <  1:350233‑M2/42‑1 (MQ=255)
                           atgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATg                                                              <  2:1471859‑M2/43‑1 (MQ=255)
                             gtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCa                                                            <  2:3510473‑M2/45‑1 (MQ=255)
                              tgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCAt                                                           >  1:4735870‑M2/48‑93 (MQ=255)
                              tgatgccgagatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCAt                                                           <  2:4222991‑M2/46‑1 (MQ=255)
                               gatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATa                                                          >  2:4298564‑M2/47‑93 (MQ=255)
                                atgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATAc                                                         <  2:1339054‑M2/48‑1 (MQ=255)
                                 tgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACg                                                        >  2:1463159‑M2/45‑93 (MQ=255)
                                  gccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACgg                                                       <  2:3644368‑M2/50‑1 (MQ=255)
                                  gccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACgg                                                       >  1:1209398‑M2/44‑93 (MQ=255)
                                    cgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATacgggt                                                     >  1:895925‑M2/42‑91 (MQ=255)
                                    cgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGtt                                                     <  1:471118‑M2/52‑1 (MQ=255)
                                     gggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGGTTTGATAGTATGCATACGGgtt                                                    >  2:3439581‑M2/41‑93 (MQ=255)
                                     gggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGttt                                                    <  2:1676367‑M2/53‑1 (MQ=255)
                                     gggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGttt                                                    <  1:687769‑M2/53‑1 (MQ=255)
                                      ggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGtttt                                                   >  1:4324531‑M2/40‑93 (MQ=255)
                                       gatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTc                                                  <  1:1581100‑M2/55‑1 (MQ=255)
                                       gatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTc                                                  <  1:189581‑M2/55‑1 (MQ=255)
                                       gatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTc                                                  <  2:2886568‑M2/55‑1 (MQ=255)
                                        atggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTcc                                                 >  2:3981383‑M2/38‑93 (MQ=255)
                                         tggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCg                                                >  1:173496‑M2/37‑93 (MQ=255)
                                           gtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGcc                                              <  2:2233467‑M2/59‑1 (MQ=255)
                                           gtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGcc                                              <  1:3593278‑M2/59‑1 (MQ=255)
                                            tttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCt                                             <  2:2136436‑M2/60‑1 (MQ=255)
                                            tttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCt                                             <  2:2456083‑M2/60‑1 (MQ=255)
                                              tgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGa                                           >  2:1120055‑M2/32‑93 (MQ=255)
                                                tattgcccggcatctttatagcgatagcaATTATTTACTGAGAGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATc                                         <  1:2971240‑M2/64‑1 (MQ=255)
                                                     cccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaaaaa                                    >  2:468925‑M2/25‑93 (MQ=255)
                                                       cggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACa                                  <  1:2784593‑M2/71‑1 (MQ=255)
                                                        ggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACaa                                 >  2:1494443‑M2/22‑93 (MQ=255)
                                                        ggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACaa                                 <  1:3879048‑M2/72‑1 (MQ=255)
                                                         gcatctttatagcgatagcaGTTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAAt                                <  2:621209‑M2/73‑1 (MQ=255)
                                                          catctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATg                               <  2:2492377‑M2/74‑1 (MQ=255)
                                                          catctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATg                               >  1:3097021‑M2/20‑93 (MQ=255)
                                                          catctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATg                               >  2:1729953‑M2/20‑93 (MQ=255)
                                                           atctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGc                              <  1:4469017‑M2/75‑1 (MQ=255)
                                                           atctttatagcgatagcaATTATTTACTGAGAGTCGCAACCATCACCGCCTTGATAGTATGCATACGGTTTTTCGCCTGATCAAAAACAATGc                              >  1:976059‑M2/19‑93 (MQ=255)
                                                            tctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCt                             <  2:3293711‑M2/76‑1 (MQ=255)
                                                            tctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCt                             <  2:549370‑M2/76‑1 (MQ=255)
                                                            tctttatagcgatagcaATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCt                             <  1:4021067‑M2/76‑1 (MQ=255)
                                                              tttatatcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTgg                           <  1:1484821‑M2/78‑1 (MQ=255)
                                                              tttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTgg                           <  1:1322041‑M2/78‑1 (MQ=255)
                                                              tttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTgg                           >  1:339579‑M2/16‑93 (MQ=255)
                                                              tttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTgg                           >  2:4597584‑M2/16‑93 (MQ=255)
                                                               ttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggc                          <  2:45366‑M2/79‑1 (MQ=255)
                                                                tatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                         <  2:1022312‑M2/80‑1 (MQ=255)
                                                                tatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                         <  1:4123244‑M2/80‑1 (MQ=255)
                                                                tatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                         >  2:3224926‑M2/14‑93 (MQ=255)
                                                                  tagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcggc                       >  1:1692221‑M2/12‑93 (MQ=255)
                                                                  tagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcggc                       >  1:398293‑M2/12‑93 (MQ=255)
                                                                  tagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcggc                       <  2:1944732‑M2/82‑1 (MQ=255)
                                                                    gcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAg                     >  2:214424‑M2/10‑93 (MQ=255)
                                                                      gatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAt                   <  2:1788296‑M2/86‑1 (MQ=255)
                                                                      gatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAt                   <  1:332152‑M2/86‑1 (MQ=255)
                                                                      gatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAt                   >  2:15593‑M2/8‑93 (MQ=255)
                                                                       atagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAtt                  <  2:3364252‑M2/87‑1 (MQ=255)
                                                                       atagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAtt                  >  1:1060185‑M2/7‑93 (MQ=255)
                                                                       atagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAtt                  <  2:2536720‑M2/87‑1 (MQ=255)
                                                                       atagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAtt                  >  1:2277448‑M2/7‑93 (MQ=255)
                                                                        tagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTc                 >  2:3768358‑M2/6‑93 (MQ=255)
                                                                        tagcaATGATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTc                 <  2:1644915‑M2/88‑1 (MQ=255)
                                                                          gcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGa               <  2:4579750‑M2/90‑1 (MQ=255)
                                                                          gcaATTATTTACTGAGCGTCGCGACCTTCACCGCTTTGATAGTATGCATACGGTTTTCAGCCTGATCAAAAACAATGCTGGCGGCAGATTCGa               >  1:382217‑M2/4‑93 (MQ=255)
                                                                             attattTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa            <  1:1958981/93‑1 (MQ=255)
                                                                             attattTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa            <  2:1345921/93‑1 (MQ=255)
                                                                             attattTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa            <  1:3796183/93‑1 (MQ=255)
                                                                              ttattTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGAc           >  1:3136332/1‑93 (MQ=255)
                                                                                   tACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcat      <  2:265232/93‑1 (MQ=255)
                                                                                     cTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcatca    >  2:2948160/1‑93 (MQ=255)
                                                                                      tGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAg   >  1:1771044/1‑93 (MQ=255)
                                                                                      tGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAg   >  2:2331656/1‑93 (MQ=255)
                                                                                      tGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAg   <  2:3197356/93‑1 (MQ=255)
                                                                                          |                                                                                         
CTTTGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGT  >  NC_000913/4474792‑4474971

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: