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,0750AG97.8% 518.6 / ‑0.8 155D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (73/78);  total (76/79)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 155 total aligned reads displayed.

CGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007980‑3008170
                                                                                               |                                                                                               
cGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACg                                                                                             <  2:2544892/100‑1 (MQ=255)
 ggTGCCGAGTGCATCGGCACGGATATGAACTAGGACGATACCGTGCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACgg                                                                                            <  2:3559130/100‑1 (MQ=255)
  gTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGc                                                                                           <  1:3284897/100‑1 (MQ=255)
  gTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGc                                                                                           >  1:135477/1‑100 (MQ=255)
  gTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGc                                                                                           <  1:273487/100‑1 (MQ=255)
     ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                        <  1:355517/100‑1 (MQ=255)
     ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                        >  1:935944/1‑100 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       <  1:3243732/100‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       >  1:817708/1‑100 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       >  2:1375666/1‑100 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       >  1:2502878/1‑100 (MQ=255)
       gAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATggg                                                                                      >  1:1751822/1‑100 (MQ=255)
           gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  <  2:1046154/100‑1 (MQ=255)
            cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                 >  1:2539959/1‑100 (MQ=255)
            cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                                 >  1:2333615/1‑100 (MQ=255)
             ntCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                                >  2:1650146/2‑100 (MQ=255)
              tCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTa                                                                               <  2:2148136/100‑1 (MQ=255)
               cgccACGGATATGAACCATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                              <  2:2979662/97‑1 (MQ=255)
                gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             <  1:3373921/100‑1 (MQ=255)
                gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             <  1:1769304/100‑1 (MQ=255)
                gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  1:1146542/1‑100 (MQ=255)
                gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             <  1:2958630/100‑1 (MQ=255)
                 tCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacan                                                                            <  2:2725369/100‑2 (MQ=255)
                  cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                           <  1:517433/100‑1 (MQ=255)
                  cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                           >  2:3256918/1‑100 (MQ=255)
                   aCGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCa                                                                          >  2:464467/1‑100 (MQ=255)
                      gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaan                                                                       <  1:1370600/100‑2 (MQ=255)
                      gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaaa                                                                       >  1:586216/1‑100 (MQ=255)
                      gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaaa                                                                       <  2:586209/100‑1 (MQ=255)
                       atatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAAt                                                                      <  1:3198295/100‑1 (MQ=255)
                       atatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAAt                                                                      <  2:1930332/100‑1 (MQ=255)
                        tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                     <  2:2537325/100‑1 (MQ=255)
                         atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                    >  1:679777/1‑100 (MQ=255)
                         atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                    >  2:1416929/1‑100 (MQ=255)
                             aCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                <  2:3218161/100‑1 (MQ=255)
                                aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                             >  1:465123/1‑100 (MQ=255)
                                aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGGGGTCCAGGGCCCGGCATGGGCAGGTTACACCCAAATCCCAACCTg                                                             >  2:281540/1‑100 (MQ=255)
                                 tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                            >  2:1181182/1‑100 (MQ=255)
                                 gatctaTACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                            <  2:3080918/95‑1 (MQ=255)
                                  gACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCAAACCTGGa                                                           <  1:672455/100‑1 (MQ=255)
                                   aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                          >  2:888520/1‑100 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                         >  1:3384297/1‑100 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                         <  1:2802887/100‑1 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                         <  1:2725780/100‑1 (MQ=255)
                                     gATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCa                                                        <  1:1062466/100‑1 (MQ=255)
                                      aTACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGagatc                                                       >  1:3080942/1‑96 (MQ=255)
                                      aTACCGTTCGCCTCACCATGCAACACGCGCAGCAGCACGGCTGGGCAGTCGTACAGGGCACGGCATGCGCAGGTTACACCTAACTCCCAACCTGGGTCAt                                                       >  2:1879133/1‑100 (MQ=255)
                                       tACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGATGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATg                                                      >  1:1249063/1‑100 (MQ=255)
                                       tACCGTGCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATg                                                      <  2:1672305/100‑1 (MQ=255)
                                        aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                     <  2:3167148/100‑1 (MQ=255)
                                        aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                     <  1:3181536/100‑1 (MQ=255)
                                         ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACATGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                    >  1:3044835/1‑100 (MQ=255)
                                         ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                    >  2:3284334/1‑100 (MQ=255)
                                         ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                    >  2:3692924/1‑100 (MQ=255)
                                          cGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCaa                                                   >  2:567825/1‑100 (MQ=255)
                                           gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  <  1:606429/100‑1 (MQ=255)
                                           gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  <  2:1436889/100‑1 (MQ=255)
                                           gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                  <  1:3692952/100‑1 (MQ=255)
                                             tCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                                <  2:2213261/100‑1 (MQ=255)
                                             tCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                                <  1:3300105/100‑1 (MQ=255)
                                               gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                              >  1:2146781/1‑100 (MQ=255)
                                               gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                              >  2:200178/1‑100 (MQ=255)
                                                ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                             <  1:1029256/100‑1 (MQ=255)
                                                  tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                           >  2:334050/1‑100 (MQ=255)
                                                   gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGTAAGGTTACAACAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          <  1:2847641/100‑1 (MQ=255)
                                                   gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                          >  2:2711510/1‑100 (MQ=255)
                                                    aCCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCaa                                         >  1:367794/1‑100 (MQ=255)
                                                       aTGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAccc                                      <  1:3256944/100‑1 (MQ=255)
                                                        tGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                     >  2:3201153/1‑100 (MQ=255)
                                                        tGCAACACGCGCAGCAGCACGGCGGGGCAGTGGTCCAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGTTCATGCAAGGCTACGCAAccc                                      >  2:1959204/1‑99 (MQ=255)
                                                         gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCCCGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                    >  2:817813/1‑100 (MQ=255)
                                                         gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                    <  1:1445009/100‑1 (MQ=255)
                                                         gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                    <  1:301633/100‑1 (MQ=255)
                                                          cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTCCACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                   >  2:205747/1‑100 (MQ=255)
                                                          cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                   >  1:1916165/1‑100 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                  >  1:1058240/1‑100 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                  <  2:208728/100‑1 (MQ=255)
                                                            acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                 <  2:474944/100‑1 (MQ=255)
                                                            acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                 <  1:2949832/100‑1 (MQ=255)
                                                             cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                <  1:3467314/100‑1 (MQ=255)
                                                             cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                >  1:226508/1‑100 (MQ=255)
                                                             cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCACAACCCGGATCATGCAAGGCTACGCAACCCTgg                                   >  2:814769/1‑97 (MQ=255)
                                                               cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                              >  2:184497/1‑100 (MQ=255)
                                                               cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATACCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                              >  2:2681458/1‑100 (MQ=255)
                                                               cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                              >  2:1523823/1‑100 (MQ=255)
                                                                    agcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                         <  1:2400869/100‑1 (MQ=255)
                                                                     gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                        <  1:2980064/100‑1 (MQ=255)
                                                                       agcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGt                      <  1:3036330/100‑1 (MQ=255)
                                                                       agcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGt                      >  2:2962498/1‑100 (MQ=255)
                                                                       agcaCGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGt                      >  1:1422195/1‑100 (MQ=255)
                                                                         caCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                    <  1:1634291/100‑1 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  <  2:1017612/100‑1 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCACGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  <  2:3158354/100‑1 (MQ=255)
                                                                              cTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaa               >  2:2626178/1‑100 (MQ=255)
                                                                               tGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaaa              <  1:3441587/100‑1 (MQ=255)
                                                                                gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt             >  1:1203598/1‑100 (MQ=255)
                                                                                 ggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATg            >  2:3395589/1‑100 (MQ=255)
                                                                                        gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  2:2896131/1‑100 (MQ=255)
                                                                                         tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt    <  2:3671160/100‑1 (MQ=255)
                                                                                           cAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATgg  >  2:1143373/1‑100 (MQ=255)
                                                                                               |                                                                                               
CGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007980‑3008170

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