Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 231,829 T→C G236G (GGT→GGC yafD → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913231,8290TC83.3% 628.6 / 120.4 293G236G (GGT→GGCyafDhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (23/26);  new base C (131/113);  total (154/139)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 293 total aligned reads displayed.

GCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAA  >  NC_000913/231739‑231918
                                                                                          |                                                                                          
  gggAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGa                                                                                        <  2:7296343/93‑1 (MQ=255)
   ggAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGaa                                                                                       <  1:4991863/93‑1 (MQ=255)
   ggAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGaa                                                                                       >  2:4035499/1‑93 (MQ=255)
      aaTGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGt                                                                                    <  2:2182216/93‑1 (MQ=255)
       aTGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTc                                                                                   >  1:3018944/1‑93 (MQ=255)
       aTGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTc                                                                                   <  2:4332287/93‑1 (MQ=255)
        tGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCa                                                                                  <  2:1478741/93‑1 (MQ=255)
         gTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAg                                                                                 <  1:3243343/93‑1 (MQ=255)
            gCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGa                                                                              <  1:5675065/93‑1 (MQ=255)
             cTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGaa                                                                             >  2:5887026/1‑93 (MQ=255)
             cTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGaa                                                                             >  2:6939156/1‑93 (MQ=255)
               gcgcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGc                                                                           >  2:7099703/1‑93 (MQ=255)
                 gcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCtt                                                                         <  1:3427108/93‑1 (MQ=255)
                 gcCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCtt                                                                         <  1:1449580/93‑1 (MQ=255)
                  ccAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTc                                                                        <  2:2352495/93‑1 (MQ=255)
                  ccAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTc                                                                        <  2:1361426/93‑1 (MQ=255)
                  ccAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTc                                                                        <  1:5077274/93‑1 (MQ=255)
                   cAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCt                                                                       <  1:5616835/93‑1 (MQ=255)
                     ggTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGt                                                                     >  2:2372980/1‑93 (MQ=255)
                     ggTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGt                                                                     >  1:3084262/1‑93 (MQ=255)
                      gTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTa                                                                    >  2:5425724/1‑93 (MQ=255)
                      gTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTa                                                                    <  2:6170476/93‑1 (MQ=255)
                      gTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTa                                                                    <  1:1127466/93‑1 (MQ=255)
                      gTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTa                                                                    >  1:2327738/1‑93 (MQ=255)
                      gTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTa                                                                    >  1:3614437/1‑93 (MQ=255)
                        gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACt                                                                  >  2:1052240/1‑93 (MQ=255)
                        gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACt                                                                  <  1:2199595/93‑1 (MQ=255)
                        gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACt                                                                  <  1:6493889/93‑1 (MQ=255)
                        gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACt                                                                  >  1:6154008/1‑93 (MQ=255)
                        gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACt                                                                  <  1:1642495/93‑1 (MQ=255)
                        gCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACt                                                                  >  1:5946468/1‑93 (MQ=255)
                          gTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTgg                                                                >  1:4113236/1‑93 (MQ=255)
                           ttttACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGt                                                               >  1:5546661/1‑93 (MQ=255)
                             ttACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTa                                                             <  2:3197777/93‑1 (MQ=255)
                             ttACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTa                                                             >  2:182583/1‑93 (MQ=255)
                             ttACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTa                                                             >  2:4293987/1‑93 (MQ=255)
                              tACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTAc                                                            >  2:1727109/1‑93 (MQ=255)
                              tACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTAc                                                            <  2:1181507/93‑1 (MQ=255)
                                ccGATGATCAGCGCCGTCGGGCGGTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgc                                                          >  1:3904573/1‑93 (MQ=255)
                                 cGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgcg                                                        >  1:2806384/1‑93 (MQ=255)
                                  gatgatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgcgc                                                       >  2:3797515/1‑93 (MQ=255)
                                   atgatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgcgcg                                                      >  2:2088777/1‑93 (MQ=255)
                                   atgatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTAcgcgcg                                                      >  2:7002694/1‑93 (MQ=255)
                                     gatCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCt                                                    >  2:5094062/1‑93 (MQ=255)
                                       tCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTc                                                  <  2:6689479/93‑1 (MQ=255)
                                        cAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTcc                                                 <  2:5402217/93‑1 (MQ=255)
                                          gcgcCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGa                                               <  1:4921079/93‑1 (MQ=255)
                                          gcgcCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGa                                               <  1:1053524/93‑1 (MQ=255)
                                          gcgcCGTCGGGCGTGTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGa                                               >  2:6442788/1‑93 (MQ=255)
                                             ccGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATca                                            <  1:6692127/93‑1 (MQ=255)
                                             ccGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATca                                            >  1:5227709/1‑93 (MQ=255)
                                              cGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATcac                                           >  1:4894928/1‑93 (MQ=255)
                                               gTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATcaca                                          >  1:4323561/1‑93 (MQ=255)
                                                tCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTGCTGGTTACGCGCGCTTCCGATCACaa                                         <  1:6045931/93‑1 (MQ=255)
                                                tCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACaa                                         <  2:2475820/93‑1 (MQ=255)
                                                tCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATcacac                                         <  2:7407289/93‑2 (MQ=255)
                                                tCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACaa                                         <  2:852215/93‑1 (MQ=255)
                                                  gggCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGGCGCTTCCGATCACAAt                                        >  1:2357476/1‑93 (MQ=255)
                                                  gggCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATc                                       <  1:5472530/93‑1 (MQ=255)
                                                  gggCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATc                                       <  2:2683907/93‑1 (MQ=255)
                                                  gggCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATc                                       <  1:206151/93‑1 (MQ=255)
                                                  gggCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCACTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATc                                       >  1:794963/1‑93 (MQ=255)
                                                      gTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCt                                   >  1:2988275/1‑93 (MQ=255)
                                                       tttGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTa                                  <  2:2410839/93‑1 (MQ=255)
                                                       tttGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTa                                  >  1:6386207/1‑93 (MQ=255)
                                                            tCGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGt                             >  2:2576242/1‑93 (MQ=255)
                                                             cGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGtt                            <  1:4667945/93‑1 (MQ=255)
                                                             cGCCCGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGtt                            >  2:5200858/1‑93 (MQ=255)
                                                                ccGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGaa                         >  2:2608551/1‑93 (MQ=255)
                                                                 cGCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAAt                        >  2:716086/1‑93 (MQ=255)
                                                                  gCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAAtt                       <  2:2784881/93‑1 (MQ=255)
                                                                  gCTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAAtt                       >  1:7314719/1‑93 (MQ=255)
                                                                   cTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTc                      <  2:1674974/93‑1 (MQ=255)
                                                                   cTCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTc                      <  1:6811164/93‑1 (MQ=255)
                                                                    tCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCa                     >  1:4412737/1‑93 (MQ=255)
                                                                    tCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCa                     <  1:470023/93‑1 (MQ=255)
                                                                    tCGATTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCa                     >  2:6342222/1‑93 (MQ=255)
                                                                     cgtttttGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAAGCCGCTACTCGTTGAATTCAg                    <  1:6280785/90‑1 (MQ=255)
                                                                       aTTTTGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTc                  >  2:2985952/1‑93 (MQ=255)
                                                                        ttttGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTcc                 <  1:4480639/93‑1 (MQ=255)
                                                                        ttttGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTcc                 <  2:4432003/93‑1 (MQ=255)
                                                                        ttttGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTcc                 <  2:7227623/93‑1 (MQ=255)
                                                                        ttttGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTcc                 <  2:6793745/93‑1 (MQ=255)
                                                                         tttGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTccc                >  2:3508840/1‑93 (MQ=255)
                                                                          ttGTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCg               >  2:5650199/1‑93 (MQ=255)
                                                                            gTTTTCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGc             >  1:3477049/1‑93 (MQ=255)
                                                                                tCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGc         <  2:1888805/93‑1 (MQ=255)
                                                                                tCTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGc         <  2:274385/93‑1 (MQ=255)
                                                                                ccTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGc         <  1:242886/92‑1 (MQ=255)
                                                                                 cTACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGcc        <  2:4859366/93‑1 (MQ=255)
                                                                                  tACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCt       >  2:4022720/1‑93 (MQ=255)
                                                                                  tACCGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCt       >  1:3981058/1‑93 (MQ=255)
                                                                                    ccGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGa     <  1:2820577/93‑1 (MQ=255)
                                                                                    ccGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGa     >  2:5846454/1‑93 (MQ=255)
                                                                                     cGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGat    <  1:2878073/93‑1 (MQ=255)
                                                                                     cGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGat    >  2:5595131/1‑93 (MQ=255)
                                                                                     cGTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGat    >  2:4671866/1‑93 (MQ=255)
                                                                                      gTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGata   <  1:1359120/93‑1 (MQ=255)
                                                                                      gTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGata   <  2:5713914/93‑1 (MQ=255)
                                                                                      gTGGCCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGata   <  2:3404295/93‑1 (MQ=255)
                                                                                          |                                                                                          
GCGGGAAATGTCGCTGCGCCAGGTGCGTTTTACCGATGATCAGCGCCGTCGGGCGTTTGGTCGCCCGCTCGATTTTGTTTTCTACCGTGGTCTGAACGTCAGTGAAGCTTCTGTACTGGTTACGCGCGCTTCCGATCACAATCCGCTACTCGTTGAATTCAGTCCCGGCAAGCCTGATAA  >  NC_000913/231739‑231918

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