Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,008,075 A→G 100% D189G (GAC→GGC)  ygeX → 2,3‑diaminopropionate ammonia lyase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,008,0750AG99.3% 461.6 / ‑5.0 137D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (71/64);  total (72/65)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

Only 100 of 139 total aligned reads displayed.

CTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007978‑3008170
                                                                                                 |                                                                                               
cTCGGTGCCGAGTGCATCGTCACGGATATAAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGgca                                                                                               >  2:683747/1‑100 (MQ=255)
 tCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGgcac                                                                                              >  1:3106882/1‑100 (MQ=255)
    gTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGc                                                                                           <  1:3318218/100‑1 (MQ=255)
     tGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACCCGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCa                                                                                          >  1:3891690/1‑100 (MQ=255)
       ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                        >  1:1106535/1‑100 (MQ=255)
        cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       <  1:73577/100‑1 (MQ=255)
        cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       >  1:3360758/1‑100 (MQ=255)
        cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       >  1:3551112/1‑100 (MQ=255)
          aGTGCATCGTCACGGATATGAATTATGACGATATCGTTCGCCTGACTATGGAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGa                                                                                     <  1:596008/100‑1 (MQ=255)
          aGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGa                                                                                     <  2:2690064/100‑1 (MQ=255)
           gTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGaa                                                                                    >  1:2392909/1‑100 (MQ=255)
            tGCATCGTCACGGATATGAACTATGACGATAACGTTCGCCTGACCATGCAACACGCGCATCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGaa                                                                                    <  1:1379271/99‑1 (MQ=255)
             gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  <  1:2722925/100‑1 (MQ=255)
             gCATCGTCACGGATATGAACTATGACGATACCGGTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGCCATGGGAAg                                                                                   >  2:2575702/1‑99 (MQ=255)
                tCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTa                                                                               <  2:408679/100‑1 (MQ=255)
                tCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTa                                                                               <  2:3263429/100‑1 (MQ=255)
                  gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             <  2:1386865/100‑1 (MQ=255)
                  gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  1:2879972/1‑100 (MQ=255)
                    cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                            >  1:2803481/1‑99 (MQ=255)
                    cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                           <  2:2858878/100‑1 (MQ=255)
                          tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                     <  2:3828834/100‑1 (MQ=255)
                            tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  2:3800276/1‑100 (MQ=255)
                            tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  2:2315822/1‑100 (MQ=255)
                            tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                    <  2:2159926/99‑1 (MQ=255)
                             gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCa                                                                  >  2:1002237/1‑100 (MQ=255)
                              aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                 <  2:3564033/100‑1 (MQ=255)
                              aaCTATGACGATACCGATCGCCTGACCATGCAACACGCGCAGCAGCACGGATGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                 >  2:739076/1‑100 (MQ=255)
                                cTATGACGATACCGTTCGCCTTACCATGCAATACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGTCACGGCATGGGAAGGTTACACCAAAATCCCAAcc                                                               <  1:1670737/100‑1 (MQ=255)
                                 tATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                <  1:1002248/98‑1 (MQ=255)
                                 tATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCt                                                              >  1:3328262/1‑100 (MQ=255)
                                  aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                             <  1:2767735/100‑1 (MQ=255)
                                  aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                             <  1:3828219/100‑1 (MQ=255)
                                  aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                             <  1:3800315/100‑1 (MQ=255)
                                   tGACGATACCTTTCGCCTGACCATGCAACGCGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAcn                                                               <  1:1796396/97‑2 (MQ=255)
                                   tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                            <  2:872095/100‑1 (MQ=255)
                                   tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                            >  2:3604650/1‑100 (MQ=255)
                                    gACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGa                                                           >  1:1029761/1‑100 (MQ=255)
                                     aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGa                                                           >  1:3073281/1‑99 (MQ=255)
                                         tACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATg                                                      <  1:3685753/100‑1 (MQ=255)
                                           ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                    >  2:1189452/1‑100 (MQ=255)
                                           ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                    >  1:3919041/1‑100 (MQ=255)
                                           ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                    >  2:1820529/1‑100 (MQ=255)
                                           ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                    <  1:2526304/100‑1 (MQ=255)
                                             cttccgatctccATGCAACCCGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  <  2:1696501/90‑1 (MQ=255)
                                              ttccgatctccATGCAATACGCGCAGCAGCACGGCTGTGAAGTTGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTAGATCATGCAAg                                                  <  2:728261/90‑1 (MQ=255)
                                               tCGCCTGACCATGCAACACGCGCAGCAGCACTGCTGGGAAGTACTACAGTGCACGGCATGTGAAGTTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  <  2:1016764/98‑1 (MQ=255)
                                               tCGCCTGACCATGCAACACGCGCAGCAGCACGGTTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                                <  1:597438/100‑1 (MQ=255)
                                                cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                               >  2:2677543/1‑100 (MQ=255)
                                                  ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                             >  1:790824/1‑100 (MQ=255)
                                                  ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCAGGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                             >  2:393058/1‑100 (MQ=255)
                                                    tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                           <  1:1043502/100‑1 (MQ=255)
                                                    tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGCCTACGc                                           >  2:3043524/1‑100 (MQ=255)
                                                    tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATACAAGGCTAcac                                           >  1:1811213/1‑98 (MQ=255)
                                                     gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          <  2:1972270/100‑1 (MQ=255)
                                                     gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          >  1:3806393/1‑100 (MQ=255)
                                                     gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          >  2:865639/1‑100 (MQ=255)
                                                     gACCATGCAACACGCGCAGCAGCACGGCTGGGAACTGGTACAGCGCACGGCATGGGAAGGTTACTCCAAACTCCCAACCTGGATCATGCAAGGCTACGc                                           <  2:3601916/99‑1 (MQ=255)
                                                     gACCATGCAACACGCGCAGCAGCACGACTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          >  2:2457558/1‑100 (MQ=255)
                                                      aCCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCaa                                         <  1:836359/100‑1 (MQ=255)
                                                       ccATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAgagatcggaag                                        >  1:1696516/1‑90 (MQ=255)
                                                       ccATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAc                                        >  1:1606583/1‑100 (MQ=255)
                                                       ccATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAc                                        >  2:588691/1‑100 (MQ=255)
                                                       ccATGAAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCGGGATTATGCAAgagatcggaag                                        >  1:728270/1‑90 (MQ=255)
                                                        cATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAcc                                       >  1:839903/1‑100 (MQ=255)
                                                           gCAACACGCGCAGCAGCACGGGTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                    <  1:624499/100‑1 (MQ=255)
                                                           gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                    <  1:3870910/100‑1 (MQ=255)
                                                           gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                    <  2:1300279/100‑1 (MQ=255)
                                                           gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCAAACCTGGATCATGCAAGGCTACGCAACCCTg                                    >  2:3187951/1‑100 (MQ=255)
                                                            cAACACGCGCAGCAGTACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                   >  2:3466388/1‑100 (MQ=255)
                                                            cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                   >  2:123499/1‑100 (MQ=255)
                                                            cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTga                                   >  1:2554067/1‑99 (MQ=255)
                                                             aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                  >  2:1746842/1‑100 (MQ=255)
                                                             aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                  >  1:350421/1‑100 (MQ=255)
                                                             aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                  >  2:2995475/1‑100 (MQ=255)
                                                              acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                 <  1:1957174/100‑1 (MQ=255)
                                                              acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                 <  1:1269670/100‑1 (MQ=255)
                                                                 cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                              >  2:2798406/1‑100 (MQ=255)
                                                                 cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                              >  1:3474708/1‑100 (MQ=255)
                                                                 cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATCGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                              >  2:2137796/1‑100 (MQ=255)
                                                                     cagcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAg                          >  1:2430010/1‑100 (MQ=255)
                                                                      agcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAAACCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                         <  1:3905299/100‑1 (MQ=255)
                                                                       gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                        >  1:1012459/1‑100 (MQ=255)
                                                                          gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCTTGGCAGATGAAGCCGTc                     <  2:609480/100‑1 (MQ=255)
                                                                           caCGGCTGGGAAGTGGTACAGGGAACGGCATGGGAAGGTTACACCAAAATCCAAACCTGGATAATGCAAGGCTACGCAACCCTGGCAGATGAAGCAGtag                    >  2:510801/1‑98 (MQ=255)
                                                                             cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  <  1:1976009/100‑1 (MQ=255)
                                                                               gatGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                <  1:2808837/98‑1 (MQ=255)
                                                                               gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                <  1:1323249/100‑1 (MQ=255)
                                                                               gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                <  2:2437954/100‑1 (MQ=255)
                                                                               gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                <  2:2021318/100‑1 (MQ=255)
                                                                                  gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt             >  2:3265858/1‑100 (MQ=255)
                                                                                  gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt             <  1:2469478/100‑1 (MQ=255)
                                                                                    gAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGc           <  1:2677567/100‑1 (MQ=255)
                                                                                      aGTGGTACAGGGCACGGCATGGGAAGGTTAAACAAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGGGCAAATGCg          >  2:2572293/1‑99 (MQ=255)
                                                                                       gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg        <  2:2992925/100‑1 (MQ=255)
                                                                                          gTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  2:1372005/1‑100 (MQ=255)
                                                                                           tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     <  2:2732789/99‑1 (MQ=255)
                                                                                           tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt    >  2:3627302/1‑100 (MQ=255)
                                                                                           tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt    >  1:3434766/1‑100 (MQ=255)
                                                                                           tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt    >  2:1677094/1‑100 (MQ=255)
                                                                                             cAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATgg  >  1:2033982/1‑100 (MQ=255)
                                                                                                 |                                                                                               
CTCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007978‑3008170

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