Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,215,178 Δ1 bp 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.97.2% 464.5 / 6.5 106coding (899/930 nt)metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (62/41);  total (64/42)
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.65e-01

Only 100 of 106 total aligned reads displayed.

CCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCT  >  NC_000913/4215084‑4215271
                                                                                              |                                                                                             
ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                         >  2:2308637/1‑100 (MQ=255)
ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                         >  1:1446252/1‑100 (MQ=255)
  gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                       <  1:2951397/100‑1 (MQ=255)
   cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGca                                                                                      >  2:205121/1‑100 (MQ=255)
    aaaaTACACCGCGAGCTAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcac                                                                                     <  2:1743528/100‑1 (MQ=255)
         acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGa                                                                                 >  1:2045151/1‑99 (MQ=255)
           acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATc                                                                              <  1:449941/100‑1 (MQ=255)
            ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                             <  2:2904601/100‑1 (MQ=255)
             cgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCa                                                                            >  1:1375003/1‑100 (MQ=255)
               cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAAc                                                                          >  1:1916970/1‑100 (MQ=255)
                   cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTg                                                                      <  2:769325/100‑1 (MQ=255)
                    gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTgg                                                                     <  2:1274728/100‑1 (MQ=255)
                    gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTg                                                                      >  2:579445/1‑99 (MQ=255)
                     cgatctcGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGa                                                                    <  2:2095408/94‑1 (MQ=255)
                      gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGAt                                                                   <  2:2194453/100‑1 (MQ=255)
                        tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTa                                                                 <  2:591370/100‑1 (MQ=255)
                         ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTaa                                                                >  1:791115/1‑100 (MQ=255)
                         ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTaa                                                                <  1:1501163/100‑1 (MQ=255)
                         ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTaa                                                                >  1:1660439/1‑100 (MQ=255)
                          gCGTAGTCACGGTAATTTACTGTTTACGAAATGGCTAAACTCTTAGGTCTACCAGATTACGCCATAAG‑TGTGCGGCACATGAATCTTACGCAgg                                                                     >  2:2904144/1‑94 (MQ=255)
                          gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAt                                                               >  1:1233212/1‑100 (MQ=255)
                           cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCAGACG‑TCTACGGCACATGAATCCAACGCTGGaagatcg                                                              >  1:2095420/1‑94 (MQ=255)
                             tAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                            >  2:1151681/1‑100 (MQ=255)
                              aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                           >  2:2406768/1‑100 (MQ=255)
                               gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                          >  1:2004439/1‑100 (MQ=255)
                                tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:1675391/100‑1 (MQ=255)
                                tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  2:791112/100‑1 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  2:2557452/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:259795/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:85888/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:1316820/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:1316448/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:319917/1‑100 (MQ=255)
                                  aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                       <  1:2046059/100‑1 (MQ=255)
                                  aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                       >  2:703123/1‑100 (MQ=255)
                                   cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                      <  2:799506/100‑1 (MQ=255)
                                   cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                      <  1:44540/100‑1 (MQ=255)
                                    ggTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                     <  2:1365386/100‑1 (MQ=255)
                                      tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                   >  1:2902318/1‑100 (MQ=255)
                                       aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGt                                                  >  2:105536/1‑100 (MQ=255)
                                        aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 <  2:2644765/100‑1 (MQ=255)
                                         tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                >  1:526072/1‑100 (MQ=255)
                                         tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                >  2:1439451/1‑100 (MQ=255)
                                         tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACAGCACATGAATCCAAAGCTGGATTAATCTTCTGTGATAg                                                   >  2:2499356/1‑97 (MQ=255)
                                          ttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGa                                               >  1:2129462/1‑100 (MQ=255)
                                           tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                              <  1:1704785/100‑1 (MQ=255)
                                             cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCgtt                                              >  1:1817755/1‑97 (MQ=255)
                                             cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  1:2495592/1‑100 (MQ=255)
                                             cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  2:1082071/1‑100 (MQ=255)
                                              tGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                           <  1:422313/100‑1 (MQ=255)
                                               gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          >  1:2236248/1‑100 (MQ=255)
                                               gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          <  2:1431511/100‑1 (MQ=255)
                                                 ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                        >  2:1296513/1‑100 (MQ=255)
                                                  tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       >  1:2985857/1‑100 (MQ=255)
                                                  tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       <  2:2408314/100‑1 (MQ=255)
                                                  tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       >  2:21775/1‑100 (MQ=255)
                                                   aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       <  1:209616/100‑1 (MQ=255)
                                                   aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  2:2910997/100‑1 (MQ=255)
                                                   aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      >  2:665764/1‑100 (MQ=255)
                                                   aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  1:539453/100‑1 (MQ=255)
                                                    ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:8841/1‑100 (MQ=255)
                                                    ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     <  1:2670971/100‑1 (MQ=255)
                                                      aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   <  2:259793/100‑1 (MQ=255)
                                                      aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   <  2:2158763/100‑1 (MQ=255)
                                                        cTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 <  2:1455766/100‑1 (MQ=255)
                                                          ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  1:2608936/1‑100 (MQ=255)
                                                          ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  2:1644569/1‑100 (MQ=255)
                                                             tCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                            <  1:547220/100‑1 (MQ=255)
                                                             tCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                            <  2:2608921/100‑1 (MQ=255)
                                                                 cTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCAAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCtt                         >  2:438827/1‑99 (MQ=255)
                                                                  tATTACGTCTACCAGATCACGCCATACG‑GCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                       <  2:1243984/100‑1 (MQ=255)
                                                                   aTTACGTCTACCAGATCACGCCATACG‑TCTACGGCACGTGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCTCCTTAcc                      >  2:1607678/1‑100 (MQ=255)
                                                                    ttACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCt                     >  1:1439780/1‑100 (MQ=255)
                                                                     tACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCt                     >  1:1679269/1‑100 (MQ=255)
                                                                     tACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATAAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    >  1:1573902/1‑100 (MQ=255)
                                                                      atctCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATCGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                   <  2:1974328/97‑1 (MQ=255)
                                                                       cGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  <  1:105538/100‑1 (MQ=255)
                                                                       cGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  <  1:1837716/100‑1 (MQ=255)
                                                                       cGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  >  2:2428536/1‑100 (MQ=255)
                                                                        gTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                 <  2:2441248/100‑1 (MQ=255)
                                                                        gTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                 >  1:1578067/1‑100 (MQ=255)
                                                                         tCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  1:2806780/100‑1 (MQ=255)
                                                                          cTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCaagat               >  1:1974340/1‑96 (MQ=255)
                                                                           tACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  1:2683802/1‑100 (MQ=255)
                                                                            aCCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc             >  1:2237947/1‑100 (MQ=255)
                                                                              cAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCtt           >  2:1863013/1‑100 (MQ=255)
                                                                              cAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCtt           >  2:608108/1‑100 (MQ=255)
                                                                               aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          <  2:1052242/100‑1 (MQ=255)
                                                                                gATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCg         <  1:1485394/100‑1 (MQ=255)
                                                                                 aTCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg        <  1:828289/100‑1 (MQ=255)
                                                                                  tCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       >  2:2251968/1‑100 (MQ=255)
                                                                                   cACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       >  2:2831227/1‑99 (MQ=255)
                                                                                    aCGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg     <  1:608111/100‑1 (MQ=255)
                                                                                     cGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       >  2:38036/1‑97 (MQ=255)
                                                                                     cGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGc    >  1:985385/1‑100 (MQ=255)
                                                                                      gCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc   >  1:391249/1‑100 (MQ=255)
                                                                                       ccATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  >  2:2867744/1‑100 (MQ=255)
                                                                                       ccATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  <  1:205123/100‑1 (MQ=255)
                                                                                       ccATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  >  2:2508749/1‑100 (MQ=255)
                                                                                       ccATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATCCAGCATCTTACCTCAGGCACCTTCGGGTGCCt  >  2:2868611/1‑100 (MQ=255)
                                                                                              |                                                                                             
CCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCT  >  NC_000913/4215084‑4215271

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑