Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,215,147 A→C 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,1470AC99.2% 445.4 / ‑4.4 126N290H (AAC→CAC) metAhomoserine O‑transsuccinylase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (64/61);  total (64/62)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 126 total aligned reads displayed.

GGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGC  >  NC_000913/4215050‑4215243
                                                                                                 |                                                                                                
ggATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCAc                                                                                                <  2:2269475/100‑1 (MQ=255)
    gTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTAtt                                                                                            <  2:1098783/100‑1 (MQ=255)
     tACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTa                                                                                           <  1:132830/100‑1 (MQ=255)
      aCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTAc                                                                                          >  1:189407/1‑100 (MQ=255)
      aCCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTAc                                                                                          <  1:2762826/100‑1 (MQ=255)
       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTAc                                                                                          >  1:1751407/1‑99 (MQ=255)
       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACg                                                                                         <  1:550861/100‑1 (MQ=255)
       ccGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTAGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACg                                                                                         <  2:664705/100‑1 (MQ=255)
         gTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTc                                                                                       <  1:939983/100‑1 (MQ=255)
         gTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTc                                                                                       <  1:2600572/100‑1 (MQ=255)
         gTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTc                                                                                       <  1:2785760/100‑1 (MQ=255)
         gTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTc                                                                                       >  2:296790/1‑100 (MQ=255)
         gTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTACCAACTGGCTCCACTATTACGTCt                                                                                      <  2:1120734/100‑1 (MQ=255)
           ataACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTa                                                                                     >  1:2142822/1‑100 (MQ=255)
            taACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTa                                                                                     >  1:1811610/1‑99 (MQ=255)
            taACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTa                                                                                     >  1:3280142/1‑99 (MQ=255)
            taACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAc                                                                                    >  1:1742418/1‑100 (MQ=255)
            taACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTAc                                                                                    >  1:2099721/1‑100 (MQ=255)
               cTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAg                                                                                 <  1:3074307/100‑1 (MQ=255)
               cTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAg                                                                                 <  2:3194585/100‑1 (MQ=255)
                tATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGa                                                                                <  1:3209359/100‑1 (MQ=255)
                 aTTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCTGAt                                                                               >  1:2070746/1‑100 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCAc                                                                            >  2:2035136/1‑99 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                           >  2:1898618/1‑100 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                           >  2:884357/1‑100 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACg                                                                           >  2:71511/1‑100 (MQ=255)
                     cccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACAGGCTCCACTATTACGTCTACCAGATCACg                                                                           >  2:2025378/1‑100 (MQ=255)
                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          >  2:498760/1‑100 (MQ=255)
                      ccGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGc                                                                          >  2:84744/1‑100 (MQ=255)
                        gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCa                                                                        >  1:2176988/1‑100 (MQ=255)
                        gCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCa                                                                        <  2:2231056/100‑1 (MQ=255)
                         cacaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCAt                                                                       <  1:1640111/100‑1 (MQ=255)
                          acaATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATa                                                                      <  2:2696186/100‑1 (MQ=255)
                           caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCTACTGGCTCCACTATTACGTCTACCGGATCACGCCATAc                                                                     <  2:2071925/100‑1 (MQ=255)
                           caATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATAc                                                                     <  2:828695/100‑1 (MQ=255)
                             aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGa                                                                   >  1:2006756/1‑100 (MQ=255)
                             aTGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGa                                                                   <  1:1366061/100‑1 (MQ=255)
                              tGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGAt                                                                  >  1:1268381/1‑100 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 >  2:1667737/1‑100 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 <  1:1444854/100‑1 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 <  1:1880639/100‑1 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 >  2:3018165/1‑100 (MQ=255)
                               gATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATc                                                                 <  2:1178655/100‑1 (MQ=255)
                                aTCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCt                                                                >  1:2354606/1‑100 (MQ=255)
                                 tCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTa                                                               >  2:1608301/1‑100 (MQ=255)
                                  ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                              >  1:1552030/1‑100 (MQ=255)
                                  ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                              >  2:1357815/1‑100 (MQ=255)
                                  ccGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTAc                                                              <  2:2588362/100‑1 (MQ=255)
                                   cGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACg                                                             <  1:884359/100‑1 (MQ=255)
                                    gCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACgg                                                            >  1:921952/1‑100 (MQ=255)
                                      aaaaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGca                                                          <  1:249618/100‑1 (MQ=255)
                                        aaTACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGcac                                                         <  1:2124341/99‑1 (MQ=255)
                                          tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                                                      <  1:1792219/100‑1 (MQ=255)
                                          tACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                                                      <  1:593617/100‑1 (MQ=255)
                                           acacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATg                                                      <  1:2948833/99‑1 (MQ=255)
                                            cacCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGaa                                                    >  1:666378/1‑100 (MQ=255)
                                             agcgcgAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                   <  1:1735788/98‑1 (MQ=255)
                                             acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                   <  1:645897/100‑1 (MQ=255)
                                             acCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                   >  1:2900062/1‑100 (MQ=255)
                                              ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAAt                                                   >  1:3329460/1‑99 (MQ=255)
                                              ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                  >  1:37585/1‑100 (MQ=255)
                                              ccGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATc                                                  >  2:341684/1‑100 (MQ=255)
                                                tcgagcgagCTGGCGTAGTCACGGTAGTTTACTGTTTACCAACTGGCTGCACTATTACGTCTACCAGATCACGCCATACGATCTACGTTACATGAATcc                                                 >  2:3221665/2‑99 (MQ=255)
                                                gcgagcgaGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATccn                                                <  1:1705163/100‑2 (MQ=255)
                                                   agcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACg                                             >  2:260679/1‑100 (MQ=255)
                                                    gcgagcTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGc                                            <  1:1886221/100‑1 (MQ=255)
                                                         cTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGAt                                       <  2:3333486/100‑1 (MQ=255)
                                                            gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTaa                                    <  2:569271/100‑1 (MQ=255)
                                                            gCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGAGCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTaa                                    <  1:463790/100‑1 (MQ=255)
                                                              gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtc                                  >  1:2106170/1‑100 (MQ=255)
                                                              gTAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtc                                  <  1:226501/100‑1 (MQ=255)
                                                               tAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtct                                 <  1:2522409/100‑1 (MQ=255)
                                                               tAGTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtct                                 <  1:2070142/100‑1 (MQ=255)
                                                                 gTCACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                               >  2:61907/1‑100 (MQ=255)
                                                                 ccgatcttAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttc                               <  1:1022933/93‑1 (MQ=255)
                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAAtcttct                              >  1:1635813/1‑99 (MQ=255)
                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                             <  1:1544455/100‑1 (MQ=255)
                                                                   cACGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtg                             <  1:1608305/100‑1 (MQ=255)
                                                                     cGGTAATTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCtgtg                           <  1:495222/100‑1 (MQ=255)
                                                                          aTTTACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGt                      >  2:3109978/1‑100 (MQ=255)
                                                                            ttccTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATATTc                     >  1:2285714/4‑99 (MQ=255)
                                                                            ttACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                    >  1:2608504/1‑100 (MQ=255)
                                                                            ttACTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCg                    >  1:2308500/1‑100 (MQ=255)
                                                                               cTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATATTCTATc                 >  1:468014/1‑100 (MQ=255)
                                                                               cTGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATc                 <  1:1042733/100‑1 (MQ=255)
                                                                                tGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                <  1:2982475/100‑1 (MQ=255)
                                                                                tGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                >  1:3002775/1‑100 (MQ=255)
                                                                                tGTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                >  1:389321/1‑100 (MQ=255)
                                                                                 gTTTACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCg                >  1:909650/1‑99 (MQ=255)
                                                                                   ttACCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTa             >  2:290199/1‑100 (MQ=255)
                                                                                     aCCAACTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAg           >  1:2223078/1‑100 (MQ=255)
                                                                                        aaCTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGa        >  2:2323362/1‑100 (MQ=255)
                                                                                         aCTGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAt       >  2:468492/1‑100 (MQ=255)
                                                                                           tGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc     <  2:1742413/100‑1 (MQ=255)
                                                                                           tGGCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTc     <  2:2608492/100‑1 (MQ=255)
                                                                                            ggCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGAtt      >  1:430484/1‑98 (MQ=255)
                                                                                            ggCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa    <  2:2511109/100‑1 (MQ=255)
                                                                                            ggCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCa    <  1:652598/100‑1 (MQ=255)
                                                                                             gCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg   >  2:2689862/1‑100 (MQ=255)
                                                                                             gCTCCACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAg   <  2:1215765/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GGATGTACCGTATAACTATTTCCCGCACAATGATCCGCAAAATACACCGCGAGCGAGCTGGCGTAGTCACGGTAATTTACTGTTTACCAACTGGCTCAACTATTACGTCTACCAGATCACGCCATACGATCTACGGCACATGAATCCAACGCTGGATTAATCTTCTGTGATAGTCGATCGTTAAGCGATTCAGC  >  NC_000913/4215050‑4215243

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: