Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,318,639 A→C F413C (TTT→TGT)  rcsC ← hybrid sensory kinase in two‑component regulatory system with RcsB and YojN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,318,6390AC99.3% 492.6 / ‑4.2 139F413C (TTT→TGT) rcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (76/62);  total (77/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

Only 100 of 149 total aligned reads displayed.

TCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATAT  >  NC_000913/2318543‑2318738
                                                                                                |                                                                                                   
tCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAAt                                                                                                  >  2:551296/1‑100 (MQ=255)
tCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAAt                                                                                                  <  2:1014716/100‑1 (MQ=255)
  caccaGCACACAAATGGCCACGTTTTCCATACGATAACGCGAATCGACGCTGCAGACTCACAGACTAGTATTGTTCCTGCTCAGGACATCCACACAAttt                                                                                                >  2:1806854/1‑99 (MQ=255)
  caccaGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTg                                                                                                >  1:882743/1‑100 (MQ=255)
   accaGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGa                                                                                               >  2:1090638/1‑100 (MQ=255)
     caGCACACAAATGGCCACGTTTTCATTTCTATAGCGCGAAGGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGTcc                                                                                             <  2:692379/100‑1 (MQ=255)
     caGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGAcc                                                                                             <  2:2980995/100‑1 (MQ=255)
     caGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGAcc                                                                                             <  2:1671559/100‑1 (MQ=255)
      aGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACct                                                                                            >  1:1026324/1‑100 (MQ=255)
      aGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACct                                                                                            <  2:2093799/100‑1 (MQ=255)
       gCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCGGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctg                                                                                           <  2:955141/100‑1 (MQ=255)
        cacacaAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctgc                                                                                          <  2:650812/100‑1 (MQ=255)
          cacaAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctgctg                                                                                        >  2:3174334/1‑100 (MQ=255)
          cacaAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctgctg                                                                                        >  1:2306078/1‑100 (MQ=255)
          cacaAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctgctg                                                                                        >  2:1272201/1‑100 (MQ=255)
          cacaAATGGCCACGTTTTCATTACGATAGCGAGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACctgctg                                                                                        >  1:3271000/1‑100 (MQ=255)
            gaactGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGAGGGGTATTGTCGCTGGTCACGACATGACCACAATTTGCCTGCTGcc                                                                                      >  2:347614/5‑100 (MQ=255)
                tGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCcaca                                                                                  <  1:1206066/100‑1 (MQ=255)
                tGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTCACCTGCTGCCcaca                                                                                  >  1:2999533/1‑100 (MQ=255)
                 ggCCACGTTTTCATTACGATAGCGCGACTGGACGAAGCTGATTTGCAGCTTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCcacga                                                                                 >  2:3979484/1‑98 (MQ=255)
                 ggCCACGTTTTCATTACGATAGCGCGAAGGGACGAAGCAGATTCGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAg                                                                                 <  2:3851363/100‑1 (MQ=255)
                  tccACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGa                                                                                <  2:3177475/99‑1 (MQ=255)
                  gCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGa                                                                                >  1:4064904/1‑100 (MQ=255)
                   ccACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGAt                                                                               >  2:2975565/1‑100 (MQ=255)
                   ccACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGAt                                                                               >  2:632036/1‑100 (MQ=255)
                   ccACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGAt                                                                               >  1:2196311/1‑100 (MQ=255)
                    cACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATa                                                                              >  2:2287103/1‑100 (MQ=255)
                      cGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAAt                                                                            <  2:253361/100‑1 (MQ=255)
                        ttttCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAAttt                                                                          <  2:2115204/100‑1 (MQ=255)
                            cATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGt                                                                      >  2:248949/1‑100 (MQ=255)
                             ataaCGATGGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACCATTGACCTGCTGCCCACAGATAAGTTGCGtn                                                                     <  2:3073907/97‑2 (MQ=255)
                             aTTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTc                                                                     >  1:2458100/1‑100 (MQ=255)
                             aTTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTc                                                                     >  1:258216/1‑100 (MQ=255)
                             aTTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTc                                                                     >  2:3540659/1‑100 (MQ=255)
                             aTTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTc                                                                     >  2:2944938/1‑100 (MQ=255)
                               tACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAg                                                                   >  2:2854906/1‑100 (MQ=255)
                                 cGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTc                                                                 >  2:3121883/1‑100 (MQ=255)
                                  gATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCg                                                                <  2:773926/100‑1 (MQ=255)
                                      gcgcGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTgg                                                            >  1:2955785/1‑100 (MQ=255)
                                      gcgcGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTCGCGTCAGTCGTTgg                                                            >  1:3163877/1‑100 (MQ=255)
                                       cgcgAATGGACGAAGCTGCTTTCCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCACTCGTTGGc                                                           >  2:787594/1‑100 (MQ=255)
                                       cgcgAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGc                                                           <  2:1984877/100‑1 (MQ=255)
                                       cgcgAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGc                                                           <  1:3914313/100‑1 (MQ=255)
                                       cgcgAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGc                                                           >  2:2689803/1‑100 (MQ=255)
                                         cgAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCgg                                                         >  2:1160471/1‑100 (MQ=255)
                                          gAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGt                                                        >  1:1561431/1‑100 (MQ=255)
                                          gAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGt                                                        >  1:1984743/1‑100 (MQ=255)
                                           aaTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCGCCAGGTAATTTGCGTCAGTCGTTGGCGGtg                                                       >  2:936975/1‑99 (MQ=255)
                                           aaTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTc                                                       >  1:1124954/1‑100 (MQ=255)
                                            tctgACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTcc                                                      <  1:2356340/97‑1 (MQ=255)
                                            aTGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTcc                                                      <  2:482855/100‑1 (MQ=255)
                                             tGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCt                                                     <  2:1362022/100‑1 (MQ=255)
                                             tGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCt                                                     <  1:2583608/100‑1 (MQ=255)
                                               gACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTccaga                                                   >  2:2320108/1‑97 (MQ=255)
                                                 cGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGTCCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCCCATg                                                 <  2:844917/100‑1 (MQ=255)
                                                 cGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATg                                                 >  2:360939/1‑100 (MQ=255)
                                                      cTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTaa                                            >  2:3900322/1‑100 (MQ=255)
                                                       tGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAg                                           >  1:133651/1‑100 (MQ=255)
                                                        gATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACGTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGc                                          >  1:3583367/1‑100 (MQ=255)
                                                         aTTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCa                                         <  1:144098/100‑1 (MQ=255)
                                                         aTTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCa                                         <  1:3421927/100‑1 (MQ=255)
                                                          tttGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCat                                        >  2:1990033/1‑100 (MQ=255)
                                                            tGCAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCatat                                      >  2:3713821/1‑100 (MQ=255)
                                                              cAGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCTTACGCatgagg                                    >  1:2396502/1‑96 (MQ=255)
                                                              cAGATTGGTATTGATGCTGGTCAGGACCTCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCGTGCGTAAGCATATTg                                    >  2:3744265/1‑100 (MQ=255)
                                                               aGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTga                                   <  2:3252078/100‑1 (MQ=255)
                                                               aGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTga                                   <  2:617907/100‑1 (MQ=255)
                                                               aGATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTga                                   <  2:3738814/100‑1 (MQ=255)
                                                                gATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                  <  1:1767279/100‑1 (MQ=255)
                                                                gATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                  >  1:3225911/1‑100 (MQ=255)
                                                                gATTGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTgag                                  >  1:1612457/1‑100 (MQ=255)
                                                                   tGGTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATa                               <  2:1597680/100‑1 (MQ=255)
                                                                    ggTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAg                              >  1:1609752/1‑100 (MQ=255)
                                                                    ggTATTGTTGCTGGTCAGGACAGCAACACAATTGACCTGCTGCCCACATATAATTTGCGTCAGTCGTGGGCGGTCCTCATGCGTAAGCATATTGAGATAg                              <  2:295620/100‑1 (MQ=255)
                                                                     gTATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAgg                             <  1:2493688/100‑1 (MQ=255)
                                                                      tATTGTTGCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGt                            <  1:406938/100‑1 (MQ=255)
                                                                        ttgttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTAt                          >  2:434066/1‑100 (MQ=255)
                                                                        ttgttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTAt                          >  2:3841395/1‑100 (MQ=255)
                                                                           ttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcg                       <  2:2853750/100‑1 (MQ=255)
                                                                           ttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcg                       <  2:1002114/100‑1 (MQ=255)
                                                                           ttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcg                       <  1:442466/100‑1 (MQ=255)
                                                                           ttgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcg                       <  1:3590346/100‑1 (MQ=255)
                                                                            tgCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATgcgc                      >  2:301464/1‑100 (MQ=255)
                                                                             gCTGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGcc                     >  1:1556545/1‑100 (MQ=255)
                                                                               tGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAg                   <  1:455963/100‑1 (MQ=255)
                                                                               tGGTCAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAg                   <  2:1089988/100‑1 (MQ=255)
                                                                                   cAGGACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCg               >  2:1718128/1‑100 (MQ=255)
                                                                                     ggACATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGtt             <  1:3931940/100‑1 (MQ=255)
                                                                                       aCATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTAc           <  2:3014859/100‑1 (MQ=255)
                                                                                       aCATCAACACAATTGACCGGCTGCCCACAGATAATTTTCGTCAGGCGTGGGCGTTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTAc           <  1:4079517/100‑1 (MQ=255)
                                                                                        cATCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGAGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACt          >  1:1938171/1‑100 (MQ=255)
                                                                                         aTCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACtt         >  2:3709487/1‑100 (MQ=255)
                                                                                          tCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTa        <  2:1782519/100‑1 (MQ=255)
                                                                                          tCAACACAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTa        <  1:1069175/100‑1 (MQ=255)
                                                                                           caacacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCGTATTGAGATAGGTATGCGCCAGTTCGTTACTTaa       >  1:2015652/1‑100 (MQ=255)
                                                                                           caacacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTaa       >  1:829325/1‑100 (MQ=255)
                                                                                           caacacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCAAGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTaa       >  2:2746498/1‑100 (MQ=255)
                                                                                            aacacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTaaa      <  1:151591/100‑1 (MQ=255)
                                                                                             acacAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTaaaa     >  2:1175184/1‑100 (MQ=255)
                                                                                                aaaTTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAatat  >  1:915684/1‑100 (MQ=255)
                                                                                                |                                                                                                   
TCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCCCACAGATAATTTGCGTCAGTCGTTGGCGGTCCTCATGCGTAAGCATATTGAGATAGGTATGCGCCAGTTCGTTACTTAAAATAT  >  NC_000913/2318543‑2318738

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