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.C98.5% 583.0 / 1.2 195coding (346/462 nt)ybcTendopeptidase‑like protein
Reads supporting (aligned to +/- strand):  ref base . (0/3);  new base C (86/106);  total (86/109)
Fisher's exact test for biased strand distribution p-value = 2.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 224 total aligned reads displayed.

ATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCC  >  NC_000913/577583‑577762
                                                                                              |                                                                                       
atgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCt                                                                                          <  1:539502/93‑1 (MQ=34)
   ttGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCGCGGCCTCCGGCGTGGATAATGCAGCCTccc                                                                                      <  2:3417317/93‑1 (MQ=18)
   ttGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCCGCCTccc                                                                                      <  1:6358001/93‑1 (MQ=18)
   ttGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccc                                                                                      >  2:4635437/1‑93 (MQ=33)
     tccGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccccc                                                                                    <  1:5791301/92‑1 (MQ=18)
       cGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccccccc                                                                                  <  1:1206778/93‑2 (MQ=18)
        cgtGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccccg                                                                                 <  2:4277569/91‑3 (MQ=18)
          tGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccccgac                                                                               <  2:4990668/93‑5 (MQ=18)
            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTg                                                                              >  2:4200149/1‑93 (MQ=16)
            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTg                                                                              >  1:1796449/1‑93 (MQ=16)
              cgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGc                                                                            <  2:5476614/93‑1 (MQ=16)
                 cgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGa                                                                         >  2:6076305/1‑93 (MQ=255)
                 cgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGa                                                                         >  2:6380052/1‑93 (MQ=16)
                 cgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGa                                                                         >  2:6436278/1‑93 (MQ=16)
                    cgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCAGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGacac                                                                      <  2:2481080/93‑1 (MQ=255)
                      gTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCg                                                                    >  1:4280193/1‑93 (MQ=16)
                       ttGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                                   >  2:513919/1‑93 (MQ=255)
                       ttGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                                   >  1:1503117/1‑93 (MQ=16)
                        tGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCt                                                                  <  1:5374896/93‑1 (MQ=16)
                        tGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCt                                                                  >  1:4572472/1‑93 (MQ=255)
                         gCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTg                                                                 >  1:2970377/1‑93 (MQ=16)
                         gCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTg                                                                 <  1:6451771/93‑1 (MQ=16)
                         gCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTg                                                                 >  2:5861118/1‑93 (MQ=255)
                           acaTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGaa                                                               >  1:2356768/1‑93 (MQ=255)
                              tCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACgg                                                            <  2:1605324/93‑1 (MQ=16)
                                aaaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGa                                                          >  1:2742479/1‑93 (MQ=16)
                                 aaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                         <  1:6152845/93‑1 (MQ=16)
                                   gCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatta                                                       >  1:1171760/1‑93 (MQ=255)
                                    cAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattat                                                      <  1:2822893/93‑1 (MQ=16)
                                    cAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattat                                                      >  2:4021997/1‑93 (MQ=16)
                                     aGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattatt                                                     <  1:2612840/93‑1 (MQ=16)
                                     aGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattatt                                                     >  1:1576809/1‑93 (MQ=16)
                                      gTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTAttt                                                    >  2:2728422/1‑93 (MQ=16)
                                      gTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTAttt                                                    <  1:5710264/93‑1 (MQ=16)
                                      gTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTAttt                                                    >  1:5934468/1‑93 (MQ=255)
                                         tGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAc                                                 <  2:2435404/93‑1 (MQ=16)
                                         tGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAc                                                 >  1:2982984/1‑93 (MQ=255)
                                          gtcagtcaGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAcc                                                >  2:3160169/1‑93 (MQ=255)
                                          gtcagtcaGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAcc                                                >  2:3605678/1‑93 (MQ=255)
                                          gtcagtcaGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAcc                                                >  1:1672360/1‑93 (MQ=255)
                                           tcagtcagTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAccc                                               <  2:864367/93‑1 (MQ=16)
                                           tcagtcagTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAccc                                               >  2:3747404/1‑93 (MQ=255)
                                           tcagtcagTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAccc                                               <  2:329661/93‑1 (MQ=16)
                                            cagtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCGCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCt                                              <  2:938308/93‑1 (MQ=255)
                                              gtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCa                                            <  2:2283316/93‑1 (MQ=255)
                                               tcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCag                                           <  1:4500861/93‑1 (MQ=255)
                                               tcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCag                                           <  2:1161211/93‑1 (MQ=255)
                                               tcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCag                                           >  1:3581631/1‑93 (MQ=255)
                                                cagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGTCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCaga                                          >  2:4520350/1‑93 (MQ=255)
                                                cagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCaga                                          <  2:5357603/93‑1 (MQ=255)
                                                   tGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagag                                       <  2:4788208/93‑1 (MQ=255)
                                                      gTGAAACCACCCCGGCCTCCGGCGAGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAgg                                    <  1:4649092/93‑1 (MQ=14)
                                                       tGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGc                                   <  2:4063666/93‑1 (MQ=255)
                                                       tGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGc                                   <  1:5129823/93‑1 (MQ=255)
                                                        gacgCCACCACGGCCTCCGTCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCCCCCTCAGAGAGAGGCt                                  <  1:2437032/90‑1 (MQ=17)
                                                        gAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                  <  1:6206240/93‑1 (MQ=255)
                                                         aaGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGCGAGGCTg                                 <  1:6130794/93‑1 (MQ=255)
                                                         aaGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTg                                 >  1:2992852/1‑93 (MQ=255)
                                                             caccacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCa                             >  1:59971/1‑93 (MQ=255)
                                                             caccacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGAATATTTCACCCTCAGAGAGAGGCTGATCa                             >  2:3373709/1‑93 (MQ=255)
                                                              cccacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCAc                            <  1:1355684/92‑1 (MQ=255)
                                                              accacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCAc                            <  1:3163784/93‑1 (MQ=255)
                                                              accacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCAc                            <  1:4245764/93‑1 (MQ=255)
                                                              accacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCAc                            >  2:2538852/1‑93 (MQ=255)
                                                               ccacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACt                           >  1:56307/1‑93 (MQ=255)
                                                               ccacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACt                           <  1:2035247/93‑1 (MQ=255)
                                                               ccacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACt                           >  1:6282377/1‑93 (MQ=255)
                                                                cacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTa                          <  1:4931036/93‑1 (MQ=255)
                                                                cacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGACCACtt                          >  2:4609168/1‑92 (MQ=17)
                                                                 acGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTAt                         <  2:1451381/93‑1 (MQ=255)
                                                                  cGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATg                        <  1:5784561/93‑1 (MQ=255)
                                                                  cGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATg                        <  1:3766399/93‑1 (MQ=255)
                                                                   ggCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGc                       >  1:3535429/1‑93 (MQ=255)
                                                                    gCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCa                      >  2:1380509/1‑93 (MQ=255)
                                                                    gCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCa                      >  2:5012107/1‑93 (MQ=255)
                                                                    gCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCa                      <  2:843378/93‑1 (MQ=255)
                                                                    gCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCa                      >  2:2228594/1‑93 (MQ=255)
                                                                     ccTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaa                     <  1:2507457/93‑1 (MQ=255)
                                                                     ccTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaa                     <  2:2863371/93‑1 (MQ=255)
                                                                      cTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCCCCCTCAGAGAGAGGCTGATCACTATGCaaa                    <  2:1630793/93‑1 (MQ=255)
                                                                      cTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaa                    >  1:779147/1‑93 (MQ=255)
                                                                       tCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaa                   >  1:967569/1‑93 (MQ=255)
                                                                       tCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaa                   >  1:4016952/1‑93 (MQ=255)
                                                                        ccGGCGTGTATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTACGCaaaaa                  <  1:5161104/93‑1 (MQ=18)
                                                                        ccGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaaa                  <  1:4895249/93‑1 (MQ=255)
                                                                         cGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaac                 <  2:2933195/93‑1 (MQ=255)
                                                                           gCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaa               <  1:4408379/93‑1 (MQ=255)
                                                                            cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCGCCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac              <  1:2999482/93‑1 (MQ=255)
                                                                            cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac              <  2:5768429/93‑1 (MQ=255)
                                                                            cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac              >  2:3458588/1‑93 (MQ=255)
                                                                            cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac              <  1:2423050/93‑1 (MQ=255)
                                                                              tGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTg            <  2:2065828/93‑1 (MQ=255)
                                                                              tGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTg            <  1:1887993/93‑1 (MQ=255)
                                                                               ggATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTgg           >  2:6010892/1‑93 (MQ=255)
                                                                                 aTAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaa         >  2:5577011/1‑93 (MQ=255)
                                                                                 aTAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaa         <  2:900232/93‑1 (MQ=255)
                                                                                  tAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaag        >  2:1009101/1‑93 (MQ=255)
                                                                                      tcAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAAc    <  1:3205252/92‑1 (MQ=255)
                                                                                      gCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAAc    >  1:2542681/1‑93 (MQ=255)
                                                                                       cAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAAcc   <  1:3910048/93‑1 (MQ=255)
                                                                                              |                                                                                       
ATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCC  >  NC_000913/577583‑577762

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