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,0750AG96.7% 286.8 / ‑1.5 89D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (41/43);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

CCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAT  >  NC_000913/3007985‑3008168
                                                                                          |                                                                                             
ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                      <  1:571681/100‑1 (MQ=255)
 cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                     <  2:2139543/100‑1 (MQ=255)
     tGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAg                                                                                 <  2:63615/100‑1 (MQ=255)
      gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                <  1:239328/100‑1 (MQ=255)
       cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                               <  1:931933/100‑1 (MQ=255)
       cATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGt                                                                               >  2:1928437/1‑100 (MQ=255)
        aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                              <  1:1086349/100‑1 (MQ=255)
        aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGtt                                                                              >  2:2169882/1‑100 (MQ=255)
           gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                           <  1:1725433/100‑1 (MQ=255)
            tCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                          <  2:1075998/100‑1 (MQ=255)
            tCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                          >  2:1799911/1‑100 (MQ=255)
             cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                         <  1:114544/100‑1 (MQ=255)
             cACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACAcc                                                                         <  1:674124/100‑1 (MQ=255)
              aCGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCa                                                                        <  2:1774509/100‑1 (MQ=255)
               cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                       >  1:841375/1‑100 (MQ=255)
               cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                       >  2:1951193/1‑100 (MQ=255)
                 gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaaa                                                                     >  1:1484345/1‑100 (MQ=255)
                   tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                   >  1:464174/1‑100 (MQ=255)
                   tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCCGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                   >  2:331003/1‑100 (MQ=255)
                    atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                  <  2:288825/100‑1 (MQ=255)
                     tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                 >  2:1185729/1‑100 (MQ=255)
                     tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                  >  1:151275/1‑99 (MQ=255)
                     tGAACTATGACGATACAGTTCGCATGAACATGCAACACGCGCAGCAGCACGGCTCGGAAGTGGTAAAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                 >  2:550143/1‑100 (MQ=255)
                     tGAAATATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                 <  2:796529/100‑1 (MQ=255)
                       aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                               >  1:941448/1‑100 (MQ=255)
                       aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                               >  1:940305/1‑100 (MQ=255)
                       aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                               >  2:1480204/1‑100 (MQ=255)
                           aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                           >  1:1449006/1‑100 (MQ=255)
                             gACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGa                                                         >  1:337003/1‑100 (MQ=255)
                              aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                        <  2:370118/100‑1 (MQ=255)
                              aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                        >  1:942799/1‑100 (MQ=255)
                              aCGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGAt                                                        <  1:1225276/100‑1 (MQ=255)
                               cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                       <  2:1244190/100‑1 (MQ=255)
                               cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                       <  1:369872/100‑1 (MQ=255)
                               cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                       <  1:429041/100‑1 (MQ=255)
                               cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGAACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                       <  1:2169897/100‑1 (MQ=255)
                                 aTACCCTTCGCCTGACCATGCAACACGCGCACGATCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAATGTTACACCAAAATCCCAACCTGGATc                                                       <  2:223616/98‑1 (MQ=255)
                                  tACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATg                                                    >  1:1955722/1‑100 (MQ=255)
                                   aCCGTTCGCCTGACCATGCAACACCCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGc                                                   <  1:1980657/100‑1 (MQ=255)
                                    ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                  >  2:441545/1‑100 (MQ=255)
                                    ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                  >  2:539721/1‑100 (MQ=255)
                                    ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                  <  1:476982/100‑1 (MQ=255)
                                    ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                  >  2:1210370/1‑100 (MQ=255)
                                     cGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCaa                                                 >  1:419910/1‑100 (MQ=255)
                                        tCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGc                                              <  2:2248489/100‑1 (MQ=255)
                                         cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                             <  1:878227/100‑1 (MQ=255)
                                         cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                             >  1:37717/1‑100 (MQ=255)
                                             tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                         <  1:368884/100‑1 (MQ=255)
                                             tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                         >  1:2211177/1‑100 (MQ=255)
                                             tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                         >  1:1826577/1‑100 (MQ=255)
                                                ccATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAc                                      >  2:1273830/1‑100 (MQ=255)
                                                  aTGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAccc                                    >  2:1948160/1‑100 (MQ=255)
                                                  aTGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAccc                                    <  1:916289/100‑1 (MQ=255)
                                                  aTGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAccc                                    >  2:1014484/1‑100 (MQ=255)
                                                  aTGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAccc                                    <  1:2239557/100‑1 (MQ=255)
                                                   tGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                   >  1:538690/1‑100 (MQ=255)
                                                   tGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                   >  2:1644659/1‑100 (MQ=255)
                                                     cAACACGCGCAGCAGCACGGCAGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                   >  2:8427/1‑98 (MQ=255)
                                                      aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                >  1:2047501/1‑100 (MQ=255)
                                                       acacGCGCAGCAGCGCGGCTGGGAAGTGGTACAGGTCACGGCATGGGAAGGTTACCCCAAAATCCCAACCTTGTTCATGTAAGGCTACGCAACCCTGGCa                               <  2:714778/100‑1 (MQ=255)
                                                       acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                               <  2:1076287/100‑1 (MQ=255)
                                                       acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                               >  2:2065517/1‑100 (MQ=255)
                                                       acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                               >  1:789898/1‑100 (MQ=255)
                                                        cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGGATGGGAAGGGTACACCAACATCCCACCCTGGATCATGCAAGGATACGCAACCCTGGCAg                              <  2:660323/100‑1 (MQ=255)
                                                        cacGCGCAGCAGCACGGCTGGGAAGTAGTACAGGGCACGGCATGTGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                              <  2:456791/100‑1 (MQ=255)
                                                         acGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCGAACCTGGATCACGCAAGGCTACGCAACCCt                                   >  2:1159322/1‑94 (MQ=255)
                                                          cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                            >  1:858226/1‑100 (MQ=255)
                                                           gcgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATg                           <  2:1535462/100‑1 (MQ=255)
                                                            cgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATg                           <  1:1799923/99‑1 (MQ=255)
                                                               agcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAg                        <  1:1156662/99‑1 (MQ=255)
                                                               agcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGc                       <  2:1955709/100‑1 (MQ=255)
                                                                  agcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGt                    >  1:1841548/1‑100 (MQ=255)
                                                                   gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCACGGCTACGCAACCCTGGCAGATGAAGcc                      >  2:1696138/1‑97 (MQ=255)
                                                                   gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTc                   <  1:2051449/100‑1 (MQ=255)
                                                                    caCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                  <  1:118513/100‑1 (MQ=255)
                                                                    caCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                  <  1:1849569/100‑1 (MQ=255)
                                                                     aCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGa                 >  1:1396368/1‑100 (MQ=255)
                                                                       ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc               <  1:1644670/100‑1 (MQ=255)
                                                                        gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc               >  2:1291497/1‑99 (MQ=255)
                                                                            ggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATg          <  2:684327/100‑1 (MQ=255)
                                                                             gAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGc         <  1:1692810/100‑1 (MQ=255)
                                                                               aGTGGTACAGGACACGGCATAGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGACAGACAAAGCCGTCAAGCAAATGCGt       >  2:736628/1‑100 (MQ=255)
                                                                                gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg      <  1:1529701/100‑1 (MQ=255)
                                                                                gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg      <  2:2197375/100‑1 (MQ=255)
                                                                                gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg      <  1:464781/100‑1 (MQ=255)
                                                                                  ggTACAGGGCACGGCATGGGAAGGTTACACCAGAATCTCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGACCAAATGCATGaa    >  1:2132439/1‑100 (MQ=255)
                                                                                   gTACAGGTCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaaa   >  2:160330/1‑100 (MQ=255)
                                                                                   gTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa    <  1:2021959/99‑1 (MQ=255)
                                                                                     aCAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGCTGAAGCCGTCGAGCAAATGCGTGAAAt  >  1:1173207/1‑99 (MQ=255)
                                                                                          |                                                                                             
CCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAAT  >  NC_000913/3007985‑3008168

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: