Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 418,805 A→G intergenic (+397/‑10) phoR → / → brnQ positive and negative sensor protein for pho regulon/branched chain amino acid transport system II carrier protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913418,8050AG90.8% 720.9 / 54.2 262intergenic (+397/‑10)phoR/brnQpositive and negative sensor protein for pho regulon/branched chain amino acid transport system II carrier protein
Reads supporting (aligned to +/- strand):  ref base A (11/13);  new base G (105/133);  total (116/146)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

Only 100 of 262 total aligned reads displayed.

TATTTCCCGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACAT  >  NC_000913/418715‑418894
                                                                                          |                                                                                         
 aTTTCCCGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTa                                                                                        <  1:2959042/93‑1 (MQ=255)
  tttCCCGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTAt                                                                                       >  1:1028426/1‑93 (MQ=255)
    tCCCGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGa                                                                                     <  1:27996/93‑1 (MQ=255)
     cccGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGAt                                                                                    <  2:4470757/93‑1 (MQ=255)
      ccGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGAtt                                                                                   <  1:1621253/93‑1 (MQ=255)
       cGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGAttt                                                                                  >  1:6440700/1‑93 (MQ=255)
        gCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTa                                                                                 <  2:6021884/93‑1 (MQ=255)
        gCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTa                                                                                 <  2:1729661/93‑1 (MQ=255)
         cTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACGGGCGGTATGATTTAt                                                                                <  2:5512343/93‑1 (MQ=255)
           aaaaTTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGa                                                                              <  2:208622/93‑1 (MQ=255)
           aaaaTTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGa                                                                              >  2:6379394/1‑93 (MQ=255)
           aaaaTTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGa                                                                              >  1:3370627/1‑93 (MQ=255)
           aaaaTTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGa                                                                              >  2:7033487/1‑93 (MQ=255)
             aaTTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGAcc                                                                            <  2:5418707/93‑1 (MQ=255)
             aaTTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGAcc                                                                            >  2:1514018/1‑93 (MQ=255)
               ttGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCa                                                                          <  1:5778460/93‑1 (MQ=255)
                tGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCAt                                                                         >  2:64987/1‑93 (MQ=255)
                tGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCAt                                                                         >  1:764011/1‑93 (MQ=255)
                 gTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATc                                                                        >  1:5905344/1‑93 (MQ=255)
                 gTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATc                                                                        >  1:5658238/1‑93 (MQ=255)
                 gTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATc                                                                        <  1:5075544/93‑1 (MQ=255)
                  tttAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCa                                                                       >  2:4297235/1‑93 (MQ=255)
                  tttAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCa                                                                       >  1:5370641/1‑93 (MQ=255)
                     aaaTATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAAtt                                                                    >  1:2599061/1‑93 (MQ=255)
                     aaaTATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAAtt                                                                    <  1:2467837/93‑1 (MQ=255)
                     aaaTATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAAtt                                                                    >  1:5528325/1‑93 (MQ=255)
                      aaTATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTa                                                                   >  1:5550417/1‑93 (MQ=255)
                       atatACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTaa                                                                  <  2:4940525/93‑1 (MQ=255)
                       atatACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTaa                                                                  >  2:6591712/1‑93 (MQ=255)
                         ataCCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGa                                                                >  1:5951163/1‑93 (MQ=255)
                          taCCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGAt                                                               >  2:3034664/1‑93 (MQ=255)
                          taCCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGAt                                                               >  2:5234507/1‑93 (MQ=255)
                           aCCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATc                                                              <  2:6272587/93‑1 (MQ=255)
                           aCCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATc                                                              <  1:4985589/93‑1 (MQ=255)
                            ccGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATcg                                                             <  1:6519852/93‑1 (MQ=255)
                            ccGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATcg                                                             >  2:3536956/1‑93 (MQ=255)
                             cGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATcgc                                                            >  2:5085766/1‑93 (MQ=255)
                             cGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATcgc                                                            >  2:3412426/1‑93 (MQ=255)
                             cGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATcgc                                                            <  2:912368/93‑1 (MQ=255)
                              gCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATcgcg                                                           >  1:4228783/1‑93 (MQ=255)
                              gCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATcgcg                                                           >  2:1022284/1‑93 (MQ=255)
                              gCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATcgcg                                                           >  1:7470370/1‑93 (MQ=255)
                               cTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATcgcgc                                                          <  1:3867573/93‑1 (MQ=255)
                                 gTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGa                                                        <  2:573898/93‑1 (MQ=255)
                                 gTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGa                                                        <  2:5890584/93‑1 (MQ=255)
                                   atcatcCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGata                                                      >  1:3212743/1‑93 (MQ=255)
                                    tcatcCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGatat                                                     <  1:6531401/93‑1 (MQ=255)
                                     catcCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATatc                                                    <  2:2661394/93‑1 (MQ=255)
                                      atcCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATatca                                                   <  2:2741912/93‑1 (MQ=255)
                                      atcCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATatca                                                   <  2:1370818/93‑1 (MQ=255)
                                        ccccAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTCTGATTTATGACCCATCAATTAAGATCGCGCGATatcatc                                                 <  1:6217591/93‑1 (MQ=255)
                                          ccAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGc                                               >  2:3891670/1‑93 (MQ=255)
                                          ccAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGc                                               <  2:1625439/93‑1 (MQ=255)
                                           cAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGct                                              >  2:3843433/1‑93 (MQ=255)
                                            aGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGcac                                             >  1:837397/1‑91 (MQ=255)
                                             gggATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGctct                                            <  2:3628387/93‑1 (MQ=255)
                                             gggATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGctct                                            <  1:4057864/93‑1 (MQ=255)
                                             gggATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGctct                                            <  1:3505119/93‑1 (MQ=255)
                                              ggATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTg                                           >  1:3742779/1‑93 (MQ=255)
                                              ggATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTg                                           <  1:3070005/93‑1 (MQ=255)
                                              ggATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTg                                           <  1:1648071/93‑1 (MQ=255)
                                               gATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTgg                                          <  2:1455269/93‑1 (MQ=255)
                                                 ttGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGc                                        <  1:2725981/93‑1 (MQ=255)
                                                     cacaAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTa                                    <  2:1941536/93‑1 (MQ=255)
                                                     cacaAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTa                                    <  2:1740958/93‑1 (MQ=255)
                                                     cacaAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTa                                    <  2:4075237/93‑1 (MQ=255)
                                                     cacaAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTa                                    <  1:7267088/93‑1 (MQ=255)
                                                      acaAAAATTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTAt                                   <  2:7270641/93‑1 (MQ=255)
                                                        aaaaaTTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGa                                 <  1:3963797/93‑1 (MQ=255)
                                                        aaaaaTTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGa                                 <  2:2375352/93‑1 (MQ=255)
                                                        aaaaaTTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGa                                 <  1:6287908/93‑1 (MQ=255)
                                                        aaaaaTTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGa                                 >  2:1862768/1‑93 (MQ=255)
                                                        aaaaaTTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGa                                 <  1:4295856/93‑1 (MQ=255)
                                                          aaaTTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACa                               <  1:5267912/93‑1 (MQ=255)
                                                           aaTTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACAt                              <  2:641739/93‑1 (MQ=255)
                                                            aTTTAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACAtt                             <  2:4344582/93‑1 (MQ=255)
                                                               tAACGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGc                          >  1:7187206/1‑93 (MQ=255)
                                                                aaCGTNACAACACCACATCCACAGGCGGTATGATTTATGACCCNNNNATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCg                         <  2:3722193/93‑1 (MQ=255)
                                                                 aCGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCgt                        <  2:5129961/93‑1 (MQ=255)
                                                                  cGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCgtt                       <  2:2935827/93‑1 (MQ=255)
                                                                  cGTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCgtt                       <  1:4779119/93‑1 (MQ=255)
                                                                   gTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCgttg                      <  2:1194902/93‑1 (MQ=255)
                                                                   gTTACAACACCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCgttg                      >  2:2563888/1‑93 (MQ=255)
                                                                      acaacaCCACATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTtc                   >  1:6783658/1‑93 (MQ=255)
                                                                           accacATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCgg              >  2:4206388/1‑93 (MQ=255)
                                                                           accacATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCgg              >  2:2226822/1‑93 (MQ=255)
                                                                           accacATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCgg              <  1:6950335/93‑1 (MQ=255)
                                                                           accacATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCgg              <  1:2116478/93‑1 (MQ=255)
                                                                            ccacATCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGgc             <  1:1940381/93‑1 (MQ=255)
                                                                             cacaTCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGgcg            >  2:5809187/1‑93 (MQ=255)
                                                                              acaTCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGgcgc           <  2:1807254/93‑1 (MQ=255)
                                                                               caTCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCa          <  1:4606668/93‑1 (MQ=255)
                                                                               caTCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCa          >  2:776970/1‑93 (MQ=255)
                                                                               caTCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCa          <  1:67018/93‑1 (MQ=255)
                                                                                aTCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAg         <  1:4388469/93‑1 (MQ=255)
                                                                                 tCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAgg        >  2:6757390/1‑93 (MQ=255)
                                                                                 tCCACAGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAgg        >  1:1320799/1‑93 (MQ=255)
                                                                                   cacaGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTa      <  1:2280184/93‑1 (MQ=255)
                                                                                    acaGGCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTaa     >  2:4499291/1‑93 (MQ=255)
                                                                                       ggCGGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACat  >  1:2891720/1‑93 (MQ=255)
                                                                                          |                                                                                         
TATTTCCCGCTAAAATTGTTTAAATATACCGCTGTATCATCCCCAGGGATTGGCACAAAAATTTAACGTTACAACACCACATCCACAGGCAGTATGATTTATGACCCATCAATTAAGATCGCGCGATATCATCGCTCTGGGCTTTATGACATTTGCGTTGTTCGTCGGCGCAGGTAACAT  >  NC_000913/418715‑418894

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