Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,215,147 A→C 100% N290H (AAC→CAC)  metA → homoserine O‑transsuccinylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,215,1470AC97.6% 581.4 / 3.1 167N290H (AAC→CAC) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (3/1);  new base C (90/73);  total (93/74)
Fisher's exact test for biased strand distribution p-value = 6.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.98e-01

Only 100 of 167 total aligned reads displayed.

GGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGC  >  NC_000913/4215050‑4215243
                                                                                                 |                                                                                                
ggATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCAc                                                                                                <  2:1808313/100‑1 (MQ=255)
 gATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACt                                                                                               <  2:2807671/100‑1 (MQ=255)
 gATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACt                                                                                               <  1:420776/100‑1 (MQ=255)
   tGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTAt                                                                                             >  1:3638080/1‑100 (MQ=255)
   tGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTAt                                                                                             >  2:516736/1‑100 (MQ=255)
   tGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACAAACTGGCTCCACTAt                                                                                             >  2:2067286/1‑100 (MQ=255)
     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                           <  2:3414648/100‑1 (MQ=255)
     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                           <  1:2379336/100‑1 (MQ=255)
      aCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTAc                                                                                          <  1:3682758/100‑1 (MQ=255)
      aCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTAc                                                                                          >  1:3652019/1‑100 (MQ=255)
       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACg                                                                                         <  1:3490623/100‑1 (MQ=255)
       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACg                                                                                         <  1:1926563/100‑1 (MQ=255)
        cGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGt                                                                                        <  2:2834768/100‑1 (MQ=255)
        cGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGt                                                                                        >  2:1782946/1‑100 (MQ=255)
          tataACTATTTCCCGCACTATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCt                                                                                      <  2:2750030/100‑1 (MQ=255)
            taACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTCTTACGTCTAc                                                                                    >  2:592190/1‑100 (MQ=255)
             aaCTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAcc                                                                                   <  1:1315203/100‑1 (MQ=255)
             aaCTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAcc                                                                                   <  1:335058/100‑1 (MQ=255)
             aaCTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGCCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAcc                                                                                   <  2:2117898/100‑1 (MQ=255)
                tATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAg                                                                                 >  1:2745763/1‑99 (MQ=255)
                tATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGa                                                                                >  2:2941504/1‑100 (MQ=255)
                 aTTTCCCGCACAATGATCCGCAAAATACAACGAGAGCCAGCTCACGTAGTCACGGTAATTTACTGTTTACCAACTGTCTCCACCATTAAGTCTACCa                                                                                  >  2:2233984/1‑97 (MQ=255)
                     cccGCACAATGATCCGCGAAATACTCCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATc                                                                              <  2:1296777/97‑1 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                           >  1:1801789/1‑100 (MQ=255)
                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          <  2:3792655/100‑1 (MQ=255)
                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          <  1:1011509/100‑1 (MQ=255)
                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          <  2:3222976/100‑1 (MQ=255)
                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          >  1:2494684/1‑100 (MQ=255)
                       cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCa                                                                             >  2:1589114/1‑96 (MQ=255)
                       cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGcc                                                                         <  1:3011632/100‑1 (MQ=255)
                       cGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGcc                                                                         <  1:451819/100‑1 (MQ=255)
                        gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCa                                                                        >  2:416239/1‑100 (MQ=255)
                           caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATAc                                                                     >  2:1195827/1‑100 (MQ=255)
                           caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATAc                                                                     >  1:3543345/1‑100 (MQ=255)
                           caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATAc                                                                     <  1:3102821/100‑1 (MQ=255)
                           caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCGCCACTATTACGTCTACCAGATCACGCCATAc                                                                     >  1:2506624/1‑100 (MQ=255)
                             aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGa                                                                   >  1:3602054/1‑100 (MQ=255)
                              tGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGAt                                                                  <  2:748188/100‑1 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 >  2:2907465/1‑100 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 >  2:2903961/1‑100 (MQ=255)
                                 tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTa                                                               >  1:486322/1‑100 (MQ=255)
                                  ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCtcc                                                              <  1:3541782/100‑3 (MQ=255)
                                  ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                              >  1:881246/1‑100 (MQ=255)
                                   cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGCTCTAc                                                              >  1:2346136/1‑99 (MQ=255)
                                   cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                             >  1:2109245/1‑100 (MQ=255)
                                   cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                             >  2:624667/1‑100 (MQ=255)
                                     cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGc                                                           <  1:524421/100‑1 (MQ=255)
                                     cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGc                                                           <  1:2326614/100‑1 (MQ=255)
                                     cAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGc                                                           <  1:2479990/100‑1 (MQ=255)
                                       aaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGcac                                                         >  2:756024/1‑100 (MQ=255)
                                          tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                                                      <  1:1604318/100‑1 (MQ=255)
                                          tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                                                      >  1:2192956/1‑100 (MQ=255)
                                            cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTTCGTCTACCAGATAACGCCATACGATCTACGTCACATGa                                                     >  2:2670394/1‑99 (MQ=255)
                                              ncgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                  >  1:107175/2‑100 (MQ=255)
                                              ccGCGAGCGAGCTAGCGTAGTCACGGTAATTTACTGTTGACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                  >  2:2357635/1‑100 (MQ=255)
                                               cggtaccgagcTGGCGTAGTCACCGTAATGTACTGTTTACCAACTGGCTCCACTATTACCTCTACCAGATCACCCCATACGATCTACGGCACATGAATcc                                                 <  2:3220421/94‑1 (MQ=255)
                                               cgcgAGCGAGCTGGCGTATTCACAGTAATTTACTGTTTACCAACTGGCTCCACTATTACTTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                  <  2:1481216/99‑1 (MQ=255)
                                                     cgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGc                                            >  1:782409/1‑99 (MQ=255)
                                                      gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg                                          <  2:1603254/100‑1 (MQ=255)
                                                      gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg                                          <  2:2144052/100‑1 (MQ=255)
                                                      gagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTg                                          <  1:1779948/100‑1 (MQ=255)
                                                       agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCt                                           >  1:1339238/1‑98 (MQ=255)
                                                       agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTgg                                         <  2:232864/100‑1 (MQ=255)
                                                       agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTgg                                         <  2:91546/100‑1 (MQ=255)
                                                       agcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTAAGTCTACCAGATCACGCCATACGCTCTACGGCACATGATTCCAACGCTg                                          >  2:3597920/1‑99 (MQ=255)
                                                        gcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGa                                        >  2:3501045/1‑100 (MQ=255)
                                                         cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGa                                        >  2:3085935/1‑99 (MQ=255)
                                                         cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt                                       <  1:2971614/100‑1 (MQ=255)
                                                         cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATAACGCCATACGCTCTACGGCACATGAATCCAACGCTgatt                                       >  1:1267987/1‑97 (MQ=255)
                                                          tGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                      <  2:1714401/100‑1 (MQ=255)
                                                           ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTa                                     >  2:2987986/1‑100 (MQ=255)
                                                           ggCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAtt                                      >  1:873031/1‑99 (MQ=255)
                                                                aGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtctt                                >  2:2024661/1‑100 (MQ=255)
                                                                 gTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                               >  1:3021922/1‑100 (MQ=255)
                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                             <  2:3453899/100‑1 (MQ=255)
                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                             >  1:3786054/1‑100 (MQ=255)
                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                             >  2:1140071/1‑100 (MQ=255)
                                                                     cGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                           >  2:3722264/1‑100 (MQ=255)
                                                                     cGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgt                            >  1:2082347/1‑99 (MQ=255)
                                                                      ggTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGa                          <  1:601475/100‑1 (MQ=255)
                                                                       gTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGAt                         <  2:832485/100‑1 (MQ=255)
                                                                         aaTTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACACCATACGATCCACGGCACATGAATCCAACGCTGGATTAATTTTCTGCGATa                        >  2:1890112/1‑99 (MQ=255)
                                                                           tttACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCATCTGTGATAGTc                     <  2:1324087/100‑1 (MQ=255)
                                                                             tACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTTAtttt                      >  1:2455243/1‑94 (MQ=255)
                                                                             tACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCCGGATTAATCTTCTGTGATAGTc                     >  1:2433060/1‑98 (MQ=255)
                                                                               cTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGAt                  >  1:3056905/1‑99 (MQ=255)
                                                                               cTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATc                 <  1:3186633/100‑1 (MQ=255)
                                                                                tGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                >  1:546836/1‑100 (MQ=255)
                                                                                   ttACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa             <  2:3429565/100‑1 (MQ=255)
                                                                                   ttACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa             <  2:3802345/100‑1 (MQ=255)
                                                                                   ttACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa             <  1:2952589/100‑1 (MQ=255)
                                                                                    tACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGtt              >  1:3731116/1‑98 (MQ=255)
                                                                                     aCCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg           >  2:2727000/1‑100 (MQ=255)
                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc          >  1:3133742/1‑100 (MQ=255)
                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc          >  1:3422413/1‑100 (MQ=255)
                                                                                      ccAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGc          >  2:1965335/1‑100 (MQ=255)
                                                                                       cAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCg         >  1:974780/1‑100 (MQ=255)
                                                                                            ggCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa    >  1:2938110/1‑100 (MQ=255)
                                                                                             gCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg   >  2:1159668/1‑100 (MQ=255)
                                                                                              cTCCACTATTGCGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAACCTTCTGTGATAGTCGATCGTTAAGCGATTCAGc  <  2:2506608/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGC  >  NC_000913/4215050‑4215243

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