Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,512,256 A→G S664P (TCC→CCC)  fecA ← FecA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,512,2560AG88.4% 777.5 / 79.9 278S664P (TCC→CCC) fecAFecA
Reads supporting (aligned to +/- strand):  ref base A (16/16);  new base G (101/143);  total (117/161)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.16e-01

Only 100 of 294 total aligned reads displayed.

CGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCG  >  NC_000913/4512164‑4512348
                                                                                            |                                                                                            
   aaaCTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGaa                                                                                           >  1:6112964/1‑93 (MQ=255)
    aaCTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAAt                                                                                          <  2:709983/93‑1 (MQ=255)
      ctggctggACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGGATgg                                                                                        >  1:796643/1‑93 (MQ=255)
      ctggctggACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATgg                                                                                        <  1:2584163/93‑1 (MQ=255)
       tggctggACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGt                                                                                       >  1:3092228/1‑93 (MQ=255)
        ggctggACTGGAAATCGCTGGTCAGATTGAACGTCCAGATTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTa                                                                                      >  1:3239371/1‑93 (MQ=255)
         gctggACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTAc                                                                                     <  1:3658495/93‑1 (MQ=255)
            ggactggaAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGGATGGTACCAg                                                                                  >  2:1286148/1‑93 (MQ=255)
            ggactggaAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGGATGGTACCAg                                                                                  >  1:5957311/1‑93 (MQ=255)
            ggactggaAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAg                                                                                  <  2:5357936/93‑1 (MQ=255)
             gactggaAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGa                                                                                 <  1:2523474/93‑1 (MQ=255)
             gactggaAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGa                                                                                 <  1:686531/93‑1 (MQ=255)
               ctggaAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGAtt                                                                               >  2:5527499/1‑93 (MQ=255)
                tggaAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTg                                                                              >  1:4623293/1‑93 (MQ=255)
                tggaAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTg                                                                              <  1:3133929/93‑1 (MQ=255)
                   aaaTCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGACg                                                                           >  2:1536166/1‑93 (MQ=255)
                      tCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAg                                                                        <  1:5779158/93‑1 (MQ=255)
                       cGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAgg                                                                       <  2:5924023/93‑1 (MQ=255)
                        gCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAggg                                                                      >  2:3485097/1‑92 (MQ=255)
                        gCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGgt                                                                      <  1:1573028/93‑1 (MQ=255)
                         cTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGgtg                                                                     <  1:6604745/93‑1 (MQ=255)
                          tGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGgtgt                                                                    <  2:1052581/93‑1 (MQ=255)
                            ttCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCg                                                                  >  2:2013155/1‑93 (MQ=255)
                             tCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGc                                                                 <  2:2094285/93‑1 (MQ=255)
                               aGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGGATGGTACCAGATTGCCGTAGGTGTCGCCt                                                               >  2:811120/1‑93 (MQ=255)
                                gATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCtt                                                              <  1:3436539/93‑1 (MQ=255)
                                 aTTGAACGTCCAGTTGCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCttt                                                             <  1:3270116/93‑1 (MQ=255)
                                   tGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTtct                                                           <  1:2887632/93‑1 (MQ=255)
                                   tGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTtct                                                           >  1:4539513/1‑93 (MQ=255)
                                     aaCGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcg                                                         >  2:630798/1‑93 (MQ=255)
                                     aaCGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcg                                                         >  2:4588301/1‑93 (MQ=255)
                                     aaCGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcg                                                         <  1:4323452/93‑1 (MQ=255)
                                      aCGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcgc                                                        >  1:4452212/1‑93 (MQ=255)
                                      aCGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcgc                                                        <  1:6553450/93‑1 (MQ=255)
                                      aCGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcgc                                                        <  1:4943671/93‑1 (MQ=255)
                                       cGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcgcg                                                       <  2:3393400/93‑1 (MQ=255)
                                       cGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTcgcg                                                       >  1:5244381/1‑93 (MQ=255)
                                         tCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGa                                                     <  2:1427638/93‑1 (MQ=255)
                                         gccAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGa                                                     <  1:923128/92‑1 (MQ=255)
                                           cAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGAtt                                                   <  2:4259979/93‑1 (MQ=255)
                                            aGTTTCCTGGCTTGTAGTCCACGCCCAGCCTGCCCTTTTGTTTCGGGGGGGATGGTACCCGATTGGCGTAGGTGTCGCCTTTTTCGCGGAttt                                                  >  2:3864707/1‑93 (MQ=255)
                                              tttCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTcc                                                <  2:4206549/93‑1 (MQ=255)
                                              tttCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTcc                                                >  2:3188099/1‑93 (MQ=255)
                                               ttCCGGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCcg                                               <  2:588781/93‑1 (MQ=255)
                                                  cTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGt                                            >  2:6396551/1‑93 (MQ=255)
                                                  cTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGt                                            <  1:5552961/93‑1 (MQ=255)
                                                  cTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGt                                            <  2:2683146/93‑1 (MQ=255)
                                                   tGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGtt                                           <  1:5165768/93‑1 (MQ=255)
                                                   tGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGtt                                           >  1:5927226/1‑93 (MQ=255)
                                                    ggCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTc                                          >  2:3001840/1‑93 (MQ=255)
                                                     gCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTca                                         <  1:5679482/93‑1 (MQ=255)
                                                         gTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATa                                     <  1:4554080/93‑1 (MQ=255)
                                                         gTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATa                                     <  1:5901357/93‑1 (MQ=255)
                                                          tAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATAc                                    <  2:1984663/93‑1 (MQ=255)
                                                          tAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATAc                                    >  2:3775752/1‑93 (MQ=255)
                                                           aGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACg                                   >  1:4126366/1‑93 (MQ=255)
                                                            gTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGc                                  >  2:5050753/1‑93 (MQ=255)
                                                            gTCCACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGc                                  >  2:5812645/1‑93 (MQ=255)
                                                              ccACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCAt                                <  1:2153703/93‑1 (MQ=255)
                                                              ccACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCAt                                <  1:1614976/93‑1 (MQ=255)
                                                               cACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATa                               <  1:6267289/93‑1 (MQ=255)
                                                               cACGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATa                               >  1:1434363/1‑93 (MQ=255)
                                                                aCGCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAg                              <  1:1743408/93‑1 (MQ=255)
                                                                aCGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAg                              <  1:3810450/93‑1 (MQ=255)
                                                                  gCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCt                            <  2:5604729/93‑1 (MQ=255)
                                                                  gCCCAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCt                            <  2:6132929/93‑1 (MQ=255)
                                                                   cccAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTg                           <  2:5166878/93‑1 (MQ=255)
                                                                   cccAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTg                           >  1:5593896/1‑93 (MQ=255)
                                                                    ccAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTgg                          <  1:233353/93‑1 (MQ=255)
                                                                     cAGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGc                         <  2:1154630/93‑1 (MQ=255)
                                                                      aGCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCg                        <  2:4415643/93‑1 (MQ=255)
                                                                       gCGTGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGt                       >  1:3388560/1‑93 (MQ=255)
                                                                       gCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGt                       >  1:3939912/1‑93 (MQ=255)
                                                                          tGCCTTTATGTTTCGGGGGTAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGa                    <  1:1348647/93‑1 (MQ=255)
                                                                          tGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGa                    <  2:1720411/93‑1 (MQ=255)
                                                                          tGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGa                    >  1:4690025/1‑93 (MQ=255)
                                                                          tGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGa                    >  2:4772394/1‑93 (MQ=255)
                                                                          tGCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGa                    <  2:3110806/93‑1 (MQ=255)
                                                                           gCCTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGAt                   <  1:3326740/93‑1 (MQ=255)
                                                                            ccTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTATGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATg                  <  1:1539119/93‑1 (MQ=255)
                                                                            ccTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACTTACGCATAGCTGGCTTAGATg                  <  1:4226336/93‑1 (MQ=255)
                                                                            ccTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATg                  >  2:3315955/1‑93 (MQ=255)
                                                                            ccTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATg                  >  2:4948156/1‑93 (MQ=255)
                                                                             cTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATgg                 >  1:5527242/1‑93 (MQ=255)
                                                                             cTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATgg                 <  1:4813205/93‑1 (MQ=255)
                                                                             cTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATgg                 <  1:1223485/93‑1 (MQ=255)
                                                                             cTTTATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATgg                 <  1:3048814/93‑1 (MQ=255)
                                                                              tttATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGa                >  1:3496806/1‑93 (MQ=255)
                                                                               ttATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGaa               <  1:3714087/93‑1 (MQ=255)
                                                                               ttATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGaa               >  2:2536026/1‑93 (MQ=255)
                                                                               ttATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGaa               <  2:1805613/93‑1 (MQ=255)
                                                                                tATGTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGaaa              <  1:5442323/93‑1 (MQ=255)
                                                                                   gTTTCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGt           <  1:1505196/93‑1 (MQ=255)
                                                                                    tttCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGGTTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGtt          >  1:1997829/1‑93 (MQ=255)
                                                                                      tCGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTAt        <  2:2626720/93‑1 (MQ=255)
                                                                                       cGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATc       <  2:4711258/93‑1 (MQ=255)
                                                                                       cGGGGGGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATc       <  1:1086757/93‑1 (MQ=255)
                                                                                         ggggGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCaa     >  1:1680988/1‑93 (MQ=255)
                                                                                         gggAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCaa     <  2:4039762/93‑1 (MQ=255)
                                                                                            agaATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCg  <  2:21774/93‑1 (MQ=255)
                                                                                            |                                                                                            
CGCAAACTGGCTGGACTGGAAATCGCTGTTCAGATTGAACGTCCAGTTTCCTGGCTTGTAGTCCACGCCCAGCGTGCCTTTATGTTTCGGGGAGAATGGTACCAGATTGCCGTAGGTGTCGCCTTTCTCGCGGATTTCCGCGTTCACATACGCATAGCTGGCGTAGATGGAAACGTTATCAAGCG  >  NC_000913/4512164‑4512348

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑