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
                                                                                              |                                                                                             
ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                         >  2:2308637/1‑100 (MQ=255)
ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                         >  1:1446252/1‑100 (MQ=255)
  gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGc                                                                                       <  1:2951397/100‑1 (MQ=255)
   cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGca                                                                                      >  2:205121/1‑100 (MQ=255)
    aaaaTACACCGCGAGCTAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcac                                                                                     <  2:1743528/100‑1 (MQ=255)
         acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGa                                                                                 >  1:2045151/1‑99 (MQ=255)
           acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATc                                                                              <  1:449941/100‑1 (MQ=255)
            ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                             <  2:2904601/100‑1 (MQ=255)
             cgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCa                                                                            >  1:1375003/1‑100 (MQ=255)
               cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAAc                                                                          >  1:1916970/1‑100 (MQ=255)
                   cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTg                                                                      <  2:769325/100‑1 (MQ=255)
                    gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTgg                                                                     <  2:1274728/100‑1 (MQ=255)
                    gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTg                                                                      >  2:579445/1‑99 (MQ=255)
                     cgatctcGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGa                                                                    <  2:2095408/94‑1 (MQ=255)
                      gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAt                                                                   <  2:2194453/100‑1 (MQ=255)
                        tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTa                                                                 <  2:591370/100‑1 (MQ=255)
                         ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTaa                                                                >  1:791115/1‑100 (MQ=255)
                         ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTaa                                                                <  1:1501163/100‑1 (MQ=255)
                         ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTaa                                                                >  1:1660439/1‑100 (MQ=255)
                          gCGTAGTCACGGTAATTTACTGTTTACGAAATGGCTAAACTCTTAGGTCTACCAGATTACGCCATAAGTGTGCGGCACATGAATCTTACGCAgg                                                                     >  2:2904144/1‑94 (MQ=255)
                          gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAt                                                               >  1:1233212/1‑100 (MQ=255)
                           cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCAGACGTCTACGGCACATGAATCCAACGCTGGaagatcg                                                              >  1:2095420/1‑94 (MQ=255)
                             tAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtctt                                                            >  2:1151681/1‑100 (MQ=255)
                              aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                           >  2:2406768/1‑100 (MQ=255)
                               gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                          >  1:2004439/1‑100 (MQ=255)
                                tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  1:1675391/100‑1 (MQ=255)
                                tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                         <  2:791112/100‑1 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  2:2557452/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:259795/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:85888/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:1316820/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:1316448/1‑100 (MQ=255)
                                 cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  1:319917/1‑100 (MQ=255)
                                  aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                       <  1:2046059/100‑1 (MQ=255)
                                  aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                       >  2:703123/1‑100 (MQ=255)
                                   cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                      <  2:799506/100‑1 (MQ=255)
                                   cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                      <  1:44540/100‑1 (MQ=255)
                                    ggTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                     <  2:1365386/100‑1 (MQ=255)
                                      tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                   >  1:2902318/1‑100 (MQ=255)
                                       aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGt                                                  >  2:105536/1‑100 (MQ=255)
                                        aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 <  2:2644765/100‑1 (MQ=255)
                                         tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                >  1:526072/1‑100 (MQ=255)
                                         tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                >  2:1439451/1‑100 (MQ=255)
                                         tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACAGCACATGAATCCAAAGCTGGATTAATCTTCTGTGATAg                                                   >  2:2499356/1‑97 (MQ=255)
                                          ttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGa                                               >  1:2129462/1‑100 (MQ=255)
                                           tACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                                              <  1:1704785/100‑1 (MQ=255)
                                             cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCgtt                                              >  1:1817755/1‑97 (MQ=255)
                                             cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  1:2495592/1‑100 (MQ=255)
                                             cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  2:1082071/1‑100 (MQ=255)
                                              tGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                           <  1:422313/100‑1 (MQ=255)
                                               gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          >  1:2236248/1‑100 (MQ=255)
                                               gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          <  2:1431511/100‑1 (MQ=255)
                                                 ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                        >  2:1296513/1‑100 (MQ=255)
                                                  tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       >  1:2985857/1‑100 (MQ=255)
                                                  tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       <  2:2408314/100‑1 (MQ=255)
                                                  tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       >  2:21775/1‑100 (MQ=255)
                                                   aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       <  1:209616/100‑1 (MQ=255)
                                                   aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  2:2910997/100‑1 (MQ=255)
                                                   aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      >  2:665764/1‑100 (MQ=255)
                                                   aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  1:539453/100‑1 (MQ=255)
                                                    ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:8841/1‑100 (MQ=255)
                                                    ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     <  1:2670971/100‑1 (MQ=255)
                                                      aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   <  2:259793/100‑1 (MQ=255)
                                                      aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   <  2:2158763/100‑1 (MQ=255)
                                                        cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 <  2:1455766/100‑1 (MQ=255)
                                                          ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  1:2608936/1‑100 (MQ=255)
                                                          ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  2:1644569/1‑100 (MQ=255)
                                                             tCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                            <  1:547220/100‑1 (MQ=255)
                                                             tCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                            <  2:2608921/100‑1 (MQ=255)
                                                                 cTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCAAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCtt                         >  2:438827/1‑99 (MQ=255)
                                                                  tATTACGTCTACCAGATCACGCCATACGGCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                       <  2:1243984/100‑1 (MQ=255)
                                                                   aTTACGTCTACCAGATCACGCCATACGTCTACGGCACGTGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCTCCTTAcc                      >  2:1607678/1‑100 (MQ=255)
                                                                    ttACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCt                     >  1:1439780/1‑100 (MQ=255)
                                                                     tACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCt                     >  1:1679269/1‑100 (MQ=255)
                                                                     tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATAAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    >  1:1573902/1‑100 (MQ=255)
                                                                      atctCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATCGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                   <  2:1974328/97‑1 (MQ=255)
                                                                       cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  <  1:105538/100‑1 (MQ=255)
                                                                       cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  <  1:1837716/100‑1 (MQ=255)
                                                                       cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  >  2:2428536/1‑100 (MQ=255)
                                                                        gTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                 <  2:2441248/100‑1 (MQ=255)
                                                                        gTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                 >  1:1578067/1‑100 (MQ=255)
                                                                         tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  1:2806780/100‑1 (MQ=255)
                                                                          cTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCaagat               >  1:1974340/1‑96 (MQ=255)
                                                                           tACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  1:2683802/1‑100 (MQ=255)
                                                                            aCCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc             >  1:2237947/1‑100 (MQ=255)
                                                                              cAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCtt           >  2:1863013/1‑100 (MQ=255)
                                                                              cAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCtt           >  2:608108/1‑100 (MQ=255)
                                                                               aGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          <  2:1052242/100‑1 (MQ=255)
                                                                                gATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCg         <  1:1485394/100‑1 (MQ=255)
                                                                                 aTCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCgg        <  1:828289/100‑1 (MQ=255)
                                                                                  tCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       >  2:2251968/1‑100 (MQ=255)
                                                                                   cACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       >  2:2831227/1‑99 (MQ=255)
                                                                                    aCGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg     <  1:608111/100‑1 (MQ=255)
                                                                                     cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       >  2:38036/1‑97 (MQ=255)
                                                                                     cGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGc    >  1:985385/1‑100 (MQ=255)
                                                                                      gCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc   >  1:391249/1‑100 (MQ=255)
                                                                                       ccATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  >  2:2867744/1‑100 (MQ=255)
                                                                                       ccATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  <  1:205123/100‑1 (MQ=255)
                                                                                       ccATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  >  2:2508749/1‑100 (MQ=255)
                                                                                       ccATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATCCAGCATCTTACCTCAGGCACCTTCGGGTGCCt  >  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: