Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,718,074 T→C S59P (TCC→CCC)  slyB → putative outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,718,0740TC99.4% 514.1 / ‑4.7 160S59P (TCC→CCC) slyBputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (62/97);  total (62/98)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 170 total aligned reads displayed.

TATACCGCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATTCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACT  >  NC_000913/1717984‑1718163
                                                                                          |                                                                                         
tataCCGCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATccc                                                                                         <  2:1801623/93‑1 (MQ=255)
  taCCGCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGGTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCaa                                                                                       <  1:2707920/93‑1 (MQ=255)
    ccGCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACg                                                                                     >  1:971631/1‑93 (MQ=255)
    ccGCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACg                                                                                     >  1:670174/1‑93 (MQ=255)
    ccGCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACg                                                                                     <  1:2075485/93‑1 (MQ=255)
      gCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGtt                                                                                   <  1:4405451/93‑1 (MQ=255)
      gCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGtt                                                                                   >  1:2198552/1‑93 (MQ=255)
        ttCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTAt                                                                                 <  1:89937/93‑1 (MQ=255)
         tccgAAGAGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATc                                                                                >  2:4653875/4‑93 (MQ=255)
         tCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATc                                                                                <  2:4542829/93‑1 (MQ=255)
           tGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCgg                                                                              <  1:194664/93‑1 (MQ=255)
            gAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGt                                                                             <  1:3204706/93‑1 (MQ=255)
               gCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCa                                                                          >  1:1626531/1‑93 (MQ=255)
                 gAAACAAGTACAGAATGTCAGCTATGGCACCGTCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAAt                                                                        <  1:4234535/93‑1 (MQ=255)
                 gAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAAt                                                                        >  2:1477802/1‑93 (MQ=255)
                 gAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGACGTTCCCAACGTTATCGGTGCAAt                                                                        <  1:3928577/93‑1 (MQ=255)
                 gAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGAACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCCAt                                                                        >  1:414326/1‑93 (MQ=255)
                  aaaCAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAAtt                                                                       <  1:3513407/93‑1 (MQ=255)
                  aaaCAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAAtt                                                                       <  2:1266794/93‑1 (MQ=255)
                   aaCAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTg                                                                      >  1:4157093/1‑93 (MQ=255)
                     cAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGc                                                                    <  1:3137182/93‑1 (MQ=255)
                      aaGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCg                                                                   <  2:4668349/93‑1 (MQ=255)
                      aaGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCg                                                                   <  1:657159/93‑1 (MQ=255)
                       aGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCgg                                                                  >  2:3720090/1‑93 (MQ=255)
                        gTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCggg                                                                 >  2:4467938/1‑92 (MQ=255)
                         tACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTg                                                                <  2:551064/93‑1 (MQ=255)
                          aCAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGc                                                               <  2:1928571/93‑1 (MQ=255)
                          aCAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGc                                                               <  1:1514220/93‑1 (MQ=255)
                          aCAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGATACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGc                                                               <  1:1981259/93‑1 (MQ=255)
                          aCAGAATGTCAGCTATGGCACCATCCTTAACGTTCGTCCGGTTCAGATTCAGGGCGGTGATGATCCCAACGTTTTCGGTGCAATTGGCGGTGc                                                               <  2:835283/93‑1 (MQ=255)
                           cAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCt                                                              <  1:1727403/93‑1 (MQ=255)
                            aGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTg                                                             <  2:2422454/93‑1 (MQ=255)
                               aTGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTc                                                          <  2:1295273/93‑1 (MQ=255)
                               aTGTCAGCTATGGCACCATCGTTAAAGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTc                                                          <  2:3318149/93‑1 (MQ=255)
                                 gTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCtt                                                        <  1:3039561/93‑1 (MQ=255)
                                  tCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTtg                                                       <  2:1236701/93‑1 (MQ=255)
                                    aGCTATGGCACCATCGTTAACGTACGCCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTtggt                                                     <  2:1290594/93‑1 (MQ=255)
                                     gCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTtggtg                                                    <  2:1343716/93‑1 (MQ=255)
                                     gCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTtggtg                                                    <  1:2442828/93‑1 (MQ=255)
                                     gCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTtggtg                                                    >  2:2548585/1‑93 (MQ=255)
                                     gCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTtggtg                                                    >  1:3529817/1‑93 (MQ=255)
                                     gCTATGGCACCATCGTTAACGTACGTCCGGGACAGATTCAAGGCGGTGATGATACCAACGGTATCGGTGCGATTGGCGGTGCTGGTTTtggtg                                                    >  2:1458299/1‑93 (MQ=255)
                                       tATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTggtggt                                                  <  1:3823514/93‑1 (MQ=255)
                                       tATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTggtggt                                                  <  1:4392354/93‑1 (MQ=255)
                                        aTGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTtggtgggt                                                 >  2:1754444/1‑91 (MQ=255)
                                        aTGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGtt                                                 >  1:1985486/1‑93 (MQ=255)
                                         tGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGttt                                                <  1:1203958/93‑1 (MQ=255)
                                         tGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGttt                                                >  2:103976/1‑93 (MQ=255)
                                          ggCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTc                                               <  2:1453854/93‑1 (MQ=255)
                                           gCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTcc                                              <  2:4389441/93‑1 (MQ=255)
                                             aCCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTg                                            <  1:1192020/93‑1 (MQ=255)
                                              ccATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTgg                                           >  2:1899800/1‑93 (MQ=255)
                                              ccATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTgg                                           <  2:2268123/93‑1 (MQ=255)
                                              ccATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTgg                                           <  1:1331018/93‑1 (MQ=255)
                                               cATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTggg                                          >  1:1103067/1‑93 (MQ=255)
                                                   gTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAAt                                      <  2:505814/93‑1 (MQ=255)
                                                   gTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAAt                                      >  2:4303483/1‑93 (MQ=255)
                                                      aaCGTACGTCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACt                                   >  1:3592423/1‑93 (MQ=255)
                                                       acgtacgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTg                                  >  1:3121379/1‑93 (MQ=255)
                                                          tacgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtg                               <  1:4001987/93‑1 (MQ=255)
                                                          tacgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtg                               <  1:3958210/93‑1 (MQ=255)
                                                          tacgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtg                               >  1:3891311/1‑93 (MQ=255)
                                                          tacgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtg                               >  1:2784493/1‑93 (MQ=255)
                                                           tcgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtgg                              <  2:1000963/92‑1 (MQ=255)
                                                           acgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtgg                              <  1:2313655/93‑1 (MQ=255)
                                                           acgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtgg                              <  2:2158299/93‑1 (MQ=255)
                                                            cgtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtggt                             <  2:1164902/93‑1 (MQ=255)
                                                             gtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtggtg                            <  2:2625258/93‑1 (MQ=255)
                                                             gtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtggtg                            <  1:1115011/93‑1 (MQ=255)
                                                             gtCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtggtg                            <  2:1115934/93‑1 (MQ=255)
                                                              tCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTtggtgg                           >  1:60728/1‑93 (MQ=255)
                                                              tCCGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGGGGTTTCCTGGGGAATACTGTtggtgg                           >  1:1488571/1‑93 (MQ=255)
                                                                cGGTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCg                         >  1:716445/1‑93 (MQ=255)
                                                                  gTACAGATTCAGGGCGGTGATGATCCCACCGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGa                       <  1:257971/93‑1 (MQ=255)
                                                                  gTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGa                       <  2:1022812/93‑1 (MQ=255)
                                                                  gTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGa                       >  1:616495/1‑93 (MQ=255)
                                                                  gTACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGa                       >  1:4479329/1‑93 (MQ=255)
                                                                   tACAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGTGGaa                      >  1:2814405/1‑93 (MQ=255)
                                                                    aCAGATTCAGGGCGGTGATGTTCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAAc                     <  1:2919026/93‑1 (MQ=255)
                                                                    aCAGATTCAGGGCGGTGATGTTCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAAc                     <  1:2919027/93‑1 (MQ=255)
                                                                    aCAGATTCAGGGCGGTGATGATTCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAAc                     <  2:2013517/93‑1 (MQ=255)
                                                                    aCAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAAc                     <  1:2109501/93‑1 (MQ=255)
                                                                    aCAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAAc                     <  1:4161816/93‑1 (MQ=255)
                                                                     cAGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAAcc                    <  2:3695223/93‑1 (MQ=255)
                                                                      aGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCg                   >  1:3244744/1‑93 (MQ=255)
                                                                      aGATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCg                   <  1:1147556/93‑1 (MQ=255)
                                                                       gATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCgg                  <  2:3431550/93‑1 (MQ=255)
                                                                       gATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCgg                  <  2:387588/93‑1 (MQ=255)
                                                                       gATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCgg                  >  2:3252637/1‑93 (MQ=255)
                                                                       gATTCAGGGCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCgg                  <  2:2819657/93‑1 (MQ=255)
                                                                             gggCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTtc            >  2:3347926/1‑93 (MQ=255)
                                                                             gggCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTtc            <  2:2357160/93‑1 (MQ=255)
                                                                               gCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTtctc          <  1:2544164/93‑1 (MQ=255)
                                                                               gCGGTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTtctc          <  2:693036/93‑1 (MQ=255)
                                                                                 ggTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTg        <  2:3888622/93‑1 (MQ=255)
                                                                                 ggTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTg        <  2:4080319/93‑1 (MQ=255)
                                                                                  gTGATGATCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTgg       <  1:3842021/93‑1 (MQ=255)
                                                                                    gatgatCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCt     >  2:4574326/1‑93 (MQ=255)
                                                                                     atgatCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTa    >  2:3060588/1‑93 (MQ=255)
                                                                                     atgatCCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTa    >  2:2722034/1‑93 (MQ=255)
                                                                                          |                                                                                         
TATACCGCTTCTGAAGCGAAACAAGTACAGAATGTCAGCTATGGCACCATCGTTAACGTACGTCCGGTACAGATTCAGGGCGGTGATGATTCCAACGTTATCGGTGCAATTGGCGGTGCTGTTCTTGGTGGTTTCCTGGGGAATACTGTTGGTGGCGGAACCGGGCGTTCTCTGGCTACT  >  NC_000913/1717984‑1718163

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: