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
                                                                                               |                                                                                             
tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTAc                                                                                          <  1:2443285/100‑1 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                         <  2:2520692/100‑1 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                         <  2:2505992/100‑1 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                         >  2:2566877/1‑100 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                         >  1:2223637/1‑100 (MQ=255)
 ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACg                                                                                         >  1:2328216/1‑100 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACgg                                                                                        >  2:2799067/1‑100 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACgg                                                                                        >  1:958297/1‑100 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGc                                                                                       <  1:2508682/100‑1 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGc                                                                                       >  1:871318/1‑100 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGc                                                                                       >  2:835992/1‑100 (MQ=255)
     aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcac                                                                                     <  1:1282504/100‑1 (MQ=255)
      aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcaca                                                                                    >  1:1150385/1‑100 (MQ=255)
        aTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATg                                                                                  >  1:1271452/1‑100 (MQ=255)
            acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATc                                                                              <  2:1246330/100‑1 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                             >  2:1474954/1‑100 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                             <  2:2328229/100‑1 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                             <  2:2237139/100‑1 (MQ=255)
                cgagcgagCTGGCGTAGTCACGGTAATTGACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAAc                                                                          <  2:710950/100‑1 (MQ=255)
                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGc                                                                        >  2:2696849/1‑100 (MQ=255)
                     gagtttgCGTAGTCACGGTAATTTACTGGTTACCAACTGGCGCAACTATTACATCTACCAGATCACGCCATACGTCTACGGCACATGAACCCAACGCTgg                                                                     <  1:1038239/94‑1 (MQ=255)
                     gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTgg                                                                     >  2:3011117/1‑100 (MQ=255)
                      agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGa                                                                    >  1:1676240/1‑100 (MQ=255)
                         tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTa                                                                 <  1:689046/100‑1 (MQ=255)
                            ngTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              >  2:2719715/2‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              >  1:1426172/1‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              <  1:80527/100‑1 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              >  2:1518170/1‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                              <  1:731301/100‑1 (MQ=255)
                             gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtct                                                             <  1:532/100‑1 (MQ=255)
                               aGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttc                                                           >  2:1801324/1‑100 (MQ=255)
                                gTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtcttct                                                          >  1:765693/1‑100 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        <  1:1435414/100‑1 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                        >  2:853038/1‑100 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGAACATGAATCCAACGCTGGATTACTCTTCtgt                                                        >  2:751314/1‑100 (MQ=255)
                                   aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                       <  2:393550/100‑1 (MQ=255)
                                    cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                                                      >  1:1192573/1‑100 (MQ=255)
                                    cGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCAATACGTCTACGGCACATGAATCCAACGCTGGATTCATCTTCTGTGa                                                      >  2:1906854/1‑100 (MQ=255)
                                       tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                   >  2:1739145/1‑100 (MQ=255)
                                        aaTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                   >  1:2203521/1‑100 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 <  2:1243446/100‑1 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 <  2:284829/100‑1 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                 >  1:267269/1‑100 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                >  1:2570490/1‑100 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAtgg                                            >  1:1853284/1‑98 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                            >  2:759034/1‑100 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          >  2:2985924/1‑100 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                          >  1:1757703/1‑100 (MQ=255)
                                                 tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCTAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                         <  2:2158065/100‑1 (MQ=255)
                                                 tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                         <  2:2158958/100‑1 (MQ=255)
                                                 tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                         <  1:1351522/100‑1 (MQ=255)
                                                  ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                        <  1:1058787/100‑1 (MQ=255)
                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       <  1:1154968/100‑1 (MQ=255)
                                                   cttccgatctcTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                       <  1:2951595/90‑1 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                      <  2:1599727/100‑1 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  2:1984184/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                     >  1:322417/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACACCATACGTATACGGCACACGAATCCAACGATGGATTAATCCTCTGTGATAGACGGTCGttctgcg                                     >  1:1937732/1‑95 (MQ=255)
                                                      cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                    >  1:977016/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   >  2:244667/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   >  1:3039541/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                   >  2:1076211/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                 >  1:1731514/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAGCGTTAAGCGATTc                                 <  2:2053069/100‑1 (MQ=255)
                                                          tGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                >  2:910766/1‑100 (MQ=255)
                                                          tGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa                                <  2:2411930/100‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               <  1:49911/100‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                               >  1:2424838/1‑100 (MQ=255)
                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc                              <  1:2031040/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             <  1:883442/100‑1 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                             >  2:2435415/1‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGAGAtcggaag                             >  2:2951614/1‑94 (MQ=255)
                                                               cAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                           >  2:2118811/1‑100 (MQ=255)
                                                                aaCTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCt                          >  2:1169972/1‑100 (MQ=255)
                                                                  cTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTa                        >  1:169142/1‑100 (MQ=255)
                                                                   tATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                       <  2:1150391/100‑1 (MQ=255)
                                                                   tATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                       >  1:1136955/1‑100 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    >  2:1887757/1‑100 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    <  2:339258/100‑1 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                    <  2:337196/100‑1 (MQ=255)
                                                                       aCGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                   <  2:751677/100‑1 (MQ=255)
                                                                       aCGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                   >  1:880450/1‑100 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  <  1:1724547/100‑1 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                  >  2:1916512/1‑100 (MQ=255)
                                                                         gTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                 <  2:948312/100‑1 (MQ=255)
                                                                          tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  1:853033/100‑1 (MQ=255)
                                                                          tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                >  2:1263189/1‑100 (MQ=255)
                                                                          tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  1:2512387/100‑1 (MQ=255)
                                                                          tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGc                <  2:1965326/100‑1 (MQ=255)
                                                                           cTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa               <  2:1026910/100‑1 (MQ=255)
                                                                           cTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa               >  1:2901362/1‑100 (MQ=255)
                                                                            tACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGATGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAc              >  1:1672507/1‑100 (MQ=255)
                                                                                aGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          >  2:1016423/1‑100 (MQ=255)
                                                                                aGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          >  2:560437/1‑100 (MQ=255)
                                                                                aGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTc          >  2:916720/1‑100 (MQ=255)
                                                                                  aTCACGCCATACGTCTACGGCACATGACTCCAACTCTGGATTAATCTTCTGTGCTAGTCGATCGTTAAGCGATTCAGCACCTTTCCTCAGGCACCTTCg         >  1:564098/1‑99 (MQ=255)
                                                                                   tCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg       >  1:834998/1‑100 (MQ=255)
                                                                                    cACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt      >  2:2274652/1‑100 (MQ=255)
                                                                                       gCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGcc   <  2:2044126/100‑1 (MQ=255)
                                                                                        ccATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCt  >  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: