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,3980AG100.0% 836.8 / NA 240D68G (GAC→GGC) metF5,10‑methylenetetrahydrofolate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (126/114);  total (126/114)

Only 100 of 246 total aligned reads displayed.

AGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACG  >  NC_000913/4130308‑4130487
                                                                                          |                                                                                         
 gCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCg                                                                                        >  2:1533908/1‑93 (MQ=255)
   aGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTa                                                                                      <  2:6159472/93‑1 (MQ=255)
    gACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTAc                                                                                     >  2:4026131/1‑93 (MQ=255)
    gACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTAc                                                                                     >  1:6060053/1‑93 (MQ=255)
      cccTGTGGAACTCCATCGATCGCCTTAGCAGCCTGACACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGc                                                                                   <  2:2879686/93‑1 (MQ=255)
       ccTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGca                                                                                  <  1:323936/93‑1 (MQ=255)
        cTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGcac                                                                                 <  2:887680/93‑1 (MQ=255)
        cTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGcac                                                                                 <  1:5054785/93‑1 (MQ=255)
           tgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGc                                                                              <  1:5679014/93‑1 (MQ=255)
           tgGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGc                                                                              <  2:377445/93‑1 (MQ=255)
            ggAACTCCATCGATCTCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCa                                                                             <  1:4729761/93‑1 (MQ=255)
            ggAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCa                                                                             >  1:5960088/1‑93 (MQ=255)
            ggAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCa                                                                             <  2:3131992/93‑1 (MQ=255)
            ggAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCa                                                                             >  1:4109180/1‑93 (MQ=255)
             gaccTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTTACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCat                                                                            <  1:1810864/90‑1 (MQ=255)
             gAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCat                                                                            <  1:3877133/93‑1 (MQ=255)
              aaCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCatt                                                                           <  1:4893041/93‑1 (MQ=255)
               aCTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCatta                                                                          >  1:2430813/1‑93 (MQ=255)
                cTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCattat                                                                         <  2:2130513/93‑1 (MQ=255)
                 tCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCattatt                                                                        >  2:4617930/1‑93 (MQ=255)
                  ccATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCAttatta                                                                       >  2:1850328/1‑93 (MQ=255)
                   cATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTaa                                                                      <  1:3391084/93‑1 (MQ=255)
                   cATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTaa                                                                      >  2:4559089/1‑93 (MQ=255)
                    atcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTaaa                                                                     <  2:4702368/93‑1 (MQ=255)
                     tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                    >  1:378245/1‑93 (MQ=255)
                     tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                    <  2:2180009/93‑1 (MQ=255)
                     tcgatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAg                                                                    >  1:2225012/1‑93 (MQ=255)
                       gatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGc                                                                  >  1:1783514/1‑93 (MQ=255)
                       gatcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGc                                                                  <  2:4444040/93‑1 (MQ=255)
                        atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                 >  2:2546841/1‑93 (MQ=255)
                        atcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCa                                                                 >  1:4646410/1‑93 (MQ=255)
                         tcgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAt                                                                <  1:5199929/93‑1 (MQ=255)
                          cgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAtt                                                               <  1:5431535/93‑1 (MQ=255)
                          cgCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCAtt                                                               >  1:2571834/1‑93 (MQ=255)
                             cTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTaaa                                                            >  2:1976694/1‑93 (MQ=255)
                              ttAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                           >  2:1928827/1‑93 (MQ=255)
                              ttAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                           >  1:582567/1‑93 (MQ=255)
                              ttAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAg                                                           >  2:1330990/1‑93 (MQ=255)
                               tAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGa                                                          <  1:4397062/93‑1 (MQ=255)
                                 gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                        <  2:3997631/93‑1 (MQ=255)
                                 gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                        <  2:2450505/93‑1 (MQ=255)
                                 gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                        <  1:3091395/93‑1 (MQ=255)
                                 gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                        >  1:1171949/1‑93 (MQ=255)
                                 gcagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATc                                                        >  1:6192474/1‑93 (MQ=255)
                                  cagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCg                                                       <  1:2519288/93‑1 (MQ=255)
                                  cagcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCg                                                       >  2:4087600/1‑93 (MQ=255)
                                   agcCTGAAGCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTCCGCACAGCATTATTAAAGGCATTAAAGATCGc                                                      <  2:5207827/93‑1 (MQ=255)
                                   agcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGc                                                      >  1:975710/1‑93 (MQ=255)
                                   agcCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGc                                                      >  2:4804380/1‑93 (MQ=255)
                                      cTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACt                                                   >  1:6077137/1‑93 (MQ=255)
                                         aaaCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGt                                                <  1:4226360/93‑1 (MQ=255)
                                         aaaCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGt                                                >  2:6362375/1‑93 (MQ=255)
                                           aCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCt                                              >  2:3717426/1‑93 (MQ=255)
                                           aCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCt                                              >  2:49084/1‑93 (MQ=255)
                                           aCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCt                                              <  1:6300699/93‑1 (MQ=255)
                                           aCCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCt                                              <  1:292439/93‑1 (MQ=255)
                                            ccGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGAGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTg                                             >  1:1066316/1‑93 (MQ=255)
                                              gAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                           <  2:3594961/93‑1 (MQ=255)
                                              gAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGa                                           >  2:5319460/1‑93 (MQ=255)
                                               aaGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGaa                                          >  2:5952593/1‑93 (MQ=255)
                                               aaGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGaa                                          >  2:3415601/1‑93 (MQ=255)
                                                aGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAg                                         >  2:5912165/1‑93 (MQ=255)
                                                aGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAg                                         >  2:4254223/1‑93 (MQ=255)
                                                 gTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGc                                        >  2:5467504/1‑93 (MQ=255)
                                                  tttGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCg                                       <  2:1896952/93‑1 (MQ=255)
                                                  tttGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCg                                       <  1:2832768/93‑1 (MQ=255)
                                                    tGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGc                                     <  1:4085028/93‑1 (MQ=255)
                                                    tGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGc                                     >  1:1364978/1‑93 (MQ=255)
                                                      tATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAc                                   <  1:4576571/93‑1 (MQ=255)
                                                      tATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAc                                   <  2:2744164/93‑1 (MQ=255)
                                                       aTCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAcc                                  >  2:5136066/1‑93 (MQ=255)
                                                       aTCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCAcc                                  <  2:363999/93‑1 (MQ=255)
                                                        tCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCg                                 >  2:3525517/1‑93 (MQ=255)
                                                        tCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCg                                 <  1:553271/93‑1 (MQ=255)
                                                        tCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCg                                 <  2:5273534/93‑1 (MQ=255)
                                                           gTGACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCAt                              >  2:5986714/1‑93 (MQ=255)
                                                             gACCTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCt                            >  2:1702115/1‑93 (MQ=255)
                                                               ccTATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTa                          >  1:917540/1‑93 (MQ=255)
                                                                 tATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACt                        <  1:4347994/93‑1 (MQ=255)
                                                                 tATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACt                        >  2:4334436/1‑93 (MQ=255)
                                                                 tATGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACt                        <  1:170938/93‑1 (MQ=255)
                                                                  aTGGCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACtt                       >  2:5771682/1‑93 (MQ=255)
                                                                    ggCGCGAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGc                     >  2:3291492/1‑93 (MQ=255)
                                                                       gcgAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCAtt                  >  2:6134169/1‑93 (MQ=255)
                                                                         gAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGa                <  2:1312206/93‑1 (MQ=255)
                                                                         gAACTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGa                <  1:3647962/93‑1 (MQ=255)
                                                                          aaCTCCGGCGAGCTCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACGGCATCTTACTTGCATTGAt               >  2:929779/1‑93 (MQ=255)
                                                                          aaCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTACAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt               <  1:2182508/93‑1 (MQ=255)
                                                                          aaCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt               <  1:3579006/93‑1 (MQ=255)
                                                                          aaCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt               >  1:3431872/1‑93 (MQ=255)
                                                                          aaCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGAt               <  1:1957596/93‑1 (MQ=255)
                                                                          aaCTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGCCTGGAAGCGGCACCGCATCTTACTTGCATTGAt               <  2:2268398/93‑1 (MQ=255)
                                                                            cTCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGc             <  1:5859723/93‑1 (MQ=255)
                                                                             tCCGGCGAGCGCGGCCGTACGCATAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCg            <  2:753152/93‑1 (MQ=255)
                                                                             tCCGGCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCg            >  2:2967626/1‑93 (MQ=255)
                                                                                ggCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGCTCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACg         <  1:6377624/93‑1 (MQ=255)
                                                                                 gCGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGc        <  1:691770/93‑1 (MQ=255)
                                                                                  cGAGCGCGGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGcc       <  1:5898586/93‑1 (MQ=255)
                                                                                      cgcgGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcgac   >  1:1998537/1‑93 (MQ=255)
                                                                                       gcgGCCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCcgacg  >  1:101890/1‑93 (MQ=255)
                                                                                          |                                                                                         
AGCAGACCCTGTGGAACTCCATCGATCGCCTTAGCAGCCTGAAACCGAAGTTTGTATCGGTGACCTATGGCGCGAACTCCGGCGAGCGCGACCGTACGCACAGCATTATTAAAGGCATTAAAGATCGCACTGGTCTGGAAGCGGCACCGCATCTTACTTGCATTGATGCGACGCCCGACG  >  NC_000913/4130308‑4130487

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