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.C97.0% 861.2 / 17.9 295coding (346/462 nt)ybcTendopeptidase‑like protein
Reads supporting (aligned to +/- strand):  ref base . (4/5);  new base C (123/163);  total (127/168)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 347 total aligned reads displayed.

ATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCC  >  NC_000913/577583‑577762
                                                                                             |                                                                                        
atgTTGCCGCTGGTCGTCGTCGGTTGCATATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCt                                                                                           >  2:6707136/1‑93 (MQ=18)
atgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCt                                                                                           >  1:6560271/1‑93 (MQ=33)
 tgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTc                                                                                         >  2:8040633/1‑93 (MQ=33)
   ttGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccc                                                                                      >  1:2342322/1‑93 (MQ=33)
      gcGCTGGTCGACGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccc                                                                                   <  1:2655560/92‑1 (MQ=18)
      ccGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccc                                                                                   >  2:394332/1‑93 (MQ=33)
       cGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTccccccc                                                                                  <  1:8460640/93‑2 (MQ=18)
       cGCTGGCAGGCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCGccccccc                                                                                  <  2:7848057/93‑2 (MQ=11)
        gCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTAAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTcccccccg                                                                                 >  2:2338997/1‑91 (MQ=16)
            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTg                                                                              >  1:1562981/1‑93 (MQ=16)
            gtcgtcgtcGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGCGGTTAATGCAGCCTCCCCCCCGACTg                                                                              >  2:5978967/1‑93 (MQ=11)
             tcgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTgg                                                                             >  2:5323428/1‑93 (MQ=16)
              cgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGc                                                                            <  1:7377847/93‑1 (MQ=16)
              cgtcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGc                                                                            <  2:8178107/93‑1 (MQ=16)
                tcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAg                                                                          >  2:472061/1‑93 (MQ=16)
                    cgGTTGCAGATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGacac                                                                      <  2:2718821/93‑1 (MQ=255)
                    cgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGacac                                                                      <  2:6417690/93‑1 (MQ=16)
                    cgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGacac                                                                      <  1:8065406/93‑1 (MQ=16)
                     ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACAcc                                                                     <  1:3722776/93‑1 (MQ=16)
                     ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACAcc                                                                     >  2:5472505/1‑93 (MQ=16)
                     ggTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACAcc                                                                     >  2:3860687/1‑93 (MQ=16)
                      gTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCg                                                                    <  2:8622856/93‑1 (MQ=16)
                      gTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACGGGCAGacacgg                                                                    >  1:1661887/1‑91 (MQ=11)
                       ttGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                                   <  2:625242/93‑1 (MQ=16)
                       ttGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                                   <  2:7055921/93‑1 (MQ=16)
                       tggCACATCAAAGCAGTCTGTCAGTCAGTTCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGc                                                                   <  1:7439528/91‑1 (MQ=11)
                        tGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCt                                                                  >  1:2760486/1‑93 (MQ=16)
                        tGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCt                                                                  <  1:8129324/93‑1 (MQ=16)
                          cacaTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGa                                                                <  1:8719746/93‑1 (MQ=16)
                           acaTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGaa                                                               >  1:1819624/1‑93 (MQ=255)
                           acaTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGaa                                                               <  1:5760939/93‑1 (MQ=16)
                            caTCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTAGCAGACACCACTGAAc                                                              >  2:3222327/1‑93 (MQ=11)
                              tCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACgg                                                            <  2:4937627/93‑1 (MQ=16)
                               taaaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTTGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACggg                                                           <  1:7109009/92‑1 (MQ=11)
                                 aaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                         >  1:4827521/1‑93 (MQ=255)
                                 aaGCAGTCTGTCAGGCAGTGCGTGAAGCCCCCTCGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGat                                                         <  1:6795613/93‑1 (MQ=11)
                                  aGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatt                                                        >  2:1062136/1‑93 (MQ=16)
                                  aGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatt                                                        >  1:4680253/1‑93 (MQ=255)
                                  aGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatt                                                        >  2:4051688/1‑93 (MQ=255)
                                   gCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGatta                                                       <  1:1147426/93‑1 (MQ=16)
                                    cAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattat                                                      <  1:3786811/93‑1 (MQ=16)
                                    cAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattat                                                      <  2:7430919/93‑1 (MQ=16)
                                    cAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattat                                                      <  2:5213823/93‑1 (MQ=16)
                                     aGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGattatt                                                     <  1:2259982/93‑1 (MQ=16)
                                     aGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCCCCCCCCGACTGGCAGACACCGCTGAACGGGATTAttt                                                    <  2:2956514/93‑1 (MQ=18)
                                      gTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTAttt                                                    >  2:4532796/1‑93 (MQ=255)
                                       gcTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTc                                                   <  1:4250889/92‑1 (MQ=255)
                                        cTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCa                                                  >  1:2523413/1‑93 (MQ=255)
                                         tGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAc                                                 <  1:4384223/93‑1 (MQ=16)
                                         tGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAc                                                 <  1:3391104/93‑1 (MQ=16)
                                           tcagtcagTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAccc                                               <  1:6303110/93‑1 (MQ=16)
                                           tcagtcagTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAccc                                               >  2:2220597/1‑93 (MQ=255)
                                            cagtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCt                                              >  1:8229235/1‑93 (MQ=16)
                                             agtcagtGCGTGAAGCCCCCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                             <  1:2828810/93‑1 (MQ=255)
                                             agtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTc                                             <  1:189836/93‑1 (MQ=16)
                                             agtcagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTAACCCTc                                             >  2:1312746/1‑93 (MQ=255)
                                                cagtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCaga                                          <  1:1038475/93‑1 (MQ=255)
                                                  gtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagaga                                        <  1:6337655/93‑1 (MQ=255)
                                                    gCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagaga                                      <  2:1635454/93‑1 (MQ=255)
                                                    gCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagaga                                      <  1:3630764/93‑1 (MQ=255)
                                                     cGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagagag                                     <  2:2001521/93‑1 (MQ=255)
                                                     cGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATAATTTCACCCTCagagagag                                     <  1:8535200/93‑1 (MQ=255)
                                                      gTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTCTTTCACCCTCAGAGAGAgg                                    >  1:7700626/1‑93 (MQ=255)
                                                      gTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAgg                                    <  1:64620/93‑1 (MQ=255)
                                                       tGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagagagagtc                                   <  2:8300878/93‑3 (MQ=255)
                                                       tGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGc                                   <  1:6690952/93‑1 (MQ=255)
                                                        gAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                  <  2:2551819/93‑1 (MQ=255)
                                                        gAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                  <  1:6003481/93‑1 (MQ=255)
                                                        gAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                  <  1:397513/93‑1 (MQ=255)
                                                        gAAGCCACCACGGCCTCCGGCATGCATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCt                                  <  2:2898197/93‑1 (MQ=17)
                                                          aGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGGGAGGCTGa                                >  2:1825937/1‑93 (MQ=255)
                                                           gCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGAt                               >  2:881963/1‑93 (MQ=255)
                                                           gCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGAt                               <  2:3310433/93‑1 (MQ=255)
                                                            ccaccaCGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATc                              <  1:5883453/93‑1 (MQ=255)
                                                            ccaccaCGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATc                              >  1:5605908/1‑93 (MQ=255)
                                                              accacGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCAc                            >  2:5764438/1‑93 (MQ=255)
                                                                 acGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTAt                         <  2:837837/93‑1 (MQ=255)
                                                                 acGGCCTCCGGCGCGGATAATGCAGACTCCCCCCCGACCGGCCGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGACTGACCACTAt                         >  2:1530590/1‑93 (MQ=17)
                                                                   ggCCTCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGc                       <  2:3696145/93‑1 (MQ=255)
                                                                       tCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaa                   >  2:3616569/1‑93 (MQ=255)
                                                                       tCCGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaa                   >  1:5150095/1‑93 (MQ=255)
                                                                        ccGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaaa                  >  1:1300017/1‑93 (MQ=255)
                                                                        ccGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaaa                  >  1:8578873/1‑93 (MQ=255)
                                                                        ccGGCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaaaaa                  <  1:3323758/93‑1 (MQ=255)
                                                                          ggCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaaca                >  2:5691668/1‑93 (MQ=255)
                                                                           gCGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaa               <  1:4429481/93‑1 (MQ=255)
                                                                            cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac              >  2:4739284/1‑93 (MQ=255)
                                                                            cGTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac              >  1:8000016/1‑93 (MQ=255)
                                                                             gTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACt             <  2:6793709/93‑1 (MQ=255)
                                                                             gTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACt             <  2:147829/93‑1 (MQ=255)
                                                                             gTGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACt             >  1:6887070/1‑93 (MQ=255)
                                                                              tGGATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTg            <  1:5038415/93‑1 (MQ=255)
                                                                               ggATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTgg           <  2:34890/93‑1 (MQ=255)
                                                                                gATAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTgga          >  2:7167157/1‑93 (MQ=255)
                                                                                 aTAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGATAGAGGCTGATCACTATGCAAAAACAACTggaa         >  2:7379840/1‑93 (MQ=255)
                                                                                 aTAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaa         <  2:4449202/93‑1 (MQ=255)
                                                                                 aTAATGCAGCCACCACCCCGACTGGCAGACACCGCAGAACGGGATTATTACAACCTCAAAGAGAGAAAGATCACTATGCAAAAACAACTggaa         <  1:7835121/93‑1 (MQ=18)
                                                                                  tAATGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaag        <  2:4226478/93‑1 (MQ=255)
                                                                                    aTGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaagga      <  2:787423/93‑1 (MQ=255)
                                                                                     tGCAGCCTCCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTggaaggaa     >  1:71843/1‑93 (MQ=255)
                                                                                             |                                                                                        
ATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCC  >  NC_000913/577583‑577762

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