Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,930,097 (A)7→6 intergenic (+14/+42) purT → / ← eda phosphoribosylglycinamide formyltransferase 2/2‑keto‑3‑deoxygluconate 6‑phosphate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,930,0910A.99.5% 829.4 / ‑3.2 181intergenic (+8/+48)purT/edaphosphoribosylglycinamide formyltransferase 2/2‑keto‑3‑deoxygluconate 6‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base . (80/100);  total (80/101)
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 196 total aligned reads displayed.

GGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAA  >  NC_000913/1929999‑1930182
                                                                                            |                                                                                           
     cACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaa                                                                                        >  2:5606364/1‑93 (MQ=255)
      aCTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaac                                                                                       <  1:5146214/93‑2 (MQ=255)
        tGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCaaaaaacgc                                                                                     >  1:5212558/1‑90 (MQ=255)
          gCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTa                                                                                  <  1:1065963/93‑1 (MQ=255)
           cTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTAc                                                                                 <  1:5673345/93‑1 (MQ=255)
           cTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTAc                                                                                 <  1:2639957/93‑1 (MQ=255)
           cTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTAc                                                                                 <  1:1870159/93‑1 (MQ=255)
              cTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaa                                                                              >  1:4050792/1‑93 (MQ=255)
               tGCAGAGAGTGTTGTTGACGCCATGGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaaa                                                                             <  2:3554821/93‑1 (MQ=255)
                gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaaaa                                                                            <  1:5592191/93‑1 (MQ=255)
                gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaaaa                                                                            <  1:3070607/93‑1 (MQ=255)
                gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaaaa                                                                            <  1:2279548/93‑1 (MQ=255)
                gCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACaaaaa                                                                            <  2:3587389/93‑1 (MQ=255)
                  agagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATg                                                                          >  2:1352238/1‑93 (MQ=255)
                  agagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATg                                                                          <  1:4612252/93‑1 (MQ=255)
                   gagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAAttc                                                                         >  2:271865/1‑91 (MQ=255)
                   gagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGc                                                                         >  1:3791042/1‑93 (MQ=255)
                   gagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGc                                                                         >  1:2720668/1‑93 (MQ=255)
                   gagagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGc                                                                         <  1:1078581/93‑1 (MQ=255)
                    agcgtggtgTTGACGCCATTGAACGCGCGAAGCACGCCGCCGTACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGcc                                                                        <  2:219676/90‑1 (MQ=255)
                    agagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGcc                                                                        <  2:5717495/93‑1 (MQ=255)
                    agagTGTTGTCGACGCCATTGAACGCGCGAAGCACGCCGCCGGAGAGATAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGcc                                                                        <  1:2189389/93‑1 (MQ=255)
                     gagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGccc                                                                       <  2:5505314/93‑1 (MQ=255)
                     gagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGccc                                                                       <  2:5131365/93‑1 (MQ=255)
                     gagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGccc                                                                       <  1:933959/93‑1 (MQ=255)
                     gagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGccc                                                                       <  2:1932247/93‑1 (MQ=255)
                     gagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGccc                                                                       >  2:439051/1‑93 (MQ=255)
                     gagTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGccc                                                                       <  2:402268/93‑1 (MQ=255)
                      agTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCg                                                                      <  1:1617946/93‑1 (MQ=255)
                        tgttgtTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGAt                                                                    <  1:6168062/93‑1 (MQ=255)
                         gttgttGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATc                                                                   >  1:639569/1‑93 (MQ=255)
                          ttgttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATcc                                                                  <  2:2922997/93‑1 (MQ=255)
                            gttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTc                                                                <  2:520789/93‑1 (MQ=255)
                             ttgACGCCATTGAACGCGGGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCg                                                               <  1:2329456/93‑1 (MQ=255)
                             ttgACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCg                                                               >  2:4987550/1‑93 (MQ=255)
                               gACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGAt                                                             <  1:1806448/93‑1 (MQ=255)
                                aCGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATc                                                            <  2:5490589/93‑1 (MQ=255)
                                  gCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATACTCCATCgg                                                          >  2:3435060/1‑93 (MQ=255)
                                   ccATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATACCCGATCCTCGATCggg                                                         >  2:302744/1‑93 (MQ=255)
                                   ccATTGAAAGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCggg                                                         >  2:6119175/1‑93 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAATCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                       <  1:3219229/93‑1 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAACCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                       <  1:3698483/93‑1 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                       >  2:1877837/1‑93 (MQ=255)
                                     aTTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCa                                                       >  1:820908/1‑93 (MQ=255)
                                      ttGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAt                                                      <  2:2483288/93‑1 (MQ=255)
                                        gAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttt                                                    <  2:5686781/93‑1 (MQ=255)
                                        gAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAttt                                                    <  1:1754033/93‑1 (MQ=255)
                                         aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAACCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                   <  2:2762313/93‑1 (MQ=255)
                                         aaCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCAtttt                                                   <  1:3080269/93‑1 (MQ=255)
                                          aCGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTg                                                  >  1:4878514/1‑93 (MQ=255)
                                            gcgcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGAc                                                <  1:788173/93‑1 (MQ=255)
                                              gcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtt                                              >  1:5828456/1‑93 (MQ=255)
                                              gcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtt                                              <  1:2520726/93‑1 (MQ=255)
                                              gcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtt                                              <  2:3208777/93‑1 (MQ=255)
                                              gcgAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACtt                                              <  1:1573837/93‑1 (MQ=255)
                                                 aaGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                           <  1:1007272/93‑1 (MQ=255)
                                                 aaGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTa                                           <  1:6293889/93‑1 (MQ=255)
                                                  aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                          >  1:1368733/1‑93 (MQ=255)
                                                  aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTAc                                          <  2:5034182/93‑1 (MQ=255)
                                                  aGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATACTCGATCGGGCATTTTGACTTTTAc                                          >  1:4258976/1‑93 (MQ=255)
                                                   gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACa                                         >  2:450312/1‑93 (MQ=255)
                                                   gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACa                                         <  1:1447852/93‑1 (MQ=255)
                                                   gCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACa                                         <  1:4745681/93‑1 (MQ=255)
                                                       gccgccGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCtt                                     <  1:1848827/93‑1 (MQ=255)
                                                        ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTa                                    >  1:1337468/1‑93 (MQ=255)
                                                        ccgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTa                                    <  2:1760752/93‑1 (MQ=255)
                                                        ccgccgAACAGGTAAAAGTACAGGGTTAACCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTa                                    <  2:2343826/93‑1 (MQ=255)
                                                         cgccgGACAGGTAAAATTACAGGGTTAAACCAGGCAAAAAACGCTACAAAGATGCCGGATCCTCGATCGGGCGTTTTGACTTTTACAGCTTAg                                   <  1:3388078/93‑1 (MQ=255)
                                                         cgccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAg                                   >  2:2101923/1‑93 (MQ=255)
                                                          gccgGACAGGTAAAAGTACAGGTTTAACCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgc                                  <  1:5712025/93‑1 (MQ=255)
                                                          gccgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTTTCTTTTTTATCTTAgc                                  >  1:4381397/1‑93 (MQ=255)
                                                            cgGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAgcgc                                >  2:4089241/1‑93 (MQ=255)
                                                             ggACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGcc                               >  2:651519/1‑93 (MQ=255)
                                                              gACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCt                              <  2:1601320/93‑1 (MQ=255)
                                                               aCAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCtt                             >  1:5031942/1‑93 (MQ=255)
                                                                 aGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCctct                           >  2:4630508/1‑93 (MQ=255)
                                                                 aGGTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCt                           <  2:838865/93‑1 (MQ=255)
                                                                  ggTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTa                          >  2:3789525/1‑93 (MQ=255)
                                                                  ggTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTa                          <  1:1028766/93‑1 (MQ=255)
                                                                   gTAACAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTTCAGCTTAGAGCCTTCtcc                         >  2:5985743/1‑91 (MQ=255)
                                                                   gTAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTAc                         >  2:687169/1‑93 (MQ=255)
                                                                    tAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                        <  1:5411155/93‑1 (MQ=255)
                                                                    tAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACa                        <  2:4894150/93‑1 (MQ=255)
                                                                    tAAAAGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCCTCTACa                        >  2:3599874/1‑93 (MQ=255)
                                                                     aaaaGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAg                       <  2:1638530/93‑1 (MQ=255)
                                                                      aaaGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGc                      <  1:3628551/93‑1 (MQ=255)
                                                                      aaaGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGc                      <  1:2861623/93‑1 (MQ=255)
                                                                       aaGTACAGGGTTAAACCCGGCAAAAACCGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCt                     <  1:5329970/93‑1 (MQ=255)
                                                                        aGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCtt                    >  1:5871710/1‑93 (MQ=255)
                                                                        aGTACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCCGGCACTTTTACTTTTTCAGCTTAGCGCCTTCTACAGCtt                    >  1:2460811/1‑93 (MQ=255)
                                                                          tACAGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCa                  <  1:4201398/93‑1 (MQ=255)
                                                                             aGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgc               >  1:4939076/1‑93 (MQ=255)
                                                                             aGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgc               <  2:5657709/93‑1 (MQ=255)
                                                                             aGGGTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCAcgc               <  2:6206505/93‑1 (MQ=255)
                                                                              gggTTAAACCCGGCAAAAAACGCTACAAAAATGCCCCATCCCCCATCGGGGCTTTTTACTTTTTTAACTTAGCGCCTTCTTCAACTTCAcccg              >  2:2282436/1‑93 (MQ=255)
                                                                               ggTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGCGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGATTCACgcgc             >  1:2983062/1‑93 (MQ=255)
                                                                                gTTAAACCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGcc            >  1:4764327/1‑93 (MQ=255)
                                                                                    aaCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCt        >  2:2051061/1‑93 (MQ=255)
                                                                                     aCCCGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCtt       >  1:2943082/1‑93 (MQ=255)
                                                                                      cccGGCAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTa      >  2:3095040/1‑93 (MQ=255)
                                                                                            |                                                                                           
GGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCCCGATCCTCGATCGGGCATTTTGACTTTTACAGCTTAGCGCCTTCTACAGCTTCACGCGCCAGCTTAGTAA  >  NC_000913/1929999‑1930182

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