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,8820CT99.8% 1826.3 / ‑6.4 501T331T (ACG→ACAargIornithine carbamoyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (229/271);  total (230/271)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

Only 100 of 512 total aligned reads displayed.

TGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGT  >  NC_000913/4474795‑4474971
                                                                                       |                                                                                         
atagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTcgc                                                                                      >  1:10016482‑M2/75‑93 (MQ=255)
 tagacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTcgcg                                                                                     >  1:7430918‑M2/74‑93 (MQ=255)
  agacctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGa                                                                                    <  2:7113964‑M2/21‑1 (MQ=255)
    acctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGAcc                                                                                  >  1:4937730‑M2/71‑93 (MQ=255)
    acctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGAcc                                                                                  <  1:4292194‑M2/23‑1 (MQ=255)
     cctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCa                                                                                 >  2:457917‑M2/70‑93 (MQ=255)
      ctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCAt                                                                                <  2:1787952‑M2/25‑1 (MQ=255)
      ctatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCAt                                                                                <  2:9118257‑M2/25‑1 (MQ=255)
       tatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATc                                                                               >  1:1460614‑M2/68‑93 (MQ=255)
       tatcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATc                                                                               <  2:5785833‑M2/26‑1 (MQ=255)
        atcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCa                                                                              <  1:6099790‑M2/27‑1 (MQ=255)
         tcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCAc                                                                             <  2:1326482‑M2/28‑1 (MQ=255)
         tcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCAc                                                                             <  1:7906277‑M2/28‑1 (MQ=255)
         tcactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCAc                                                                             >  2:4705634‑M2/66‑93 (MQ=255)
          cactgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCAcc                                                                            >  2:3473653‑M2/65‑93 (MQ=255)
           actgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCg                                                                           >  2:1180547‑M2/64‑93 (MQ=255)
            ctgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGc                                                                          >  2:9188724‑M2/63‑93 (MQ=255)
             tgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCt                                                                         >  1:3989134‑M2/62‑93 (MQ=255)
             tgcataaaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCt                                                                         >  2:9754818‑M2/62‑93 (MQ=255)
                   aaataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATa                                                                   <  1:1315477‑M2/38‑1 (MQ=255)
                   aaataatgtgatcccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATa                                                                   <  2:2102444‑M2/38‑1 (MQ=255)
                    aataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAg                                                                  <  1:7080249‑M2/39‑1 (MQ=255)
                    aataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAg                                                                  >  2:5388691‑M2/55‑93 (MQ=255)
                    aataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAg                                                                  <  1:8030309‑M2/39‑1 (MQ=255)
                     ataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGt                                                                 <  2:8886656‑M2/40‑1 (MQ=255)
                     ataatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGt                                                                 >  2:2360930‑M2/54‑93 (MQ=255)
                      taatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTa                                                                <  2:4180285‑M2/41‑1 (MQ=255)
                      taatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTa                                                                <  1:9345829‑M2/41‑1 (MQ=255)
                       aatgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTAt                                                               <  1:8657057‑M2/42‑1 (MQ=255)
                        atgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATg                                                              >  2:1260233‑M2/51‑93 (MQ=255)
                         tgtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGc                                                             >  2:5246260‑M2/50‑93 (MQ=255)
                          gtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCa                                                            >  2:10014987‑M2/49‑93 (MQ=255)
                          gtgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCa                                                            <  2:790006‑M2/45‑1 (MQ=255)
                           tgatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCAt                                                           <  1:6333992‑M2/46‑1 (MQ=255)
                            gatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGATTGTCGCGACCATCACCGCTTTGATAGTATGCATa                                                          <  2:6376787‑M2/47‑1 (MQ=255)
                            gatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATa                                                          <  1:6886228‑M2/47‑1 (MQ=255)
                            gatgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATa                                                          <  2:6610908‑M2/47‑1 (MQ=255)
                              tgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACg                                                        <  1:3697182‑M2/49‑1 (MQ=255)
                              tgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACg                                                        <  1:962992‑M2/49‑1 (MQ=255)
                              tgccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACg                                                        <  1:8159830‑M2/49‑1 (MQ=255)
                              agccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACg                                                        <  2:2454470‑M2/49‑1 (MQ=255)
                               gccgggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACgg                                                       <  1:9903524‑M2/50‑1 (MQ=255)
                                  gggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGttt                                                    >  1:1946214‑M2/41‑93 (MQ=255)
                                  gggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGttt                                                    >  1:438332‑M2/41‑93 (MQ=255)
                                  aggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGttt                                                    >  1:7521130‑M2/41‑93 (MQ=255)
                                   ggatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGtttt                                                   >  1:8230038‑M2/40‑93 (MQ=255)
                                    gatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTc                                                  <  1:6532437‑M2/55‑1 (MQ=255)
                                    gatggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTc                                                  >  1:2021718‑M2/39‑93 (MQ=255)
                                     atggtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTcc                                                 <  1:5468864‑M2/56‑1 (MQ=255)
                                       tgtttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGc                                               <  2:1603449‑M2/58‑1 (MQ=255)
                                         tttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCt                                             <  1:9730591‑M2/60‑1 (MQ=255)
                                         tttgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCt                                             <  2:5670585‑M2/60‑1 (MQ=255)
                                           tgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGa                                           <  1:8839443‑M2/62‑1 (MQ=255)
                                           tgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGa                                           >  2:3675868‑M2/32‑93 (MQ=255)
                                           tgtatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGa                                           <  2:2825705‑M2/62‑1 (MQ=255)
                                            gtattgtccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGAt                                          <  1:3088139‑M2/63‑1 (MQ=255)
                                             tatttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATc                                         >  2:3303341‑M2/30‑93 (MQ=255)
                                                ttcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaaa                                      <  2:3888137‑M2/67‑1 (MQ=255)
                                                 tcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaaaa                                     <  1:3188279‑M2/68‑1 (MQ=255)
                                                 tcccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaaaa                                     >  1:5615882‑M2/26‑93 (MQ=255)
                                                 tcccggcatctttatagcgatagcaATTATTTACTGAGTGGCGCGACCATCACCGCTTTGATAGTATGCATACGGTGTTCCGCCTGATCaaaa                                     <  2:6049239‑M2/68‑1 (MQ=255)
                                                  cccggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCaaaaa                                    <  2:9815136‑M2/69‑1 (MQ=255)
                                                     ggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACaa                                 <  1:7309180‑M2/72‑1 (MQ=255)
                                                     ggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACaa                                 <  2:3033457‑M2/72‑1 (MQ=255)
                                                     ggcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACaa                                 <  2:9374575‑M2/72‑1 (MQ=255)
                                                      gcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAAt                                >  2:8333250‑M2/21‑93 (MQ=255)
                                                      gcatctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAAt                                >  2:143041‑M2/21‑93 (MQ=255)
                                                        atctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGc                              >  1:5087236‑M2/19‑93 (MQ=255)
                                                        atctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGc                              >  1:5950542‑M2/19‑93 (MQ=255)
                                                        atctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGc                              <  2:3159952‑M2/75‑1 (MQ=255)
                                                         tctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCt                             >  1:2193890‑M2/18‑93 (MQ=255)
                                                         tctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCt                             <  2:5924289‑M2/76‑1 (MQ=255)
                                                          ctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                            <  1:2286061‑M2/77‑1 (MQ=255)
                                                          ctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                            <  1:9002126‑M2/77‑1 (MQ=255)
                                                          ctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                            >  1:359483‑M2/17‑93 (MQ=255)
                                                          ctttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTg                            <  1:5259939‑M2/77‑1 (MQ=255)
                                                           tttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTgg                           <  2:5148124‑M2/78‑1 (MQ=255)
                                                           tttatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTgg                           <  2:8518032‑M2/78‑1 (MQ=255)
                                                             tatagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcg                         >  1:5730633‑M2/14‑93 (MQ=255)
                                                              atagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCtggcgt                        <  2:6164210‑M2/81‑2 (MQ=255)
                                                               tagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcggc                       >  2:2615456‑M2/12‑93 (MQ=255)
                                                               tagcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTggcggc                       >  2:5073725‑M2/12‑93 (MQ=255)
                                                                agcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCa                      <  2:2866702‑M2/83‑1 (MQ=255)
                                                                 gcgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAg                     >  1:18640‑M2/10‑93 (MQ=255)
                                                                  cgatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGa                    >  2:4366931‑M2/9‑93 (MQ=255)
                                                                   gatagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAt                   <  1:2555791‑M2/86‑1 (MQ=255)
                                                                    atagcaATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGAtt                  >  1:287105‑M2/7‑93 (MQ=255)
                                                                        caATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGaa              <  1:6446820‑M2/91‑1 (MQ=255)
                                                                        caATTATTTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGaa              <  2:6943712‑M2/91‑1 (MQ=255)
                                                                          attattTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa            <  1:4722416/93‑1 (MQ=255)
                                                                          attattTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGa            <  1:4770160/93‑1 (MQ=255)
                                                                           ttattTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGAc           >  1:423102/1‑93 (MQ=255)
                                                                           ttattTACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGAc           >  2:8138814/1‑93 (MQ=255)
                                                                               ttACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtca       <  2:7454471/93‑1 (MQ=255)
                                                                                tACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcat      <  2:4024297/93‑1 (MQ=255)
                                                                                tACTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcat      >  1:949517/1‑93 (MQ=255)
                                                                                 aCTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcatc     <  2:5523508/93‑1 (MQ=255)
                                                                                  cTGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCtcatca    <  1:8397163/93‑1 (MQ=255)
                                                                                   tGAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAg   <  1:2963799/93‑1 (MQ=255)
                                                                                    gAGTGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGt  >  1:3289984/1‑93 (MQ=255)
                                                                                       |                                                                                         
TGGGTGGAAGAATAGACCTATCACTGCATAAATAATGTGATGCCGGGATGGTTTCCGCATCTTATAGCGATCAGATTATTTACTGAGCGTCGCGACCATCACCGCTTTGATAGTATGCATACGGTTTTCCGCCTGATCAAAAACAATGCTGGCGGCAGATTCGAAGACCTCATCAGT  >  NC_000913/4474795‑4474971

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