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,0750AG98.2% 366.5 / ‑2.9 109D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (57/50);  total (58/51)
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.52e-01

Only 100 of 109 total aligned reads displayed.

GGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGG  >  NC_000913/3007981‑3008171
                                                                                              |                                                                                                
ggTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACgg                                                                                             >  2:136975/1‑100 (MQ=255)
ggTGCCGAGTGCATCATCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACgg                                                                                             <  1:268521/100‑1 (MQ=255)
ggTGCCGAGGGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGTGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACgg                                                                                             <  2:811780/100‑1 (MQ=255)
   gCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCAt                                                                                          >  1:2198675/1‑100 (MQ=255)
    ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                         <  1:750267/100‑1 (MQ=255)
    ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                         >  2:2301907/1‑100 (MQ=255)
     cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                        >  2:749958/1‑100 (MQ=255)
     cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                        >  2:10855/1‑100 (MQ=255)
      gAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATggg                                                                                       >  1:816520/1‑100 (MQ=255)
      gAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATggg                                                                                       <  1:121137/100‑1 (MQ=255)
      gAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATggg                                                                                       >  1:1926440/1‑100 (MQ=255)
         tGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAg                                                                                    <  2:534015/100‑1 (MQ=255)
          gCATCGTCACGGATATGAACTATGACGATACTGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                   <  1:1923231/100‑1 (MQ=255)
          gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                   <  2:975834/100‑1 (MQ=255)
          gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                   >  2:20314/1‑100 (MQ=255)
          gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                   <  1:1588896/100‑1 (MQ=255)
          gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                   >  2:2406730/1‑100 (MQ=255)
            aTCGTCACGGATATGAACTATGACGATACCGTTCCCCCGACCATGCACCACGCACAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                                 <  2:2064492/100‑1 (MQ=255)
             tCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTa                                                                                >  1:1211353/1‑100 (MQ=255)
             tCGTCACGCATATGAACCATGACGATACCGTTCGCCTGACCATGCGACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGACGGtt                                                                                 <  1:2593689/99‑1 (MQ=255)
              cGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                               <  1:1325664/100‑1 (MQ=255)
               gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                              >  2:1928892/1‑100 (MQ=255)
                tCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAGGGTTacac                                                                             <  2:958978/100‑1 (MQ=255)
                tCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                             <  1:443594/100‑1 (MQ=255)
                tCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                             <  1:2549735/100‑1 (MQ=255)
                gcACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                             >  2:869311/2‑100 (MQ=255)
                  aCGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCa                                                                           <  2:587395/100‑1 (MQ=255)
                   cGGATATGAACTATGGCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGCGAAGGGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                          <  1:2450483/100‑1 (MQ=255)
                   cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                          >  1:1651035/1‑100 (MQ=255)
                   cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                          <  2:732564/100‑1 (MQ=255)
                     gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACAaaaa                                                                        >  2:1263801/1‑100 (MQ=255)
                       tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                      >  1:2012610/1‑100 (MQ=255)
                        atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                     <  2:1630399/100‑1 (MQ=255)
                        atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                     >  2:1380541/1‑100 (MQ=255)
                        atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                     <  2:947856/100‑1 (MQ=255)
                        atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATcc                                                                     >  1:1130416/1‑100 (MQ=255)
                        atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATcc                                                                     <  2:1130420/100‑1 (MQ=255)
                         tGAACTATGACGATCCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCCGGGCAGTGGTACAGGGCACGGCCTGGGAAGGTTACACCAAAATcac                                                                    >  1:2395908/1‑98 (MQ=255)
                         tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGGAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCa                                                                           <  1:371861/93‑1 (MQ=255)
                         tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                    >  1:1172214/1‑100 (MQ=255)
                          gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCa                                                                   <  2:2407636/100‑1 (MQ=255)
                            aCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                 >  1:2246167/1‑100 (MQ=255)
                             cTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAcc                                                                <  1:2335286/100‑1 (MQ=255)
                             cTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAcc                                                                >  1:1459511/1‑100 (MQ=255)
                                tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                             >  2:324962/1‑100 (MQ=255)
                                tGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTgg                                                             <  2:2476985/100‑1 (MQ=255)
                                 gACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGa                                                            >  1:2517574/1‑100 (MQ=255)
                                 gACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGa                                                            >  2:955698/1‑100 (MQ=255)
                                  aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCCCGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                           >  1:1173256/1‑100 (MQ=255)
                                   cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                          <  2:1139764/100‑1 (MQ=255)
                                       aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                      <  1:2454652/100‑1 (MQ=255)
                                       aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                      <  1:1223813/100‑1 (MQ=255)
                                        ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                     >  2:1632614/1‑100 (MQ=255)
                                        ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                     >  2:84429/1‑100 (MQ=255)
                                          gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                   >  2:1463126/1‑100 (MQ=255)
                                             cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                                >  2:1187872/1‑100 (MQ=255)
                                              gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                               >  1:1055378/1‑100 (MQ=255)
                                              gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                               >  1:567233/1‑100 (MQ=255)
                                              gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                               >  2:337241/1‑100 (MQ=255)
                                              gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                               >  1:2360806/1‑100 (MQ=255)
                                              gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                               <  2:2244950/100‑1 (MQ=255)
                                              ccgatctcATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                               <  2:1937555/93‑1 (MQ=255)
                                               ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                              >  1:1116161/1‑100 (MQ=255)
                                               ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                              >  2:1115800/1‑100 (MQ=255)
                                                 tGACCATGCAACACGCGGAGCAGCACGGCTGGGAAGAGGTACAGGGCAAGGCATGGGAAGGTTACAACAAAATCCCAACCTGGATCATGCAAGGCTACGc                                            >  1:2195474/1‑100 (MQ=255)
                                                 tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                            <  1:1187868/100‑1 (MQ=255)
                                                  gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                           >  2:1926881/1‑100 (MQ=255)
                                                    ccATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAc                                         >  2:2112668/1‑100 (MQ=255)
                                                    ccATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAc                                         >  2:2297127/1‑100 (MQ=255)
                                                     cATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTaagatcgg                                        >  1:1937545/1‑93 (MQ=255)
                                                     cATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAcc                                        >  2:656681/1‑100 (MQ=255)
                                                        gCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTg                                     <  2:175641/100‑1 (MQ=255)
                                                         cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                    >  2:1106660/1‑100 (MQ=255)
                                                           acacGCGCAGCAGCACGGCTGGGAAGTTGTACAGGGCACGGCATGGGAAGGTTACACCAAAGTCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                  <  1:1084024/100‑1 (MQ=255)
                                                           acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                  <  2:2069280/100‑1 (MQ=255)
                                                           acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGAGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAATGCTACGCAACCCTGGCa                                  <  2:1171286/100‑1 (MQ=255)
                                                            cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                  >  2:67982/1‑99 (MQ=255)
                                                            cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                 <  2:837768/100‑1 (MQ=255)
                                                            cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                 >  2:651993/1‑100 (MQ=255)
                                                               gcgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATg                              <  1:455046/100‑1 (MQ=255)
                                                                  cagcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAg                           <  2:709781/100‑1 (MQ=255)
                                                                  cagcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAg                           <  2:1655896/100‑1 (MQ=255)
                                                                    gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                         >  2:2558318/1‑100 (MQ=255)
                                                                     cagcaCGGCTGGGAAGTAGTACAGGGCCCGGCATTGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                          <  1:1101715/98‑1 (MQ=255)
                                                                      agcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCg                        >  2:250536/1‑99 (MQ=255)
                                                                       gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTc                      <  1:2123304/100‑1 (MQ=255)
                                                                         aCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGGTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                     <  1:2464141/99‑1 (MQ=255)
                                                                          cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                   <  2:712810/100‑1 (MQ=255)
                                                                          cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                   <  1:924894/100‑1 (MQ=255)
                                                                          cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGGTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                   >  2:764406/1‑100 (MQ=255)
                                                                           ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc                  <  2:2073634/100‑1 (MQ=255)
                                                                            gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                 <  2:1201253/100‑1 (MQ=255)
                                                                            gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                 <  1:816670/100‑1 (MQ=255)
                                                                               gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGTAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaaa               >  2:2369178/1‑99 (MQ=255)
                                                                                 gAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGc            >  2:2214114/1‑100 (MQ=255)
                                                                                 gAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGc            >  2:2213809/1‑100 (MQ=255)
                                                                                 gAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGc            >  2:85735/1‑100 (MQ=255)
                                                                                      ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa       >  1:689119/1‑100 (MQ=255)
                                                                                        tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt     >  2:46779/1‑100 (MQ=255)
                                                                                           aGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATggg  <  2:1921638/100‑1 (MQ=255)
                                                                                              |                                                                                                
GGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGG  >  NC_000913/3007981‑3008171

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: