Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,318,506 C→T 5.2% Q457Q (CAG→CAA 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,5060CT 5.2% 251.0 / 3.0 77Q457Q (CAG→CAArcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base C (39/34);  new base T (2/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

GCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCC  >  NC_000913/2318358‑2318654
                                                                                                                                                    |                                                                                                                                                    
gCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTcc                                                                                                                                                      >  2:622325/1‑149 (MQ=255)
 ccGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCt                                                                                                                                                     <  3:141246/149‑1 (MQ=255)
     tCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCaa                                                                                                                                                 >  3:208254/1‑149 (MQ=255)
     tCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCaa                                                                                                                                                 <  4:230335/149‑1 (MQ=255)
      cGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAAc                                                                                                                                                <  3:526945/149‑1 (MQ=255)
       gACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACg                                                                                                                                               >  4:87689/1‑149 (MQ=255)
        aCGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTACAACGa                                                                                                                                              >  4:543307/1‑149 (MQ=255)
           ccTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGActc                                                                                                                                           >  2:658707/1‑149 (MQ=255)
            cTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGActct                                                                                                                                          >  3:10267/1‑149 (MQ=255)
            cTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGActct                                                                                                                                          >  1:381019/1‑149 (MQ=255)
             tttCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCtt                                                                                                                                         >  3:444406/1‑149 (MQ=255)
             tttCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCtt                                                                                                                                         >  3:301103/1‑149 (MQ=255)
                cGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCa                                                                                                                                      <  1:679052/149‑1 (MQ=255)
                 ggTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCAt                                                                                                                                     >  2:20428/1‑149 (MQ=255)
                   tAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCt                                                                                                                                   >  4:208230/1‑149 (MQ=255)
                    aaCTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCtt                                                                                                                                  <  1:45998/149‑1 (MQ=255)
                    aaCTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCtt                                                                                                                                  >  2:45998/1‑149 (MQ=255)
                       tctTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCAc                                                                                                                               >  2:209537/1‑149 (MQ=255)
                       tctTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCAc                                                                                                                               <  2:209073/149‑1 (MQ=255)
                       tctTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCAc                                                                                                                               >  2:188356/1‑149 (MQ=255)
                              tttGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCgaa                                                                                                                        <  3:29731/149‑1 (MQ=255)
                                   aaCAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAAc                                                                                                                   <  3:281234/149‑1 (MQ=255)
                                    aCAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACa                                                                                                                  >  1:432321/1‑149 (MQ=255)
                                    aCAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACa                                                                                                                  >  1:409789/1‑149 (MQ=255)
                                     cAGATCCAGGTTACGGATAATGCTATTCAGCGGCGTTCGCAGCTCATGACTGCCGGTGGCAAGGAACAGCGATTCTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCCCCTGCAACGACTTTTCCATCTTCACGCGCGAAGAAACAt                                                                                                                 <  2:515455/149‑1 (MQ=255)
                                       gATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCGCGCGTgaagatggaaga                                                                                                               <  2:104063/149‑8 (MQ=255)
                                           cAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATccacca                                                                                                           <  3:172165/149‑1 (MQ=255)
                                               ttACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGcac                                                                                                       >  1:148943/1‑149 (MQ=255)
                                                tACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGcaca                                                                                                      >  4:116877/1‑149 (MQ=255)
                                                tACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGcaca                                                                                                      >  1:331750/1‑149 (MQ=255)
                                                     aTAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAAt                                                                                                 <  2:239674/149‑1 (MQ=255)
                                                      tAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATg                                                                                                <  2:96281/149‑1 (MQ=255)
                                                        aTGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGc                                                                                              <  1:20428/149‑1 (MQ=255)
                                                        aTGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGc                                                                                              <  4:208254/149‑1 (MQ=255)
                                                             aTACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGAGAAGAAACATCCACCAGCACACAAATGGCCACGt                                                                                         >  1:487387/1‑149 (MQ=255)
                                                                 aGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTc                                                                                     >  2:65462/1‑149 (MQ=255)
                                                                  gcggcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTTGCTCGCCTGCTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCGACGTTTTCa                                                                                    <  4:80984/149‑1 (MQ=255)
                                                                  gcggcgTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCa                                                                                    <  1:277114/149‑1 (MQ=255)
                                                                        ttCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGa                                                                              <  1:553268/149‑1 (MQ=255)
                                                                          cGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATa                                                                            <  2:381019/149‑1 (MQ=255)
                                                                           gCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAg                                                                           >  3:639837/1‑149 (MQ=255)
                                                                            cAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAgc                                                                          <  4:294592/149‑1 (MQ=255)
                                                                              gCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCTTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAgcgc                                                                        >  3:410120/1‑149 (MQ=255)
                                                                              gCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCTTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAgcgc                                                                        >  1:57866/1‑149 (MQ=255)
                                                                                 cATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGaa                                                                     <  4:192496/149‑1 (MQ=255)
                                                                                 cATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGaa                                                                     <  4:301103/149‑1 (MQ=255)
                                                                                     actgacGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGa                                                                 <  2:148943/149‑1 (MQ=255)
                                                                                     actgacGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGa                                                                 <  1:285044/149‑1 (MQ=255)
                                                                                       tgacGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACg                                                               >  2:548403/1‑149 (MQ=255)
                                                                                                 aaGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGAttt                                                     >  1:703553/1‑149 (MQ=255)
                                                                                                  aGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTg                                                    >  4:353103/1‑149 (MQ=255)
                                                                                                   ggAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGc                                                   <  3:533629/149‑1 (MQ=255)
                                                                                                        aTCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCTTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGAtt                                              <  4:410120/149‑1 (MQ=255)
                                                                                                        aTCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCTTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGAtt                                              <  2:57866/149‑1 (MQ=255)
                                                                                                        aTCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGAtt                                              >  2:698385/1‑149 (MQ=255)
                                                                                                        aTCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGAtt                                              >  2:261680/1‑149 (MQ=255)
                                                                                                        aTCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGAtt                                              >  4:107030/1‑149 (MQ=255)
                                                                                                         tCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTg                                             >  3:219599/1‑149 (MQ=255)
                                                                                                         tCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTg                                             >  1:480252/1‑149 (MQ=255)
                                                                                                            aTTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTa                                          <  4:639837/149‑1 (MQ=255)
                                                                                                            aTTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTa                                          <  3:208230/149‑1 (MQ=255)
                                                                                                             ttttGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTGATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTAt                                         >  1:228176/1‑149 (MQ=255)
                                                                                                              tttGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTAtt                                        <  1:661938/149‑1 (MQ=255)
                                                                                                                    tGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCt                                  <  1:37294/149‑1 (MQ=255)
                                                                                                                            cTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTg                                 >  2:464307/1‑142 (MQ=255)
                                                                                                                            cTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTg                                 <  1:464307/142‑1 (MQ=255)
                                                                                                                            cTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGa                          >  1:295217/1‑149 (MQ=255)
                                                                                                                              gTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACa                        >  1:204181/1‑149 (MQ=255)
                                                                                                                                 ccGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATca                     <  4:444406/149‑1 (MQ=255)
                                                                                                                                   gctgctTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATcaac                   >  4:112612/1‑149 (MQ=255)
                                                                                                                                        ttGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAAt              >  2:119629/1‑149 (MQ=255)
                                                                                                                                         tgtgCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAAtt             <  3:518664/149‑1 (MQ=255)
                                                                                                                                           tgCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGa           <  2:331750/149‑1 (MQ=255)
                                                                                                                                             ccATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGAcc         >  1:350939/1‑149 (MQ=255)
                                                                                                                                             ccATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGCCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGAcc         >  3:347787/1‑149 (MQ=255)
                                                                                                                                                tctcCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACaaaa               <  1:3866/140‑1 (MQ=255)
                                                                                                                                                tctcCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACaaaa               >  2:3866/1‑140 (MQ=255)
                                                                                                                                                tctcCGGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGTACGAAGCTGATTTGCAGCTTCGTATTGTTTCTGGTCAGGACATCAACAAAATTGACctgc      >  2:334825/1‑149 (MQ=255)
                                                                                                                                                    cTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGcc  >  4:114867/1‑149 (MQ=255)
                                                                                                                                                    cTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGcc  <  2:204181/149‑1 (MQ=255)
                                                                                                                                                    cTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGcc  >  3:376033/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCGATCGACGCCTTTCGGTAACTCTTTGGTTTGCAACAGATCCAGGTTACCGATAATGCCATACAGCGGCGTTCGCAGCTCATGACTGACGGTGGCAAGGAACATCGATTTTGACTGGCTCGCCTGTTCCGCTGCTTGTGCCATCTCCTGCAACGACTCTTCCATCTTCACGCGCGAAGAAACATCCACCAGCACACAAATGGCCACGTTTTCATTACGATAGCGCGAATGGACGAAGCTGATTTGCAGATTGGTATTGTTGCTGGTCAGGACATCAACAAAATTGACCTGCTGCC  >  NC_000913/2318358‑2318654

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