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.99.3% 634.5 / ‑2.6 137coding (899/930 nt)metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (78/58);  total (79/58)
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.02e-01

Only 100 of 137 total aligned reads displayed.

TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCT  >  NC_000913/4215083‑4215271
                                                                                               |                                                                                             
tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTAc                                                                                          <  1:2443285/100‑1 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                         <  2:2520692/100‑1 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                         <  2:2505992/100‑1 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                         >  2:2566877/1‑100 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                         >  1:2223637/1‑100 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACg                                                                                         >  1:2328216/1‑100 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACgg                                                                                        >  2:2799067/1‑100 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACgg                                                                                        >  1:958297/1‑100 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                       <  1:2508682/100‑1 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                       >  1:871318/1‑100 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGc                                                                                       >  2:835992/1‑100 (MQ=255)
     aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcac                                                                                     <  1:1282504/100‑1 (MQ=255)
      aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGcaca                                                                                    >  1:1150385/1‑100 (MQ=255)
        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATg                                                                                  >  1:1271452/1‑100 (MQ=255)
            acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATc                                                                              <  2:1246330/100‑1 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                             >  2:1474954/1‑100 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                             <  2:2328229/100‑1 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATcc                                                                             <  2:2237139/100‑1 (MQ=255)
                cgagcgagCTGGCGTAGTCACGGTAATTGACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAAc                                                                          <  2:710950/100‑1 (MQ=255)
                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGc                                                                        >  2:2696849/1‑100 (MQ=255)
                     gagtttgCGTAGTCACGGTAATTTACTGGTTACCAACTGGCGCAACTATTACATCTACCAGATCACGCCATACG‑TCTACGGCACATGAACCCAACGCTgg                                                                     <  1:1038239/94‑1 (MQ=255)
                     gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTgg                                                                     >  2:3011117/1‑100 (MQ=255)
                      agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGa                                                                    >  1:1676240/1‑100 (MQ=255)
                         tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTa                                                                 <  1:689046/100‑1 (MQ=255)
                            ngTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              >  2:2719715/2‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              >  1:1426172/1‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              <  1:80527/100‑1 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              >  2:1518170/1‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              <  1:731301/100‑1 (MQ=255)
                             gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtct                                                             <  1:532/100‑1 (MQ=255)
                               aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                           >  2:1801324/1‑100 (MQ=255)
                                gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                          >  1:765693/1‑100 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        <  1:1435414/100‑1 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  2:853038/1‑100 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGAACATGAATCCAACGCTGGATTACTCTTCtgt                                                        >  2:751314/1‑100 (MQ=255)
                                   aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                       <  2:393550/100‑1 (MQ=255)
                                    cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                      >  1:1192573/1‑100 (MQ=255)
                                    cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCAATACG‑TCTACGGCACATGAATCCAACGCTGGATTCATCTTCTGTGa                                                      >  2:1906854/1‑100 (MQ=255)
                                       tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                   >  2:1739145/1‑100 (MQ=255)
                                        aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                   >  1:2203521/1‑100 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 <  2:1243446/100‑1 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 <  2:284829/100‑1 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 >  1:267269/1‑100 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                >  1:2570490/1‑100 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAtgg                                            >  1:1853284/1‑98 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  2:759034/1‑100 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          >  2:2985924/1‑100 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          >  1:1757703/1‑100 (MQ=255)
                                                 tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCTAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                         <  2:2158065/100‑1 (MQ=255)
                                                 tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                         <  2:2158958/100‑1 (MQ=255)
                                                 tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                         <  1:1351522/100‑1 (MQ=255)
                                                  ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                        <  1:1058787/100‑1 (MQ=255)
                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       <  1:1154968/100‑1 (MQ=255)
                                                   cttccgatctcTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       <  1:2951595/90‑1 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  2:1599727/100‑1 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  2:1984184/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:322417/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACACCATACG‑TATACGGCACACGAATCCAACGATGGATTAATCCTCTGTGATAGACGGTCGttctgcg                                     >  1:1937732/1‑95 (MQ=255)
                                                      cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                    >  1:977016/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   >  2:244667/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   >  1:3039541/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   >  2:1076211/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 >  1:1731514/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAGCGTTAAGCGATTc                                 <  2:2053069/100‑1 (MQ=255)
                                                          tGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                >  2:910766/1‑100 (MQ=255)
                                                          tGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                <  2:2411930/100‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               <  1:49911/100‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  1:2424838/1‑100 (MQ=255)
                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                              <  1:2031040/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             <  1:883442/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  2:2435415/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGAGAtcggaag                             >  2:2951614/1‑94 (MQ=255)
                                                               cAACTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                           >  2:2118811/1‑100 (MQ=255)
                                                                aaCTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCt                          >  2:1169972/1‑100 (MQ=255)
                                                                  cTATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTa                        >  1:169142/1‑100 (MQ=255)
                                                                   tATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                       <  2:1150391/100‑1 (MQ=255)
                                                                   tATTACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                       >  1:1136955/1‑100 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    >  2:1887757/1‑100 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    <  2:339258/100‑1 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    <  2:337196/100‑1 (MQ=255)
                                                                       aCGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                   <  2:751677/100‑1 (MQ=255)
                                                                       aCGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                   >  1:880450/1‑100 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  <  1:1724547/100‑1 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  >  2:1916512/1‑100 (MQ=255)
                                                                         gTCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                 <  2:948312/100‑1 (MQ=255)
                                                                          tCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  1:853033/100‑1 (MQ=255)
                                                                          tCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                >  2:1263189/1‑100 (MQ=255)
                                                                          tCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  1:2512387/100‑1 (MQ=255)
                                                                          tCTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  2:1965326/100‑1 (MQ=255)
                                                                           cTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa               <  2:1026910/100‑1 (MQ=255)
                                                                           cTACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa               >  1:2901362/1‑100 (MQ=255)
                                                                            tACCAGATCACGCCATACG‑TCTACGGCACATGAATCCAACGATGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  1:1672507/1‑100 (MQ=255)
                                                                                aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          >  2:1016423/1‑100 (MQ=255)
                                                                                aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          >  2:560437/1‑100 (MQ=255)
                                                                                aGATCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          >  2:916720/1‑100 (MQ=255)
                                                                                  aTCACGCCATACG‑TCTACGGCACATGACTCCAACTCTGGATTAATCTTCTGTGCTAGTCGATCGTTAAGCGATTCAGCACCTTTCCTCAGGCACCTTCg         >  1:564098/1‑99 (MQ=255)
                                                                                   tCACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       >  1:834998/1‑100 (MQ=255)
                                                                                    cACGCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt      >  2:2274652/1‑100 (MQ=255)
                                                                                       gCCATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc   <  2:2044126/100‑1 (MQ=255)
                                                                                        ccATACG‑TCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  >  2:2680015/1‑100 (MQ=255)
                                                                                               |                                                                                             
TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCT  >  NC_000913/4215083‑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: ‑