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.2% 432.7 / ‑4.6 124D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (65/58);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.02e-01

Only 100 of 125 total aligned reads displayed.

TCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007979‑3008170
                                                                                                |                                                                                               
tCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGgcac                                                                                              <  1:10552/100‑1 (MQ=255)
  ggTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACgg                                                                                            >  1:1125610/1‑100 (MQ=255)
    tGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCa                                                                                          >  2:142493/1‑100 (MQ=255)
     gCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCAt                                                                                         >  2:1557299/1‑100 (MQ=255)
       cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       <  2:1879337/100‑1 (MQ=255)
       cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                       <  1:2189795/100‑1 (MQ=255)
        gAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATggg                                                                                      >  1:1084477/1‑100 (MQ=255)
            gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  <  2:1222400/100‑1 (MQ=255)
            gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  <  2:1611732/100‑1 (MQ=255)
             cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGt                                                                                 <  1:1205994/100‑1 (MQ=255)
              aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                  <  1:268228/98‑1 (MQ=255)
                cGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                              <  1:811421/100‑1 (MQ=255)
                cGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                              <  2:2102100/100‑1 (MQ=255)
                 ttttcGGATATGAACTCTGACGATACCGTTCGCCTGCCCATACAACACGCGCAGCATCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  2:1006008/5‑100 (MQ=255)
                 gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             >  2:350679/1‑100 (MQ=255)
                 gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTac                                                                              <  2:52719/99‑1 (MQ=255)
                 gTCACGGATATGAACTAGGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                             <  2:115716/100‑1 (MQ=255)
                  tCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                            <  1:1329639/100‑1 (MQ=255)
                   ctcGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                           <  1:1627515/98‑1 (MQ=255)
                   cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                           <  1:947150/100‑1 (MQ=255)
                     cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCag                                                                         >  2:1627506/1‑99 (MQ=255)
                     cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                         <  2:2256472/100‑1 (MQ=255)
                     cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                         >  1:400807/1‑100 (MQ=255)
                        atgtaggcTATGACGATACCGTTCGCCTAAGCATGCAACACGCGCAGCAGCACGGCTCTGAAGTAGTACAGGGCACGGCATGGGAAGGTTACACCAAAAt                                                                      <  2:202440/93‑1 (MQ=255)
                        atatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAAt                                                                      >  1:759229/1‑100 (MQ=255)
                         tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                     >  1:2059282/1‑100 (MQ=255)
                         tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                     >  1:1394968/1‑100 (MQ=255)
                          atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                    <  2:2295792/100‑1 (MQ=255)
                          atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                    <  1:1473913/100‑1 (MQ=255)
                           tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGGAGTGGTACAGGGCACGGCCTGGGAAGGTTAGACCAAAATcc                                                                    >  2:2120886/1‑99 (MQ=255)
                           tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  1:2534577/1‑100 (MQ=255)
                           tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                   >  1:2570/1‑100 (MQ=255)
                             aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                 <  2:2505188/100‑1 (MQ=255)
                              aCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                <  2:2537507/100‑1 (MQ=255)
                              aCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                <  1:998558/100‑1 (MQ=255)
                               cTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAcc                                                               <  1:2153285/100‑1 (MQ=255)
                               cTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAcc                                                               >  1:92528/1‑100 (MQ=255)
                               cTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAcc                                                               >  1:2239078/1‑100 (MQ=255)
                                tATGACGATACCGTTCGCTTGACCATGCAACACGCGTAGCAGCACGGCAGGGAGGTGGTACAGGGCACGGCATGAGCAGGTTACACCACAATcgc                                                                   >  2:1110368/1‑93 (MQ=255)
                                 aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGGACAGGGCACGGCCTGGGAAGGTTACACCAAAATCCCAACCTg                                                             >  2:1905366/1‑100 (MQ=255)
                                     cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                         >  1:874925/1‑100 (MQ=255)
                                       aTACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCa                                                        <  2:1513263/99‑1 (MQ=255)
                                        tACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATg                                                      <  2:1970591/100‑1 (MQ=255)
                                         aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATg                                                      <  1:2546456/99‑1 (MQ=255)
                                         aCCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                     >  2:1891201/1‑100 (MQ=255)
                                          ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                    <  2:1767188/100‑1 (MQ=255)
                                          ccGTACGCCTAACCAAGCAACACGCGCAGCAGCACGGCTGGGAATTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCAt                                                       <  2:2262269/97‑1 (MQ=255)
                                           cGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAGGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCaa                                                   <  2:1827788/100‑1 (MQ=255)
                                           cGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCaa                                                   >  1:2232294/1‑100 (MQ=255)
                                             ttCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCGCGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAgg                                                 <  2:1441585/100‑1 (MQ=255)
                                              tCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                                <  2:552050/100‑1 (MQ=255)
                                               cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                               >  2:1327258/1‑100 (MQ=255)
                                               cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                               >  2:120034/1‑100 (MQ=255)
                                               cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCACAACCTGGATAATGCAAGGCt                                               >  2:2292994/1‑100 (MQ=255)
                                                gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                              <  2:748568/100‑1 (MQ=255)
                                                gCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                              >  1:404725/1‑100 (MQ=255)
                                                 ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                             >  1:2411117/1‑100 (MQ=255)
                                                 ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTAc                                             >  1:2411185/1‑100 (MQ=255)
                                                 ccTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGCAAGGTTACACCAAAATCCCAACCTGCATCATGCAAg                                                  >  2:2080922/1‑95 (MQ=255)
                                                   tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                           <  2:199619/100‑1 (MQ=255)
                                                      ccATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACCCCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAc                                        >  2:474491/1‑100 (MQ=255)
                                                       cATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAcc                                       >  2:799279/1‑100 (MQ=255)
                                                       cATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAcc                                       >  2:2083643/1‑100 (MQ=255)
                                                        aTGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAccc                                      >  1:2137176/1‑100 (MQ=255)
                                                         tGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                     <  2:874920/100‑1 (MQ=255)
                                                         tGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                     >  2:326306/1‑100 (MQ=255)
                                                         tGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                     >  1:245588/1‑100 (MQ=255)
                                                           cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                   >  1:1153119/1‑100 (MQ=255)
                                                           cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                   >  2:2153594/1‑100 (MQ=255)
                                                           cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTgg                                   >  1:2273919/1‑100 (MQ=255)
                                                            aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                  <  2:2239065/100‑1 (MQ=255)
                                                            aaCACGCGCAGCAGCAAGGCTGGGAAGTGGTACAGGGAACGGCATGGGAAGGTTACACCAAATTCCCAACCTGCATCATGCAAGGCTACGCAACCCt                                     <  2:722611/97‑1 (MQ=255)
                                                                cgcgCAGGCGAAAGGCTCGGAGGAGGTACAGGGCACGGAAGGGGAAGGGTACACCAAAAACCCACGCTGGATCATGCAAGGCTACGCAACACTGGCAGAt                              <  2:1557599/100‑1 (MQ=255)
                                                                cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                              <  1:1431132/100‑1 (MQ=255)
                                                                   gcagcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGCTCATGCAAGGCTACGCAACCCTGGCAGATGaa                           <  1:2468719/100‑1 (MQ=255)
                                                                    cagcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAg                          <  2:1223833/100‑1 (MQ=255)
                                                                      gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                        >  1:1179840/1‑100 (MQ=255)
                                                                         gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTc                     >  1:53153/1‑100 (MQ=255)
                                                                          caCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                    <  2:1426629/100‑1 (MQ=255)
                                                                            cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  >  1:2145990/1‑100 (MQ=255)
                                                                            cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  <  1:327602/100‑1 (MQ=255)
                                                                            cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                  >  2:1205936/1‑100 (MQ=255)
                                                                              gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                <  1:756199/100‑1 (MQ=255)
                                                                                 gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt             >  1:2385538/1‑100 (MQ=255)
                                                                                 gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt             <  1:1615759/100‑1 (MQ=255)
                                                                                  agAAGGGGTACAGTGCACGGCATGGGAAGGTTAAAACAAAATCCCAACATGGATCATGCAGGGCTACGCCACCATGGCAGATGGAGACGTCGAGCaca              >  2:2300829/2‑96 (MQ=255)
                                                                                      gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg        <  2:1326056/100‑1 (MQ=255)
                                                                                      gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg        <  1:50662/100‑1 (MQ=255)
                                                                                       tGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGa       >  1:294230/1‑100 (MQ=255)
                                                                                        ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa      <  2:418772/100‑1 (MQ=255)
                                                                                        ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGa       >  1:320451/1‑99 (MQ=255)
                                                                                         gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  1:1883121/1‑100 (MQ=255)
                                                                                         gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  2:1645257/1‑100 (MQ=255)
                                                                                         gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  2:389761/1‑100 (MQ=255)
                                                                                         gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  1:1099702/1‑100 (MQ=255)
                                                                                         gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa     >  1:1098955/1‑100 (MQ=255)
                                                                                          tACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAt    >  2:200012/1‑100 (MQ=255)
                                                                                            cAGGGCACGGCATGGGAAGGTTACACCAAAATCGCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGCTGAAGCCGTCGAGCAAATGCGTGAAATgg  >  2:1565242/1‑100 (MQ=255)
                                                                                            cAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATgg  >  1:175573/1‑100 (MQ=255)
                                                                                            cAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATgg  >  2:1784732/1‑100 (MQ=255)
                                                                                                |                                                                                               
TCGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGG  >  NC_000913/3007979‑3008170

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: