Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,215,178 Δ1 bp 100% coding (899/930 nt) metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,215,1780A.98.7% 677.4 / 1.2 149coding (899/930 nt)metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (80/67);  total (81/68)
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.85e-01

Only 100 of 149 total aligned reads displayed.

TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTT  >  NC_000913/4215083‑4215272
                                                                                               |                                                                                              
tCCGCAAAATTCACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAAATATTACGTCTACCAGATCACGCCATACG‑TCTAc                                                                                           <  1:2493258/100‑1 (MQ=255)
tCCGCAAAATACACAGCGAGCGAGCTGCCGTAGTCACGGTAATTTACTGTTTACCAACTTGCGCAACTATTACGTCTACCATATCATGCCCTACG‑TCTAc                                                                                           >  1:547924/1‑100 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTAc                                                                                           >  1:1231411/1‑99 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACgg                                                                                         >  1:409972/1‑100 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                        >  2:2069823/1‑100 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                        <  2:2092283/100‑1 (MQ=255)
     aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcac                                                                                      >  1:3025102/1‑100 (MQ=255)
     aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcac                                                                                      >  2:1079132/1‑100 (MQ=255)
        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATg                                                                                   <  1:1292961/100‑1 (MQ=255)
          acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGaa                                                                                 <  2:55775/100‑1 (MQ=255)
           cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAAt                                                                                >  2:2202269/1‑100 (MQ=255)
            acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATc                                                                               <  2:1178741/100‑1 (MQ=255)
             ncgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                              >  1:70774/2‑100 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                              <  1:3530761/100‑1 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                              >  2:2391649/1‑100 (MQ=255)
              cgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATcc                                                                              >  2:1384726/1‑100 (MQ=255)
              cgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCa                                                                             <  1:2628553/100‑1 (MQ=255)
               gcgagcgaGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                              <  1:2873423/98‑1 (MQ=255)
               gcgagcgaGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCaa                                                                            >  1:2781565/1‑100 (MQ=255)
                cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAAc                                                                           >  1:3139008/1‑100 (MQ=255)
                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGc                                                                         <  1:3606840/100‑1 (MQ=255)
                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGc                                                                         <  2:3194598/100‑1 (MQ=255)
                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGc                                                                         <  2:2550625/100‑1 (MQ=255)
                  agagagcTGGCGTAGTCACGGTAATTTACTGTTTACTAACTGACTGAACTATTACTTCTCCCAGCTCGCGCCATACG‑TATACGGCACATGAATCCAACGc                                                                         >  1:1828279/4‑100 (MQ=255)
                    cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTg                                                                       >  2:2891001/1‑100 (MQ=255)
                    cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTg                                                                       >  2:1559152/1‑100 (MQ=255)
                    ccagcTGCCGTAGTCACGGTAATTTACTGTTGACCAACTGTCTCAACTGTCACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCaa                                                                            <  1:3419298/93‑1 (MQ=255)
                     gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTgg                                                                      <  2:1175968/100‑1 (MQ=255)
                      agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTg                                                                       >  1:2154510/1‑98 (MQ=255)
                      agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGa                                                                     >  1:3356444/1‑100 (MQ=255)
                       gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGAt                                                                    >  1:1276631/1‑100 (MQ=255)
                       gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTAGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGAt                                                                    <  2:2669579/100‑1 (MQ=255)
                        cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt                                                                    <  2:1978485/100‑1 (MQ=255)
                         tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTa                                                                  <  2:2238383/100‑1 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTaa                                                                 <  2:709266/100‑1 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTaa                                                                 <  1:2898800/100‑1 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTaa                                                                 <  2:3800094/100‑1 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               >  1:1390771/1‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               >  2:2478188/1‑100 (MQ=255)
                              tAGTCACGGTACTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                             >  2:3176796/1‑100 (MQ=255)
                              tAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                             >  1:459334/1‑100 (MQ=255)
                               aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                            >  2:1086427/1‑100 (MQ=255)
                               aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                             >  1:1939341/1‑99 (MQ=255)
                                gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                           >  2:3241390/1‑100 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                          <  1:10405/100‑1 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                          <  1:3694994/100‑1 (MQ=255)
                                  gagtgtcggttACTGTTTACCAACTGGCTTGACTTTTACGGTGTCGATAGAACGGCATGCG‑TCTACGGCATAGTAAGCGAACGGTGGATTAATCTTCtgt                                                         <  1:3370088/91‑1 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                         <  1:496128/100‑1 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                         >  1:2415061/1‑100 (MQ=255)
                                    ctgTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                       <  2:1527349/98‑1 (MQ=255)
                                    cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                       >  2:1257191/1‑100 (MQ=255)
                                    cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACAGCACATGAATCAAACGCTGGATTAATCTTCTGTGa                                                       >  1:2845620/1‑100 (MQ=255)
                                     ggTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                      >  2:1930232/1‑100 (MQ=255)
                                     ggTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                      <  1:2293009/100‑1 (MQ=255)
                                     ggTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCACATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                        <  1:2885709/98‑1 (MQ=255)
                                      gTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTTATCATCTGTGa                                                       <  1:1989679/98‑1 (MQ=255)
                                       tAATTTACTGTTTACCAACTGGCTCAACGATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                    <  2:160016/100‑1 (MQ=255)
                                        aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGt                                                   <  2:2546027/100‑1 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                 >  1:1081158/1‑100 (MQ=255)
                                            tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                               <  1:1527629/100‑1 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                             <  2:62122/100‑1 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                             <  2:353735/100‑1 (MQ=255)
                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATGCG‑TCTACGTCACATGAATCTAACGCGGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                        >  1:2065000/1‑100 (MQ=255)
                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                        <  2:824780/100‑1 (MQ=255)
                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                        >  1:104631/1‑100 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                       >  2:106182/1‑100 (MQ=255)
                                                      cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                     >  2:413167/1‑100 (MQ=255)
                                                      cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                     >  1:318448/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                  >  1:1165901/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                  <  2:790341/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                              <  2:1759123/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                              >  1:3284827/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATAACGCCATACG‑TCTACGGCACATGAATCCACCGCAGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGCTTCAGCa                              >  1:378258/1‑100 (MQ=255)
                                                              tCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                             <  1:3479801/100‑1 (MQ=255)
                                                               cAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            <  1:3192074/100‑1 (MQ=255)
                                                                aaCTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            <  2:3369591/99‑1 (MQ=255)
                                                                aaCTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCt                           >  2:3540191/1‑100 (MQ=255)
                                                                  cTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTa                         >  2:3434227/1‑100 (MQ=255)
                                                                   tATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAAAGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                        >  2:875050/1‑100 (MQ=255)
                                                                    aTTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAcc                       <  1:264960/100‑1 (MQ=255)
                                                                    aTTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAcc                       >  2:3760556/1‑100 (MQ=255)
                                                                    aTTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAcc                       >  2:3030783/1‑100 (MQ=255)
                                                                     ttACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCt                      >  1:3740602/1‑100 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                   >  1:2891837/1‑100 (MQ=255)
                                                                         gTCTACCAGATCACGCCATACG‑TCTACGGCAGATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                  <  2:363835/100‑1 (MQ=255)
                                                                         gTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                  >  2:1018205/1‑100 (MQ=255)
                                                                         gTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                  >  2:934490/1‑100 (MQ=255)
                                                                         gTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                  <  2:3218681/100‑1 (MQ=255)
                                                                           cTACAAGATCACGCCATTCG‑TCTACGTCCCATGAATCCAACGCCGGATTAATATTCTGTGTTATTCGCTCGTTTAGCGATTCAGCACCTTACCTCGGGCa                >  1:1522667/1‑100 (MQ=255)
                                                                             aCCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  2:3512705/1‑100 (MQ=255)
                                                                              ccAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  2:3688606/1‑99 (MQ=255)
                                                                              ccAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCt             >  2:741896/1‑100 (MQ=255)
                                                                                aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc           >  2:959816/1‑100 (MQ=255)
                                                                                aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc           >  2:3094014/1‑100 (MQ=255)
                                                                                   tCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg        <  1:1185623/100‑1 (MQ=255)
                                                                                   tCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg        >  1:700923/1‑100 (MQ=255)
                                                                                     aCGCCATACG‑TCTACGGCACATGAATCCCACGCTGCATTAATCTTCTGTCATAATCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg      >  1:912073/1‑100 (MQ=255)
                                                                                     aCGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg      >  2:3316496/1‑100 (MQ=255)
                                                                                     aCGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg      <  2:1626302/100‑1 (MQ=255)
                                                                                         cATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtt  <  2:1733281/100‑1 (MQ=255)
                                                                                               |                                                                                              
TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTT  >  NC_000913/4215083‑4215272

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: ‑