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.4% 527.0 / ‑5.1 147D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (76/70);  total (77/70)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 147 total aligned reads displayed.

CGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGC  >  NC_000913/3007980‑3008172
                                                                                               |                                                                                                 
cGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACg                                                                                               <  2:2651633/100‑1 (MQ=255)
 ggTGCCGAGTGCATCGTCACGGATATGAACTATGACGCTACCGTTCGCCTGAACATGCAACACGAGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACgg                                                                                              >  2:1085362/1‑100 (MQ=255)
  gTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGc                                                                                             <  1:3754881/100‑1 (MQ=255)
    gCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCAt                                                                                           >  2:101809/1‑100 (MQ=255)
     ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                          <  2:1771030/100‑1 (MQ=255)
     ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                          <  1:3585730/100‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                         <  2:3284002/100‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                         <  1:2684588/100‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                         >  2:1157660/1‑100 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                         <  1:158448/100‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                         >  2:3190936/1‑100 (MQ=255)
         gTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGaa                                                                                      <  2:2371781/100‑1 (MQ=255)
          agCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGTCCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAg                                                                                     <  1:2664478/99‑1 (MQ=255)
           gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGCCCATGCAACACGCGCAGCAGCACGGCTGGGACGTGGTACAGGGCACGCCATGGGAAg                                                                                     >  2:1643090/1‑99 (MQ=255)
           gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                    >  1:164737/1‑100 (MQ=255)
            cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                   >  2:3048846/1‑100 (MQ=255)
             aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                                  >  1:3310043/1‑100 (MQ=255)
             aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                                  >  1:2056022/1‑100 (MQ=255)
               cGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                                >  1:2629768/1‑100 (MQ=255)
                gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                               >  2:2458145/1‑100 (MQ=255)
                 tCACGGATATGAACTATGACGATACCGGTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                              <  1:3582214/100‑1 (MQ=255)
                  cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                             >  2:2949314/1‑100 (MQ=255)
                  cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                             >  2:1445925/1‑100 (MQ=255)
                   aCGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCa                                                                            >  2:1412979/1‑100 (MQ=255)
                      gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCacaa                                                                         >  1:2612245/1‑100 (MQ=255)
                        tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                       >  2:1684514/1‑100 (MQ=255)
                         atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                      >  2:724666/1‑100 (MQ=255)
                          tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                     >  1:275015/1‑100 (MQ=255)
                          tGAACTATGACAATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTAGAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                     >  1:1996525/1‑100 (MQ=255)
                           gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCa                                                                    >  1:3452736/1‑100 (MQ=255)
                            aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                   <  1:2574705/100‑1 (MQ=255)
                            aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                   <  1:2574636/100‑1 (MQ=255)
                                aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                               >  1:336069/1‑100 (MQ=255)
                                  gACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGa                                                             >  2:3360211/1‑100 (MQ=255)
                                   aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                            <  2:2072618/100‑1 (MQ=255)
                                   aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                            <  1:938425/100‑1 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                           <  1:3247707/100‑1 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                           <  1:2996379/100‑1 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                           <  1:1501580/100‑1 (MQ=255)
                                      aTACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                           <  1:633133/98‑1 (MQ=255)
                                        aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                       <  1:569711/100‑1 (MQ=255)
                                        aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                       <  1:724672/100‑1 (MQ=255)
                                        aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                       <  1:182092/100‑1 (MQ=255)
                                         ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTCCAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                       >  2:3437486/1‑99 (MQ=255)
                                          cGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCaa                                                     >  1:2379946/1‑100 (MQ=255)
                                           gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                    >  2:439358/1‑100 (MQ=255)
                                             tCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                                  >  1:2912709/1‑100 (MQ=255)
                                               gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                                <  1:1149073/100‑1 (MQ=255)
                                               gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                                <  2:2695082/100‑1 (MQ=255)
                                                ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                               >  2:426146/1‑100 (MQ=255)
                                                ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                               >  2:3569777/1‑100 (MQ=255)
                                                ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                               >  1:1014314/1‑100 (MQ=255)
                                                 cTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACg                                              <  1:2875831/100‑1 (MQ=255)
                                                 cTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACg                                              <  2:1894406/100‑1 (MQ=255)
                                                  tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                             >  1:403502/1‑100 (MQ=255)
                                                  tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                             >  1:2172501/1‑100 (MQ=255)
                                                          cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAACATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                     >  1:414071/1‑100 (MQ=255)
                                                          cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                     <  1:2138550/100‑1 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                      <  2:590227/98‑1 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                    <  1:1821586/100‑1 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                    >  2:3046224/1‑100 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                    >  2:127736/1‑100 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                    <  2:666586/100‑1 (MQ=255)
                                                            acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                   <  2:233343/100‑1 (MQ=255)
                                                              acGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGa                                 >  1:836069/1‑100 (MQ=255)
                                                              aagcgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGa                                 >  2:2573168/3‑100 (MQ=255)
                                                               cgcgCAGCAGCACGGCTGGGAAGTGGTGCAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                                >  2:3498169/1‑100 (MQ=255)
                                                               cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                                <  2:340745/100‑1 (MQ=255)
                                                                gcgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATg                               <  2:1378793/100‑1 (MQ=255)
                                                                 cgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGa                              <  2:637853/100‑1 (MQ=255)
                                                                    agcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                           <  2:209519/100‑1 (MQ=255)
                                                                    agcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                           <  2:2947244/100‑1 (MQ=255)
                                                                     gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                          >  2:1751054/1‑100 (MQ=255)
                                                                      cagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACCCCAAAATCCCAACCTGGATCCTGCAAGGCTTCGCAAACCTGGCAGATGCAGCCg                         <  1:1503956/100‑1 (MQ=255)
                                                                      cagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCg                         >  1:2808734/1‑100 (MQ=255)
                                                                       agcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGt                        <  2:1781769/100‑1 (MQ=255)
                                                                         caCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                      >  2:854696/1‑100 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGa                     >  1:1344478/1‑99 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                    <  1:1412218/100‑1 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                    <  1:1085357/100‑1 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                    <  2:2630712/100‑1 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                    <  2:3442939/100‑1 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                    <  1:565475/100‑1 (MQ=255)
                                                                            ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc                   >  2:949430/1‑100 (MQ=255)
                                                                            ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc                   >  1:955711/1‑100 (MQ=255)
                                                                             gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                  <  1:999675/100‑1 (MQ=255)
                                                                             gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                  <  2:778547/100‑1 (MQ=255)
                                                                                gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACTTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGTCGTCGa                     >  1:831099/1‑94 (MQ=255)
                                                                                gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt               >  2:1283640/1‑100 (MQ=255)
                                                                                 ggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATg              >  2:2230173/1‑100 (MQ=255)
                                                                                  gAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGc             >  1:3476517/1‑100 (MQ=255)
                                                                                   aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg            <  2:268289/100‑1 (MQ=255)
                                                                                   aaGTCGTACAGCGCACGACATTCCAAGTCTACACCAACAGCCCCACCTGGATCATCCAAGGCTACGCAACCCTCTCCCATGAATCCGTCTAGCAAATGCg            <  1:2382813/100‑1 (MQ=255)
                                                                                     gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg          <  1:2430279/100‑1 (MQ=255)
                                                                                      tGGTACATGGCACCGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGa         <  1:3322476/100‑1 (MQ=255)
                                                                                       ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa        >  2:3104287/1‑100 (MQ=255)
                                                                                       ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa        >  2:1897994/1‑100 (MQ=255)
                                                                                        gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa       >  2:561720/1‑100 (MQ=255)
                                                                                          aCAGGGCACGGCATGGGAAGGTTACACCATAAGCCCAACCTGGATCAGGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGa         <  1:3676461/96‑1 (MQ=255)
                                                                                             ggACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGc  >  1:996574/1‑100 (MQ=255)
                                                                                               |                                                                                                 
CGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGGC  >  NC_000913/3007980‑3008172

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