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.4% 352.1 / 1.5 111D189G (GAC→GGC) ygeX2,3‑diaminopropionate ammonia lyase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base G (53/54);  total (56/55)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 111 total aligned reads displayed.

CGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGG  >  NC_000913/3007980‑3008171
                                                                                               |                                                                                                
cGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACg                                                                                              >  1:2107136/1‑100 (MQ=255)
    gCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCAt                                                                                          >  1:2107408/1‑100 (MQ=255)
    gCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCAt                                                                                          >  1:1206382/1‑100 (MQ=255)
     ccGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATg                                                                                         <  2:2024171/100‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCAt                                                                                          <  1:1635953/98‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                        <  2:1396895/100‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                        <  2:160629/100‑1 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                        >  1:560706/1‑100 (MQ=255)
      cGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATgg                                                                                        >  1:80999/1‑100 (MQ=255)
       gAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGACGTGGTACAGGGCACGGCATggg                                                                                       >  1:2458643/1‑100 (MQ=255)
           gCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAgg                                                                                   <  2:1495239/100‑1 (MQ=255)
             aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGTAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGACTGGGAGGGtt                                                                                 >  1:777344/1‑100 (MQ=255)
             aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                                 <  1:677318/100‑1 (MQ=255)
             aTCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGtt                                                                                 <  2:1645881/100‑1 (MQ=255)
                gTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTaca                                                                              >  1:2771738/1‑100 (MQ=255)
                 tCACGGATATGAACTATGACGATACCGTTCGCCTGCCCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTacac                                                                             >  1:2491655/1‑100 (MQ=255)
                    cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                          <  2:1613779/100‑1 (MQ=255)
                    cGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaa                                                                          >  1:923404/1‑100 (MQ=255)
                    cGGATATGAACTATGACGATACCGTACACCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGATACAGGGCACGGCAGGGGAAGGGTACACCa                                                                           >  1:2053406/1‑99 (MQ=255)
                      gATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCaaaa                                                                        <  1:1029769/100‑1 (MQ=255)
                        tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                      <  2:10582/100‑1 (MQ=255)
                        tatGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATc                                                                      >  1:1754273/1‑100 (MQ=255)
                         atGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATcc                                                                     <  1:1994007/100‑1 (MQ=255)
                          tGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATccc                                                                    >  1:2417580/1‑100 (MQ=255)
                          tGAACTATGACGATACAGTTCGCATGACCATGCAACAAGCGCAGCAGCACGGCTGGGCAGTGGTACAGGGAGCGGGATGGGAAGGTTAAAGCGAAAtaca                                                                    >  1:1679679/1‑97 (MQ=255)
                           gAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCa                                                                   >  1:983412/1‑100 (MQ=255)
                            aaCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCaa                                                                  <  2:2791570/100‑1 (MQ=255)
                             aCTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAAc                                                                 <  2:2501649/100‑1 (MQ=255)
                               tATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCt                                                               <  2:1206384/100‑1 (MQ=255)
                                aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                              <  2:2236671/100‑1 (MQ=255)
                                aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                              >  1:435413/1‑100 (MQ=255)
                                aTGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTg                                                              <  1:1797707/100‑1 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                          >  2:116415/1‑100 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                          >  1:65743/1‑100 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                          <  1:2509388/100‑1 (MQ=255)
                                    cGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATc                                                          <  2:435415/100‑1 (MQ=255)
                                     gATACCGTTCGCCTGCCCATGCACCACGCACAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAGGGTTACACCAAAATCCCAACCTGGATCa                                                         <  1:314745/100‑1 (MQ=255)
                                     gATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCa                                                         <  1:2395945/100‑1 (MQ=255)
                                     gATACCGTTCGCCTGAACATGCAACACAAGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACAGCATGGGAAGGTTACACCAAAATCCTAACATGGATCa                                                         >  1:2201366/1‑100 (MQ=255)
                                         ccGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                     >  1:1468071/1‑100 (MQ=255)
                                           gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                   >  2:2028417/1‑100 (MQ=255)
                                           gTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                   <  2:1353148/100‑1 (MQ=255)
                                           gTTCGCCTGACCATGCAACACGAGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCa                                                     >  1:2195306/1‑98 (MQ=255)
                                           gTTAGCCTGACCATGCAACAGGTGCAGCAGCACGGCTGGTAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAg                                                   <  1:2605854/100‑1 (MQ=255)
                                            ttCGCCTGACCACGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAgg                                                  >  2:1916298/1‑100 (MQ=255)
                                              cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                                >  2:1120869/1‑100 (MQ=255)
                                              cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                                >  2:2002561/1‑100 (MQ=255)
                                              cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCt                                                >  1:1897235/1‑100 (MQ=255)
                                              cGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCAGGCAAGGCt                                                >  1:2512100/1‑100 (MQ=255)
                                               tcttGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTa                                               <  1:1375230/97‑1 (MQ=255)
                                                  tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTaaga                                            >  2:1375232/1‑97 (MQ=255)
                                                  tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                            >  1:2065350/1‑100 (MQ=255)
                                                  tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                            <  1:1269416/100‑1 (MQ=255)
                                                  tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                            <  2:1897243/100‑1 (MQ=255)
                                                  tGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGAGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGc                                            >  1:410201/1‑100 (MQ=255)
                                                   gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                           >  1:1705851/1‑100 (MQ=255)
                                                   gACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCa                                           >  2:2503183/1‑100 (MQ=255)
                                                    aCCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCaa                                          >  2:2693014/1‑100 (MQ=255)
                                                      cATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGACATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAAcc                                        >  2:426389/1‑100 (MQ=255)
                                                          cAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCt                                      >  1:2735951/1‑98 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                   <  2:795660/100‑1 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                   <  1:114184/100‑1 (MQ=255)
                                                           aaCACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGc                                   >  1:596832/1‑100 (MQ=255)
                                                            acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGcn                                  <  1:1203459/100‑2 (MQ=255)
                                                            acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                  >  2:2836413/1‑100 (MQ=255)
                                                            acacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                  <  1:289771/100‑1 (MQ=255)
                                                             ctcgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                 <  1:901495/98‑1 (MQ=255)
                                                             cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCa                                  <  1:1256132/99‑1 (MQ=255)
                                                             cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                 >  1:1841205/1‑100 (MQ=255)
                                                             cacGCGCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAg                                 >  1:2373885/1‑100 (MQ=255)
                                                               cgcgCAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGAt                               >  2:666444/1‑100 (MQ=255)
                                                                 cgcAGCAGCACGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATACCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGa                             >  2:2724163/1‑100 (MQ=255)
                                                                   cagcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAg                           >  1:38592/1‑100 (MQ=255)
                                                                   cagcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAg                           <  2:2619979/100‑1 (MQ=255)
                                                                    agcagcaCGGCTGGGAAGTGGTACAGGGCATGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGaa                            <  1:1470011/98‑1 (MQ=255)
                                                                     gcagcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGcc                         <  1:767101/100‑1 (MQ=255)
                                                                        gcaCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTc                      >  1:2820292/1‑100 (MQ=255)
                                                                         caCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCg                     >  1:1576594/1‑100 (MQ=255)
                                                                          aCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTTGATCATGCAAGGCTACGTAACCCTGGCAGATGAAGCCGTCg                     <  1:876761/99‑1 (MQ=255)
                                                                          aCGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGa                    <  1:2169574/100‑1 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGa                    <  1:1668080/99‑1 (MQ=255)
                                                                           cGGCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAg                   >  1:1757050/1‑100 (MQ=255)
                                                                            ggCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGc                  <  2:1409255/100‑1 (MQ=255)
                                                                             gCTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCa                 <  1:688826/100‑1 (MQ=255)
                                                                              cTGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaa                >  2:1831491/1‑100 (MQ=255)
                                                                               tGGGAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCaaa               <  2:649838/100‑1 (MQ=255)
                                                                                gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt              >  2:2038895/1‑100 (MQ=255)
                                                                                gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt              >  1:2193692/1‑100 (MQ=255)
                                                                                gggAAGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAAt              >  2:603565/1‑100 (MQ=255)
                                                                                   aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg           <  2:2107415/100‑1 (MQ=255)
                                                                                   aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg           <  1:2392869/100‑1 (MQ=255)
                                                                                   aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg           <  1:323511/100‑1 (MQ=255)
                                                                                   aaGTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCg           <  2:591175/100‑1 (MQ=255)
                                                                                     gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg         <  2:1567815/100‑1 (MQ=255)
                                                                                     gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg         >  2:496785/1‑100 (MQ=255)
                                                                                     gTGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTg         <  1:666443/100‑1 (MQ=255)
                                                                                      tGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGtta        <  2:1972075/100‑3 (MQ=255)
                                                                                      tGGTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGa        >  2:2295415/1‑100 (MQ=255)
                                                                                       ggTACAGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGaa       >  2:113413/1‑100 (MQ=255)
                                                                                            aGGGCACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATggg  <  2:642926/100‑1 (MQ=255)
                                                                                               |                                                                                                
CGGTGCCGAGTGCATCGTCACGGATATGAACTATGACGATACCGTTCGCCTGACCATGCAACACGCGCAGCAGCACGGCTGGGAAGTGGTACAGGACACGGCATGGGAAGGTTACACCAAAATCCCAACCTGGATCATGCAAGGCTACGCAACCCTGGCAGATGAAGCCGTCGAGCAAATGCGTGAAATGGG  >  NC_000913/3007980‑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: