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.94.2% 556.4 / 25.4 137coding (899/930 nt)metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (4/4);  new base . (75/54);  total (79/58)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 137 total aligned reads displayed.

TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTT  >  NC_000913/4215083‑4215272
                                                                                               |                                                                                              
tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATc                                                                                              <  1:792218/98‑1 (MQ=255)
tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTAc                                                                                           >  1:1302415/1‑100 (MQ=255)
tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTAc                                                                                           <  1:739318/100‑1 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACgg                                                                                         >  1:1039155/1‑100 (MQ=255)
  cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACgg                                                                                         >  1:2351817/1‑100 (MQ=255)
   gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGc                                                                                        <  1:399891/100‑1 (MQ=255)
    cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGca                                                                                       >  1:2680775/1‑100 (MQ=255)
    cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGca                                                                                       >  1:1907516/1‑100 (MQ=255)
      aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGcaca                                                                                     >  2:1911039/1‑100 (MQ=255)
       aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACAt                                                                                    <  1:2380853/100‑1 (MQ=255)
        aTACCCCGCGGGCGAGCTGGCGTAGTCACGGTAATTTACCGTTTACCGACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATg                                                                                   <  1:1104558/100‑1 (MQ=255)
            acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATc                                                                               >  2:1317705/1‑100 (MQ=255)
            acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATc                                                                               <  1:2881635/100‑1 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                              >  2:1897702/1‑100 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                              <  1:512068/100‑1 (MQ=255)
             ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATcc                                                                              >  2:2813992/1‑100 (MQ=255)
                cgagcgagCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAAc                                                                           >  2:1774320/1‑100 (MQ=255)
                  agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGc                                                                         >  2:1457592/1‑100 (MQ=255)
                    cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTg                                                                       >  1:1411381/1‑100 (MQ=255)
                     gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTAGCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTg                                                                       <  1:1581767/99‑1 (MQ=255)
                      agcTGGCGTAGTCACGGTCATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTgca                                                                     >  1:2497418/1‑98 (MQ=255)
                       gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAt                                                                    <  2:742163/100‑1 (MQ=255)
                        cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAtt                                                                   <  1:197299/100‑1 (MQ=255)
                         tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTa                                                                  >  2:1306151/1‑100 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTaa                                                                 >  1:194333/1‑100 (MQ=255)
                          ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTaa                                                                 >  2:2781831/1‑100 (MQ=255)
                            cggaGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               <  1:295174/97‑1 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               >  2:534407/1‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               >  2:759569/1‑100 (MQ=255)
                            cGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAAtc                                                               >  2:708263/1‑100 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                          >  1:1102062/1‑100 (MQ=255)
                                 tCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                                                          <  2:1204377/100‑1 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                         <  1:2376548/100‑1 (MQ=255)
                                  cACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                                                         <  1:1405495/100‑1 (MQ=255)
                                   aCGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                                                        <  2:1022962/100‑1 (MQ=255)
                                     tcttccgatcTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                                                      <  1:1121347/91‑1 (MQ=255)
                                       tAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAg                                                    <  2:2785614/100‑1 (MQ=255)
                                         aTTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                  <  2:194335/100‑1 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtagatcgga                                                 >  2:512379/1‑92 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTc                                                  <  2:1654692/99‑1 (MQ=255)
                                          tttACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                                                 >  2:1588170/1‑100 (MQ=255)
                                             aCTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATc                                              <  2:1356501/100‑1 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                             >  2:2584147/1‑100 (MQ=255)
                                              cTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                                             <  1:958837/100‑1 (MQ=255)
                                               tGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGt                                            <  2:2147532/100‑1 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAgatcggaaga                                           >  2:1121349/1‑91 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt                                           >  1:2307450/1‑100 (MQ=255)
                                                gTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACCGCACATGAATCCAACGCCGCATTAATCTTCCGTGATAGTCGacc                                              >  1:2589185/1‑95 (MQ=255)
                                                 tttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa                                          <  2:2181865/100‑1 (MQ=255)
                                                  ttACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTaa                                         >  2:1443314/1‑100 (MQ=255)
                                                   tACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg                                        <  2:82453/100‑1 (MQ=255)
                                                    aCCAACTGGCTCAACTATTACGTCTACCAGGTCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGGGATAGTCGATCGTTAAGc                                       >  2:2613212/1‑100 (MQ=255)
                                                     ccAACTGGCTCACCTACTACGTCTACCAGATCACTCCATCCATCTGCGGCACATCAATCCAATGCTGCAGTAATCGTCTGCGATAAACGAGCGCTAAGc                                       >  1:1904016/1‑99 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc                                       >  1:460046/1‑100 (MQ=255)
                                                     ccAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg                                      >  2:2127852/1‑100 (MQ=255)
                                                      cAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa                                     >  1:706397/1‑100 (MQ=255)
                                                       aaCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt                                    >  1:427921/1‑100 (MQ=255)
                                                        aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAtt                                   >  1:1576478/1‑100 (MQ=255)
                                                        aCTGGCTCAACTATTACGTCTACCAGATCACGCCATACG‑‑‑‑ACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                                <  2:2320333/100‑1 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                  >  2:1174215/1‑100 (MQ=255)
                                                         cTGGCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc                                  <  2:685374/100‑1 (MQ=255)
                                                           ggCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg                                <  2:1481076/100‑1 (MQ=255)
                                                            gCTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTAGATCGTTAAGAGATTCAGc                               >  2:1694788/1‑100 (MQ=255)
                                                             ntCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGACCGTTAAGCGATTCAGCa                              >  1:808199/2‑100 (MQ=255)
                                                             cTCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCa                              >  1:1711622/1‑100 (MQ=255)
                                                              tCAATTATTACGTCAAACAGCTCCAGCAATACATCTAGGGCCTATGAAACCACCGCTGCATTAACATCATGTGATAGTAGATCCTTAAGCGATTCAGCa                              <  1:1694783/99‑1 (MQ=255)
                                                              tCAACTATTACGTCTACCATATCACGCCACACGTCTACGTCACATTAATCCAACGCTCCATTACTCTTCTGTGATAGTCGCTCGTTACGCGATTCAGCAc                             <  1:1008112/100‑1 (MQ=255)
                                                              tCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATGCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                             <  2:2497429/100‑1 (MQ=255)
                                                              tCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                             <  2:380300/100‑1 (MQ=255)
                                                              tCAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAc                             <  1:96716/100‑1 (MQ=255)
                                                               cAACTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCAcc                            >  2:335282/1‑100 (MQ=255)
                                                                  cTATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTa                         <  2:1802793/100‑1 (MQ=255)
                                                                   tATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                        >  1:2222230/1‑100 (MQ=255)
                                                                   tATTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                        >  1:2298309/1‑100 (MQ=255)
                                                                    aTTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                        >  1:1490791/1‑100 (MQ=255)
                                                                    aTTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAcc                       >  1:321813/1‑100 (MQ=255)
                                                                    aTTACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGAATAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTAc                        >  1:248149/1‑99 (MQ=255)
                                                                     ttACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCt                      <  2:662224/100‑1 (MQ=255)
                                                                     ttACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCt                      <  2:657048/100‑1 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                     >  2:1707147/1‑100 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                     >  1:149232/1‑100 (MQ=255)
                                                                      tACGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTc                     >  1:525032/1‑100 (MQ=255)
                                                                       aCGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCa                    <  2:337071/100‑1 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                   >  2:1559152/1‑100 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                   <  1:1739252/100‑1 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                   <  1:538577/100‑1 (MQ=255)
                                                                        cGTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAg                   <  1:538484/100‑1 (MQ=255)
                                                                         gTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                  >  2:866724/1‑100 (MQ=255)
                                                                         gTCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAgg                  >  2:2876899/1‑100 (MQ=255)
                                                                          tCTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGGTCGTTAAGCGCTTCAGCACCTTACCTCa                    >  1:895718/1‑97 (MQ=255)
                                                                           cTACCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCa                <  2:1411383/100‑1 (MQ=255)
                                                                             aCCAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCAcc              >  2:205182/1‑100 (MQ=255)
                                                                              ccAGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCt             >  1:1983172/1‑100 (MQ=255)
                                                                                aGATCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCCTCTGTGATAGTCGATCGTTAAGCGCTTCAGCACCTTACATCAGGCACCTTc           >  1:107479/1‑100 (MQ=255)
                                                                                   tCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg        <  2:2843148/100‑1 (MQ=255)
                                                                                   tCACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCggg        <  1:2881839/100‑1 (MQ=255)
                                                                                    cACGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGt       >  2:1268948/1‑100 (MQ=255)
                                                                                     aCGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg      <  1:2540145/100‑1 (MQ=255)
                                                                                     aCGCCATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTg      <  2:1821860/100‑1 (MQ=255)
                                                                                         cATACGTCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCtt  >  1:1596612/1‑100 (MQ=255)
                                                                                               |                                                                                              
TCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGCACCTTACCTCAGGCACCTTCGGGTGCCTT  >  NC_000913/4215083‑4215272

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