Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,130,398 A→G D68G (GAC→GGC)  metF → 5,10‑methylenetetrahydrofolate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,130,3980AG99.2% 889.9 / ‑3.7 265D68G (GAC→GGC) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (113/150);  total (113/152)
Fisher's exact test for biased strand distribution p-value = 5.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

Only 100 of 271 total aligned reads displayed.

AGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACG  >  NC_000913/4130308‑4130487
                                                                                          |                                                                                         
aGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGcc                                                                                         >  1:7437116/1‑93 (MQ=255)
 gCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCg                                                                                        >  2:1291472/1‑93 (MQ=255)
  cAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGt                                                                                       >  2:2058129/1‑93 (MQ=255)
   aGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTa                                                                                      <  2:8431263/93‑1 (MQ=255)
   aGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTa                                                                                      <  1:2047311/93‑1 (MQ=255)
      cccTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGc                                                                                   <  1:7237979/93‑1 (MQ=255)
        cTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGcac                                                                                 <  2:7402421/93‑1 (MQ=255)
           tgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGc                                                                              <  1:232492/93‑1 (MQ=255)
           tgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGc                                                                              <  1:2211802/93‑1 (MQ=255)
            ggAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCa                                                                             >  2:8739634/1‑93 (MQ=255)
              aaCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCatt                                                                           <  1:3116308/93‑1 (MQ=255)
               aCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCatta                                                                          >  2:3247301/1‑93 (MQ=255)
                cTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCattat                                                                         <  1:826973/93‑1 (MQ=255)
                cTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCattat                                                                         <  1:4907840/93‑1 (MQ=255)
                 tCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCattatt                                                                        >  2:6391069/1‑93 (MQ=255)
                  ccATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCAttatta                                                                       >  1:8715094/1‑93 (MQ=255)
                    atcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTaaa                                                                     >  1:3877890/1‑93 (MQ=255)
                    atcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTaaa                                                                     <  2:6876284/93‑1 (MQ=255)
                    agcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTaaa                                                                     <  1:7044635/91‑1 (MQ=255)
                     tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                    <  1:4293193/93‑1 (MQ=255)
                     tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                    >  1:4256913/1‑93 (MQ=255)
                     tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                    <  2:2939289/93‑1 (MQ=255)
                      cgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAgg                                                                   <  2:8192087/93‑1 (MQ=255)
                        atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                 >  2:5005713/1‑93 (MQ=255)
                        atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                 <  1:7778515/93‑1 (MQ=255)
                        atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                 <  1:2119941/93‑1 (MQ=255)
                        atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                 >  1:5221566/1‑93 (MQ=255)
                        atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                 >  2:675275/1‑93 (MQ=255)
                         tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAt                                                                <  1:3561762/93‑1 (MQ=255)
                         tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAt                                                                >  2:8341692/1‑93 (MQ=255)
                         tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAt                                                                >  1:6491348/1‑93 (MQ=255)
                          cgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAtt                                                               <  1:6820201/93‑1 (MQ=255)
                           gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGATCTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTa                                                              <  2:1339150/93‑1 (MQ=255)
                           gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTa                                                              <  2:1964590/93‑1 (MQ=255)
                           gCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTa                                                              >  1:8582754/1‑93 (MQ=255)
                            ccTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaa                                                             >  1:5363394/1‑93 (MQ=255)
                            ccTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaa                                                             >  1:4539992/1‑93 (MQ=255)
                             cTTAGCAGCCTGAAACCGAAGTTTGTATGGGTGGCCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaaa                                                            <  2:5562453/93‑1 (MQ=255)
                             cTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaaa                                                            <  2:7625992/93‑1 (MQ=255)
                             cTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaaa                                                            <  1:3330737/93‑1 (MQ=255)
                              ttAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                           <  2:8751401/93‑1 (MQ=255)
                              ttAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                           >  1:1079726/1‑93 (MQ=255)
                              ttAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                           >  1:491764/1‑93 (MQ=255)
                              ttAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                           <  1:5529197/93‑1 (MQ=255)
                               tAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGa                                                          <  2:7132734/93‑1 (MQ=255)
                               tAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGa                                                          <  2:5864075/93‑1 (MQ=255)
                                agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGAt                                                         <  2:2614530/93‑1 (MQ=255)
                                agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGAt                                                         <  2:4309559/93‑1 (MQ=255)
                                agcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGAt                                                         >  2:1104329/1‑93 (MQ=255)
                                 gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGTATTAAAGATc                                                        >  1:7982418/1‑93 (MQ=255)
                                  cagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCg                                                       <  2:3452229/93‑1 (MQ=255)
                                  cagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCg                                                       <  1:5586251/93‑1 (MQ=255)
                                    gcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCa                                                     >  2:2720001/1‑93 (MQ=255)
                                       tGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTg                                                  >  2:620966/1‑93 (MQ=255)
                                       tGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTg                                                  >  2:6088674/1‑93 (MQ=255)
                                        gAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTgg                                                 >  2:7174674/1‑93 (MQ=255)
                                        gAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTgg                                                 <  1:62833/93‑1 (MQ=255)
                                          aaCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTc                                               >  1:8370630/1‑93 (MQ=255)
                                          aaCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTc                                               >  2:4482922/1‑93 (MQ=255)
                                            ccGAAGTTGGTATCGGTGACCTATGGCGCGAACGCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTg                                             <  1:3407110/93‑1 (MQ=255)
                                             cGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTgg                                            >  1:6043635/1‑93 (MQ=255)
                                              gAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                           <  2:4091822/93‑1 (MQ=255)
                                               aaGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGag                                          >  1:2964596/1‑92 (MQ=255)
                                               aaGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGaa                                          >  2:6009230/1‑93 (MQ=255)
                                               aaGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGaa                                          <  1:2733789/93‑1 (MQ=255)
                                               aaGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGaa                                          >  2:3030318/1‑93 (MQ=255)
                                                aGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAg                                         <  2:382656/93‑1 (MQ=255)
                                                aGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAg                                         >  1:6618096/1‑93 (MQ=255)
                                                 gTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGc                                        <  2:3329279/93‑1 (MQ=255)
                                                          ggTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCa                               <  2:791880/93‑1 (MQ=255)
                                                          ggTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCa                               <  1:2609086/93‑1 (MQ=255)
                                                          ggTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCa                               >  1:7009735/1‑93 (MQ=255)
                                                           gTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCAt                              >  1:5745583/1‑93 (MQ=255)
                                                           gTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCAt                              >  1:7849613/1‑93 (MQ=255)
                                                             gACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCATCGCATCt                            >  1:2502599/1‑93 (MQ=255)
                                                                cTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTAc                         >  1:4065598/1‑93 (MQ=255)
                                                                 tATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCTCTGGTCTGGAAGCGGCACCGCATCTTACt                        <  2:5930736/93‑1 (MQ=255)
                                                                  aTGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATTGCACTGGTCTGGTAGCGGCACCGCATCTTACtt                       >  2:4739472/1‑93 (MQ=255)
                                                                  aTGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACtt                       <  1:8067062/93‑1 (MQ=255)
                                                                    ggCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGc                     >  2:2823754/1‑93 (MQ=255)
                                                                     gcgcGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCa                    <  2:2451167/93‑1 (MQ=255)
                                                                     gcgcGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCa                    <  1:1383083/93‑1 (MQ=255)
                                                                       gcgAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCAtt                  <  1:1421512/93‑1 (MQ=255)
                                                                        cgAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTg                 <  1:1603843/93‑1 (MQ=255)
                                                                         gAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGa                <  2:2239446/93‑1 (MQ=255)
                                                                          aaCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt               >  1:3729452/1‑93 (MQ=255)
                                                                          aaCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt               >  1:6475031/1‑93 (MQ=255)
                                                                          aaCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt               <  2:4829769/93‑1 (MQ=255)
                                                                           aCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATg              >  1:6069343/1‑93 (MQ=255)
                                                                             tCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCg            <  2:7220137/93‑1 (MQ=255)
                                                                               cGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGAc          <  2:4501740/93‑1 (MQ=255)
                                                                               cGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGAc          >  2:2695865/1‑93 (MQ=255)
                                                                                ggCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACg         >  1:2020612/1‑93 (MQ=255)
                                                                                 gCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAATCGGCACCGCATCTTACTTGCATTGATGCGACGc        <  1:2124612/93‑1 (MQ=255)
                                                                                  cGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGcc       <  1:1221732/93‑1 (MQ=255)
                                                                                   gAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGccc      <  1:3281177/93‑1 (MQ=255)
                                                                                   gAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGccc      >  1:7569146/1‑93 (MQ=255)
                                                                                    aGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcg     >  1:1827144/1‑93 (MQ=255)
                                                                                    aGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcg     >  2:3613940/1‑93 (MQ=255)
                                                                                     gcgcGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcgc    >  2:2375075/1‑92 (MQ=255)
                                                                                          |                                                                                         
AGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACG  >  NC_000913/4130308‑4130487

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