Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,737,956 2 bp→CG coding (623‑624/1254 nt) avtA → alanine‑alpha‑ketoisovalerate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,737,9560GC100.0% 692.9 / NA 220G208A (GGC→GCG) avtAalanine‑alpha‑ketoisovalerate transaminase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (109/111);  total (109/111)
*NC_0009133,737,9570CG99.6% 695.2 / NA 220G208A (GGC→GCGavtAalanine‑alpha‑ketoisovalerate transaminase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (109/110);  total (109/111)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 231 total aligned reads displayed.

AAACCGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGGCAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGCCCGCTATG  >  NC_000913/3737869‑3738043
                                                                                       ||                                                                                      
aaaCCGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATc                                                                                    >  2:859700/1‑93 (MQ=255)
  aCCGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGCTGCTGAAGCTTGACGCGCTGGCGAATCaa                                                                                  >  2:1808441/1‑93 (MQ=255)
   ccGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                                                 >  2:4522966/1‑93 (MQ=255)
   ccGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                                                 <  1:1772380/93‑1 (MQ=255)
   ccGGGATGATTTGCGCCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                                                 <  1:4645428/93‑1 (MQ=255)
   acGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAac                                                                                 <  1:159573/92‑1 (MQ=255)
    cGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAaca                                                                                <  2:530607/93‑1 (MQ=255)
    cGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAaca                                                                                <  1:1565658/93‑1 (MQ=255)
     gggATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacat                                                                               >  1:4542211/1‑92 (MQ=255)
     gggATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacac                                                                               >  2:103808/1‑93 (MQ=255)
     gggATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacac                                                                               <  2:2942999/93‑1 (MQ=255)
     gggATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAacac                                                                               <  1:4372775/93‑1 (MQ=255)
     gggATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGAGCTGGCGAATCAacag                                                                               >  1:993035/1‑92 (MQ=255)
      ggATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACg                                                                              >  2:752556/1‑93 (MQ=255)
      ggATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACg                                                                              >  1:3452435/1‑93 (MQ=255)
        atgatTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGc                                                                            <  2:4316514/93‑1 (MQ=255)
         tgatTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCa                                                                           >  2:1146564/1‑93 (MQ=255)
         tgatTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCa                                                                           <  1:2487088/93‑1 (MQ=255)
          gatTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCAt                                                                          >  1:3238869/1‑93 (MQ=255)
           atTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCAtt                                                                         <  2:2639019/93‑1 (MQ=255)
            tttGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTc                                                                        <  2:3999876/93‑1 (MQ=255)
             ttGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTcc                                                                       <  1:417998/93‑1 (MQ=255)
             ttGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTcc                                                                       <  2:2185363/93‑1 (MQ=255)
               gCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGc                                                                     <  2:2577143/93‑1 (MQ=255)
                  tctcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTgg                                                                  <  2:3668095/93‑1 (MQ=255)
                   ctcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTggg                                                                 >  2:4474938/1‑92 (MQ=255)
                    tcCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtg                                                                >  1:709347/1‑93 (MQ=255)
                      ccGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgat                                                              <  2:4148787/93‑1 (MQ=255)
                        ggCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATTAACACGGCATTCCGCTGGtgattg                                                            >  1:2671317/1‑93 (MQ=255)
                         gCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattga                                                           <  1:4074874/93‑1 (MQ=255)
                          ccGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGtgattgat                                                          >  2:3383221/1‑93 (MQ=255)
                           cgacgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATa                                                         <  1:1130176/93‑1 (MQ=255)
                            gacgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTTGCGAATCAACACGGCATTCCGCTGGGGATTGATaa                                                        >  1:2345923/1‑93 (MQ=255)
                            gacgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATaa                                                        <  1:433936/93‑1 (MQ=255)
                              cgaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACg                                                      >  1:3798435/1‑93 (MQ=255)
                               gaATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGc                                                     >  2:1956256/1‑93 (MQ=255)
                                aaTCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCt                                                    <  1:4481975/93‑1 (MQ=255)
                                  tCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTa                                                  >  2:1570961/1‑93 (MQ=255)
                                  tCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTa                                                  <  2:2650615/93‑1 (MQ=255)
                                  tCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTa                                                  >  1:4093152/1‑93 (MQ=255)
                                  tCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTa                                                  <  1:2440918/93‑1 (MQ=255)
                                   ccAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTAt                                                 <  2:3141213/93‑1 (MQ=255)
                                   ccAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTAt                                                 <  1:2143614/93‑1 (MQ=255)
                                    cAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATg                                                <  2:2166621/93‑1 (MQ=255)
                                     aaCAGGCAATGTGATTACTGACTAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATgg                                               <  1:1224641/93‑1 (MQ=255)
                                     aaCAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATgg                                               >  1:3064079/1‑93 (MQ=255)
                                     aaCAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATgg                                               <  2:2476962/93‑1 (MQ=255)
                                      aCAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGc                                              <  1:2716873/93‑1 (MQ=255)
                                          gCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTcc                                          >  1:114992/1‑93 (MQ=255)
                                           cAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTccc                                         >  2:4646937/1‑93 (MQ=255)
                                            aaTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCg                                        <  2:4437377/93‑1 (MQ=255)
                                             aTGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGt                                       >  1:475687/1‑93 (MQ=255)
                                               gtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                                     >  2:3954060/1‑93 (MQ=255)
                                               gtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                                     >  2:4058240/1‑93 (MQ=255)
                                               gtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                                     >  2:1487760/1‑93 (MQ=255)
                                               gtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                                     >  2:1072865/1‑93 (MQ=255)
                                               gtgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTc                                     >  1:1692319/1‑93 (MQ=255)
                                                tgATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTcc                                    >  2:327532/1‑93 (MQ=255)
                                                 gATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTccc                                   >  1:3488661/1‑93 (MQ=255)
                                                   ttACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCgg                                 <  1:1820158/93‑1 (MQ=255)
                                                   ttACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCgg                                 >  2:3931639/1‑93 (MQ=255)
                                                   ttACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCgg                                 <  2:2763622/93‑1 (MQ=255)
                                                     aCTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGt                               >  2:1223300/1‑93 (MQ=255)
                                                      cTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTa                              <  1:1496797/93‑1 (MQ=255)
                                                      cTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTa                              >  1:2869656/1‑93 (MQ=255)
                                                        gACGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatc                            <  2:3420772/93‑1 (MQ=255)
                                                         aCGAAGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTatca                           <  2:3619630/93‑1 (MQ=255)
                                                            aaGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCt                        >  2:4163804/1‑93 (MQ=255)
                                                            aaGAGTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGCGTATCATCt                        >  2:3785387/1‑93 (MQ=255)
                                                              gagTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCCTTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTc                      >  1:1030095/1‑93 (MQ=255)
                                                              gagTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTc                      >  2:1293471/1‑93 (MQ=255)
                                                                gTTGCTGAAGCTTGACGCGCTGGCGAATCAACTCGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAg                    <  2:3875703/93‑1 (MQ=255)
                                                                gTTGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAg                    <  2:4346562/93‑1 (MQ=255)
                                                                gTTGCTGAAGCTTGACGCGCTGCCGAATCAACACGCCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAg                    <  2:3037674/93‑1 (MQ=255)
                                                                 ttGCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGt                   <  2:1394985/93‑1 (MQ=255)
                                                                   gCTGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGa                 <  2:664848/93‑1 (MQ=255)
                                                                     tGAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGCCGTCCCGTTCCCGGGTATCATCTTCAGTGAAg               >  1:1934362/1‑93 (MQ=255)
                                                                      gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCATTGAAgc              >  2:3522010/1‑93 (MQ=255)
                                                                      gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgc              >  1:307721/1‑93 (MQ=255)
                                                                      gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgc              <  1:1292769/93‑1 (MQ=255)
                                                                      gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgc              >  1:224709/1‑93 (MQ=255)
                                                                      gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgc              >  2:2840680/1‑93 (MQ=255)
                                                                      gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgc              <  1:537478/93‑1 (MQ=255)
                                                                      gAAGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgc              >  1:740224/1‑93 (MQ=255)
                                                                       aaGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgcg             >  1:2989334/1‑93 (MQ=255)
                                                                       aaGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgcg             >  1:3498195/1‑93 (MQ=255)
                                                                        aGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgcgc            >  2:2336438/1‑93 (MQ=255)
                                                                        aGCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgcgc            <  1:1389570/93‑1 (MQ=255)
                                                                         gCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgcgcg           >  1:2560267/1‑93 (MQ=255)
                                                                         gCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgcgcg           >  2:3477710/1‑93 (MQ=255)
                                                                         gCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgcgcg           <  1:3730650/93‑1 (MQ=255)
                                                                         gCTTGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAgcgcg           >  2:3165927/1‑93 (MQ=255)
                                                                           ttGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGcc         <  2:425934/93‑1 (MQ=255)
                                                                           ttGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGcc         <  1:3534119/93‑1 (MQ=255)
                                                                           ttGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGcc         >  1:3887433/1‑93 (MQ=255)
                                                                            tGACGCGCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGccc        >  1:4433722/1‑93 (MQ=255)
                                                                               cgcgCTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGCCCGCt     <  1:2105947/93‑1 (MQ=255)
                                                                                gcgcTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGGTCCCGGGTATCATCTTCAGTGAAGCGCGCCCGCTa    >  1:3271438/1‑93 (MQ=255)
                                                                                ccgcTGGCGAATCAACACGGCATTCCGCTGTTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGCCCGCTa    <  2:4466032/92‑1 (MQ=255)
                                                                                 cgcTGGCGAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGCCCGCTAt   >  1:3268166/1‑93 (MQ=255)
                                                                                       ||                                                                                      
AAACCGGGATGATTTGCGTCTCCCGGCCGACGAATCCAACAGGCAATGTGATTACTGACGAAGAGTTGCTGAAGCTTGACGCGCTGGGCAATCAACACGGCATTCCGCTGGTGATTGATAACGCTTATGGCGTCCCGTTCCCGGGTATCATCTTCAGTGAAGCGCGCCCGCTATG  >  NC_000913/3737869‑3738043

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