Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,158,153 C→T 6.6% A435V (GCG→GTG)  argH → argininosuccinate lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,158,1530CT 6.6% 250.2 / 5.3 75A435V (GCG→GTG) argHargininosuccinate lyase
Reads supporting (aligned to +/- strand):  ref base C (34/36);  new base T (3/2);  total (37/38)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.60e-01

TGGTGAAGCGGTGGTGGAAGCCATTCGTCAGGGCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACATTTATATGTATAAATTTGAGCCTGGCTTATCGCCGGGCTTTT  >  NC_000913/4158013‑4158270
                                                                                                                                            |                                                                                                                     
tgGTGAAGCGGTGGTGGAAGCCATTCGTCAGGGCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCgg                                                                                                                     <  4:795967/143‑1 (MQ=255)
tgGTGAAGCGGTGGTGGAAGCCATTCGTCAGGGCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCgg                                                                                                                     >  3:795967/1‑143 (MQ=255)
   tGAAGCGGTGGTGGAAGCCATTCGTCAGGGCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAggc                                                                                                             >  2:743022/1‑148 (MQ=255)
      aGCGGTGGTGGAAGCCATTCGTCAGGGCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAggcggc                                                                                                          >  2:673486/1‑148 (MQ=255)
           tggtggAAGCCATTCGTCAGGGCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGtctc                                                                                                     <  1:673486/148‑1 (MQ=255)
                     cATTCGTCAGGGCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGcagcag                                                                                           <  1:743022/148‑1 (MQ=255)
                              gggCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAg                                                                                  >  3:993722/1‑148 (MQ=255)
                               ggCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCaaa                                                                                                                 >  1:820760/1‑116 (MQ=255)
                               ggCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCaaa                                                                                                                 <  2:820760/116‑1 (MQ=255)
                               ggCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAgg                                                                                 >  1:1096008/1‑148 (MQ=255)
                                gCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGc                                                                                >  4:451013/1‑148 (MQ=255)
                                 cAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCg                                                                               <  3:451013/148‑1 (MQ=255)
                                            aaGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGTGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTgc                                                                    >  4:658556/1‑148 (MQ=255)
                                            aaGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGTGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTgc                                                                    >  2:897896/1‑148 (MQ=255)
                                            aaGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGt                                                                                         >  4:545484/1‑127 (MQ=255)
                                            aaGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGt                                                                                         <  3:545484/127‑1 (MQ=255)
                                            aaGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGc                                                                          <  1:687098/142‑1 (MQ=255)
                                            aaGATCGGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGc                                                                          >  2:687098/1‑142 (MQ=255)
                                               aTCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCa                                                                 <  3:852357/148‑1 (MQ=255)
                                                tCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAg                                                                >  3:803816/1‑148 (MQ=255)
                                                   gCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCg                                                                                                                      <  2:882609/91‑1 (MQ=255)
                                                   gCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCg                                                                                                                      >  1:882609/1‑91 (MQ=255)
                                                        tCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGc                                                                                >  2:1061425/1‑124 (MQ=255)
                                                        tCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGc                                                                                <  1:1061425/124‑1 (MQ=255)
                                                        tCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGtt                                                        <  4:993722/148‑1 (MQ=255)
                                                         cAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCa                                                                                                    <  3:811122/103‑1 (MQ=255)
                                                         cAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCa                                                                                                    >  4:811122/1‑103 (MQ=255)
                                                          aGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAg                                                      <  4:803816/148‑1 (MQ=255)
                                                          aGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAg                                                      <  2:1096008/148‑1 (MQ=255)
                                                              aGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCa                                                                                                    >  3:905388/1‑98 (MQ=255)
                                                              aGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCa                                                                                                    <  4:905388/98‑1 (MQ=255)
                                                                ttGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCa                                                                                                    <  2:630173/96‑1 (MQ=255)
                                                                ttGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCa                                                                                                    >  1:630173/1‑96 (MQ=255)
                                                                   cAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGTGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAAc                                             <  1:897896/148‑1 (MQ=255)
                                                                   cAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGTGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAAc                                             <  3:658556/148‑1 (MQ=255)
                                                                       aaTTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGGAGGCTCGGTTAg                                                      <  1:144933/135‑1 (MQ=255)
                                                                       aaTTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGGAGGCTCGGTTAg                                                      >  2:144933/1‑135 (MQ=255)
                                                                         ttCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCg                                                                                                 <  1:799891/90‑1 (MQ=255)
                                                                         ttCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCg                                                                                                 >  2:799891/1‑90 (MQ=255)
                                                                              tcagGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACATTtata                                      <  2:1041831/144‑1 (MQ=255)
                                                                              tcagGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACATTtata                                      >  1:1041831/1‑144 (MQ=255)
                                                                                  gTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCttt                                                                       <  4:855931/107‑1 (MQ=255)
                                                                                  gTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCttt                                                                       >  3:855931/1‑107 (MQ=255)
                                                                                          cGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTg                                                                           <  3:191645/95‑1 (MQ=255)
                                                                                          cGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTg                                                                           >  4:191645/1‑95 (MQ=255)
                                                                                            aaGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACa                                            <  1:177491/124‑1 (MQ=255)
                                                                                            aaGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACa                                            >  2:177491/1‑124 (MQ=255)
                                                                                              gATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGTGGCAAAAGGCGGCGTCTCACCGCAGCAGGTgg                                                                                       >  3:195315/1‑79 (MQ=255)
                                                                                              gATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTgg                                                                                       <  4:195315/79‑1 (MQ=255)
                                                                                                  tCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAg                                                      <  1:569460/108‑1 (MQ=255)
                                                                                                  tCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAg                                                      >  2:569460/1‑108 (MQ=255)
                                                                                                   cTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCg                                                                               >  1:592891/1‑82 (MQ=255)
                                                                                                   cTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCg                                                                               <  2:592891/82‑1 (MQ=255)
                                                                                                     aTCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCttt                                                                       >  2:817656/1‑88 (MQ=255)
                                                                                                     aTCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCttt                                                                       <  1:817656/88‑1 (MQ=255)
                                                                                                     aTCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCttt                                                                       >  2:864502/1‑88 (MQ=255)
                                                                                                     aTCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCttt                                                                       <  1:864502/88‑1 (MQ=255)
                                                                                                     aTCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACAtt                                          <  3:85624/117‑1 (MQ=255)
                                                                                                     aTCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACAtt                                          >  4:85624/1‑117 (MQ=255)
                                                                                                          aTTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGAt                                                                             >  3:941195/1‑77 (MQ=255)
                                                                                                          aTTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGAt                                                                             <  4:941195/77‑1 (MQ=255)
                                                                                                          aTTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAgg                                                     <  3:215271/101‑1 (MQ=255)
                                                                                                          aTTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAgg                                                     >  4:215271/1‑101 (MQ=255)
                                                                                                          aTTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTa                                                  <  4:991626/104‑1 (MQ=255)
                                                                                                          aTTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTa                                                  >  3:991626/1‑104 (MQ=255)
                                                                                                                        aTCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGc                                                                          <  4:776237/66‑1 (MQ=255)
                                                                                                                        aTCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGc                                                                          >  3:776237/1‑66 (MQ=255)
                                                                                                                            tGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACAttt                                         <  4:603072/95‑1 (MQ=255)
                                                                                                                            tGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACAttt                                         >  3:603072/1‑95 (MQ=255)
                                                                                                                            tGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACATTTATATGTATAAATTTGAGCCTGGCTTATCGCCGGGCtttt  >  4:559985/1‑134 (MQ=255)
                                                                                                                            tGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACATTTATATGTATAAATTTGAGCCTGGCTTATCGCCGGGCtttt  <  3:559985/134‑1 (MQ=255)
                                                                                                                               cTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCttt                                                                       >  3:349511/1‑62 (MQ=255)
                                                                                                                               cTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCttt                                                                       <  4:349511/62‑1 (MQ=255)
                                                                                                                                   aCAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACa                                            >  3:690730/1‑85 (MQ=255)
                                                                                                                                   aCAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACa                                            <  4:690730/85‑1 (MQ=255)
                                                                                                                                            |                                                                                                                     
TGGTGAAGCGGTGGTGGAAGCCATTCGTCAGGGCAAACCGCTGGAAGATCTGCCGCTCAGTGAGTTGCAGAAATTCAGTCAGGTGATTGACGAAGATGTCTATCCGATTCTGTCGCTGCAATCGTGCCTCGACAAGCGTGCGGCAAAAGGCGGCGTCTCACCGCAGCAGGTGGCGCAGGCGATTGCTTTTGCGCAGGCTCGGTTAGGGTAAGAACATTTATATGTATAAATTTGAGCCTGGCTTATCGCCGGGCTTTT  >  NC_000913/4158013‑4158270

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