Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 577,681 (C)6→7 coding (352/462 nt) ybcT → endopeptidase‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913577,6751.C83.1% 544.1 / 115.8 213coding (346/462 nt)ybcTendopeptidase‑like protein
Reads supporting (aligned to +/- strand):  ref base . (14/22);  new base C (82/95);  total (96/117)
Fisher's exact test for biased strand distribution p-value = 4.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 236 total aligned reads displayed.

ATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAA  >  NC_000913/577583‑577766
                                                                                             |                                                                                            
 tgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑c                                                                                            >  2:841793/1‑93 (MQ=33)
 tgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑c                                                                                            >  1:2677993/1‑93 (MQ=33)
  gTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑cc                                                                                           >  1:1859538/1‑93 (MQ=33)
  gTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑cc                                                                                           <  2:1781979/93‑1 (MQ=34)
   ttGCCGCTAGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑ccc                                                                                          <  2:6990817/93‑1 (MQ=21)
    tGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑cccc                                                                                         <  1:596002/93‑1 (MQ=33)
     gCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCTT‑‑ccccc                                                                                        <  2:4771253/93‑1 (MQ=18)
      ccGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑cccccc                                                                                       >  2:4397087/1‑93 (MQ=33)
      ccGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑cccccc                                                                                       <  2:7249253/93‑1 (MQ=33)
       cGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCg                                                                                      >  1:1044275/1‑93 (MQ=33)
         cTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑cccccccga                                                                                    >  1:5107014/1‑90 (MQ=18)
         cTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑cccccccga                                                                                    >  2:1340589/1‑90 (MQ=18)
         cTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑cccccccga                                                                                    >  2:363973/1‑90 (MQ=18)
            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTg                                                                                  <  2:5088150/93‑1 (MQ=16)
            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTg                                                                                  <  1:4816997/93‑1 (MQ=16)
            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTg                                                                                  <  1:424712/93‑1 (MQ=16)
             tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTgg                                                                                 <  2:272151/93‑1 (MQ=16)
               gtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCa                                                                               <  1:6934914/93‑1 (MQ=16)
               gtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCa                                                                               >  1:6676865/1‑93 (MQ=16)
                tcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAg                                                                              >  2:3918052/1‑93 (MQ=16)
                tcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAg                                                                              >  1:2199235/1‑93 (MQ=16)
                tcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGa                                                                             <  2:6495441/93‑1 (MQ=33)
                  gtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGac                                                                            >  1:7092536/1‑93 (MQ=255)
                  gtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGac                                                                            >  1:4453341/1‑93 (MQ=255)
                   tcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGaca                                                                           <  2:3055528/93‑1 (MQ=16)
                   tcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGaca                                                                           >  1:7130715/1‑93 (MQ=255)
                   tcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGacac                                                                          <  2:2789842/93‑1 (MQ=33)
                    cgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGacac                                                                          >  2:5251712/1‑93 (MQ=255)
                     ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACAcc                                                                         <  2:3427538/93‑1 (MQ=16)
                     ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACAcc                                                                         >  1:5757795/1‑93 (MQ=255)
                      gTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCg                                                                        >  1:4033359/1‑93 (MQ=255)
                       ttGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                                       <  2:4902312/93‑1 (MQ=16)
                       ttGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                                       >  1:2569509/1‑93 (MQ=255)
                       ttGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                                       <  2:2178026/93‑1 (MQ=16)
                        tGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCt                                                                      >  1:358065/1‑93 (MQ=16)
                        tGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCt                                                                      <  1:6561963/93‑1 (MQ=16)
                         gCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTg                                                                     >  2:3948052/1‑93 (MQ=255)
                           acaTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAAc                                                                  <  1:658456/93‑1 (MQ=33)
                            caTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAAc                                                                  >  1:2368867/1‑93 (MQ=255)
                               cAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGa                                                              <  2:6386577/93‑1 (MQ=33)
                                 aaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                             >  2:3139811/1‑93 (MQ=255)
                                 aaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                             >  1:1084093/1‑93 (MQ=255)
                                 aaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                             >  1:5298727/1‑93 (MQ=255)
                                   gCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatta                                                           <  1:2124796/93‑1 (MQ=16)
                                   gCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatta                                                           <  2:1639787/93‑1 (MQ=16)
                                    cAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattat                                                          >  1:621719/1‑93 (MQ=255)
                                      gTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTc                                                       <  2:2018174/93‑1 (MQ=33)
                                       tCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTc                                                       >  2:5287662/1‑93 (MQ=255)
                                        cTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCa                                                      >  2:2496861/1‑93 (MQ=255)
                                          gtcagtcaGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAcc                                                    <  2:3559416/93‑1 (MQ=16)
                                          gtcagtcaGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAcc                                                    >  2:2556427/1‑93 (MQ=255)
                                             agtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                                 <  1:727116/93‑1 (MQ=16)
                                              gtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCa                                                <  1:1315801/93‑1 (MQ=255)
                                              gtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCa                                                <  1:4901581/93‑1 (MQ=255)
                                               tcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCag                                               >  2:419589/1‑93 (MQ=255)
                                                cagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCaga                                              <  1:4200586/93‑1 (MQ=255)
                                                cagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCaga                                              <  2:3587932/93‑1 (MQ=255)
                                                 agtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagag                                             >  1:7506484/1‑93 (MQ=255)
                                                   tGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagag                                           <  2:4847898/93‑1 (MQ=255)
                                                    gCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagaga                                          <  2:2300059/93‑1 (MQ=255)
                                                    gCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagaga                                          <  1:3731226/93‑1 (MQ=255)
                                                    gCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagaga                                          <  2:4376477/93‑1 (MQ=255)
                                                     cGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagagag                                         <  2:3623665/93‑1 (MQ=255)
                                                        gAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                      <  2:6060110/93‑1 (MQ=255)
                                                        gAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                      >  1:3820403/1‑93 (MQ=255)
                                                        gAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                      <  2:2535448/93‑1 (MQ=255)
                                                        gAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                      >  1:4940752/1‑93 (MQ=255)
                                                          aGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGa                                    >  2:6591413/1‑93 (MQ=255)
                                                           gCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGAt                                   >  2:3700271/1‑93 (MQ=255)
                                                            ccaccaCGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATc                                  >  1:7226419/1‑93 (MQ=255)
                                                             caccacGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCAc                                >  1:4669047/1‑93 (MQ=255)
                                                              accacGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACt                               <  2:6683417/93‑1 (MQ=255)
                                                                cacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTa                              <  2:5701820/93‑1 (MQ=255)
                                                                cacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTa                              >  1:1972006/1‑93 (MQ=255)
                                                                 acGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTAt                             >  2:793693/1‑93 (MQ=255)
                                                                 acGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTAt                             >  2:1495213/1‑93 (MQ=255)
                                                                 acGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTAt                             >  1:2087272/1‑93 (MQ=255)
                                                                   ggCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCa                          >  1:2268991/1‑93 (MQ=255)
                                                                    gCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCa                          <  1:7179302/93‑1 (MQ=255)
                                                                     ccTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaa                         >  2:318339/1‑93 (MQ=255)
                                                                       tCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaa                       >  2:5497031/1‑93 (MQ=255)
                                                                        ccGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaaa                      <  2:6203000/93‑1 (MQ=255)
                                                                           gCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaa                   <  1:5291859/93‑1 (MQ=255)
                                                                           gCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac                  <  1:2615230/93‑1 (MQ=255)
                                                                           gCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac                  <  1:806678/93‑1 (MQ=255)
                                                                            cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac                  <  1:3745024/93‑1 (MQ=255)
                                                                             gTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACt                 >  1:414581/1‑93 (MQ=255)
                                                                               ggATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTgg               <  1:3126149/93‑1 (MQ=255)
                                                                                gATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTgga              >  1:3473803/1‑93 (MQ=255)
                                                                                gATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTgga              <  2:1060927/93‑1 (MQ=255)
                                                                                 aTAATGCAGCCTCCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTgga              <  1:6153246/93‑1 (MQ=25)
                                                                                 aTAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaa             >  2:4782166/1‑93 (MQ=255)
                                                                                 aTAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaa             <  1:1685744/93‑1 (MQ=255)
                                                                                 aTAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTACGCAAAAACAACTggaa             <  2:5379491/93‑1 (MQ=255)
                                                                                    aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaagga          <  1:94026/93‑1 (MQ=255)
                                                                                    aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaagga          <  1:3888078/93‑1 (MQ=255)
                                                                                     tGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaaggaa         <  2:4626969/93‑1 (MQ=255)
                                                                                       cAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAAcc       <  1:555622/93‑1 (MQ=255)
                                                                                       cAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAAcc       <  1:1764030/93‑1 (MQ=255)
                                                                                        aGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAAccc      >  2:6817256/1‑93 (MQ=255)
                                                                                             |                                                                                            
ATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCT‑‑CCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAA  >  NC_000913/577583‑577766

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