Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,318,095 C→G 38.6% L594F (TTG→TTC 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,0950CG38.6% 75.6 / 114.9 88L594F (TTG→TTCrcsChybrid sensory kinase in two‑component regulatory system with RcsB and YojN
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base G (18/16);  total (43/45)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-02

GACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCAC  >  NC_000913/2317948‑2318240
                                                                                                                                                   |                                                                                                                                                 
gACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCa                                                                                                                                                  >  3:742302/1‑149 (MQ=255)
   cTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACa                                                                                                                                               <  4:743976/149‑1 (MQ=255)
    tGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACa                                                                                                                                               >  1:138295/1‑148 (MQ=255)
    tGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAg                                                                                                                                              >  4:228566/1‑149 (MQ=255)
    tGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAg                                                                                                                                              >  4:219900/1‑149 (MQ=255)
     ggAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAgg                                                                                                                                             >  3:119502/1‑149 (MQ=255)
     ggAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAgg                                                                                                                                             >  3:348578/1‑149 (MQ=255)
      gaagaaGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGt                                                                                                                                            <  4:394460/149‑1 (MQ=255)
      gaagaaGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGt                                                                                                                                            <  4:32549/149‑1 (MQ=255)
        agaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGtt                                                                                                                                           <  1:897954/148‑1 (MQ=255)
        agaagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGtt                                                                                                                                           <  2:140432/148‑1 (MQ=255)
          aagGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGga                                                                                                                                        <  3:286836/149‑1 (MQ=255)
          aagGGAGCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTgg                                                                                                                                         >  2:137725/1‑148 (MQ=255)
              gATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGgag                                                                                                                                       >  4:414058/1‑146 (MQ=255)
              gATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGgag                                                                                                                                       <  3:414058/146‑1 (MQ=255)
              gATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGAt                                                                                                                                     >  2:141756/1‑148 (MQ=255)
               aTCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATg                                                                                                                                    >  2:652567/1‑148 (MQ=255)
                    aCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACa                                                                                                                                               >  2:393652/1‑132 (MQ=255)
                    aCAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACa                                                                                                                                               <  1:393652/132‑1 (MQ=255)
                     cAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACctgc                                                                                                                              >  1:84718/1‑148 (MQ=255)
                          gcACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTaa                                                                                                                         <  2:84718/148‑1 (MQ=255)
                             ccacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACg                                                                                                                      >  2:509774/1‑148 (MQ=255)
                             ccacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAAGAGGTTGGAGATGACCTGCTGTAAACGc                                                                                                                     >  3:657931/1‑149 (MQ=255)
                              cacTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCa                                                                                                                    >  3:439057/1‑149 (MQ=255)
                               acTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCa                                                                                                                    <  1:141756/148‑1 (MQ=255)
                                cttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCAt                                                                                                                   >  1:905538/1‑148 (MQ=255)
                                cttcttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCAt                                                                                                                   >  1:367182/1‑148 (MQ=255)
                                   cttTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCgg                                                                                                                <  1:532971/148‑1 (MQ=255)
                                     ttCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTc                                                                                                             <  4:802638/149‑1 (MQ=255)
                                              aTCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTa                                                                                                     <  2:447395/148‑1 (MQ=255)
                                               tCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTaa                                                                                                    <  2:166710/148‑1 (MQ=255)
                                                ccccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAg                                                                                                   <  2:367182/148‑1 (MQ=255)
                                                ccccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATCCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATCACCTGCTGTAAACGCATCGGGTCGCCATTTAAg                                                                                                   <  1:730727/148‑1 (MQ=255)
                                                 cccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGc                                                                                                 >  4:666679/1‑149 (MQ=255)
                                                  ccACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGcc                                                                                                >  4:938594/1‑149 (MQ=255)
                                                             tCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATc                                                                                     <  4:348578/149‑1 (MQ=255)
                                                               gcgAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCg                                                                                   >  3:174922/1‑149 (MQ=255)
                                                                  aaCACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATCTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATGCGGtt                                                                                <  4:815085/149‑1 (MQ=255)
                                                                   acacGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACAGCCGGtac                                                                               <  4:115502/149‑3 (MQ=255)
                                                                     acGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCaa                                                                             <  4:917624/149‑1 (MQ=255)
                                                                     acGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCaa                                                                             <  4:294259/149‑1 (MQ=255)
                                                                     acGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCa                                                                              <  1:132188/148‑1 (MQ=255)
                                                                     acGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCa                                                                              <  1:835170/148‑1 (MQ=255)
                                                                          tAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAg                                                                        <  3:91613/149‑1 (MQ=255)
                                                                             agagaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAg                                                                     >  3:318323/1‑149 (MQ=255)
                                                                               agaTAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTa                                                                   >  3:319065/1‑149 (MQ=255)
                                                                                    aTCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAg                                                               >  2:10234/1‑148 (MQ=255)
                                                                                    aTCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGc                                                              <  3:938594/149‑1 (MQ=255)
                                                                                     tCGCTATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATaa                                                                          >  1:808679/1‑136 (MQ=255)
                                                                                     tCGCTATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATaa                                                                          <  2:808679/136‑1 (MQ=255)
                                                                                      cGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGTCACTGGCACATCCGGTTCAATAAAGCATGACAAGCCt                                                            <  4:684637/149‑1 (MQ=255)
                                                                                        ccATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTa                                                           >  2:418929/1‑148 (MQ=255)
                                                                                         cATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTaa                                                          <  2:817055/148‑1 (MQ=255)
                                                                                              gcgcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCt                                                            <  3:259125/141‑1 (MQ=255)
                                                                                              gcgcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCt                                                            >  4:259125/1‑141 (MQ=255)
                                                                                               cgcgAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTCCAAGCCTAACTGCtt                                                    <  2:200981/148‑1 (MQ=255)
                                                                                                        gCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTaccac                                           >  1:273258/1‑148 (MQ=255)
                                                                                                         cAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACaa                                                                >  3:306916/1‑126 (MQ=255)
                                                                                                         cAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAAGAAAGCAGTACaa                                                                <  4:306916/126‑1 (MQ=255)
                                                                                                         cAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCa                                                                              >  2:686501/1‑112 (MQ=255)
                                                                                                         cAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCa                                                                              <  1:686501/112‑1 (MQ=255)
                                                                                                          aaaaCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTaa                                                                                                    >  3:623247/1‑89 (MQ=255)
                                                                                                          aaaaCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTaa                                                                                                    <  4:623247/89‑1 (MQ=255)
                                                                                                           aaaCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAg                                        >  1:787335/1‑148 (MQ=255)
                                                                                                           aaaCTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAg                                        >  1:711725/1‑148 (MQ=255)
                                                                                                              cTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCgg                                     <  2:336429/148‑1 (MQ=255)
                                                                                                                ataCAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTa                                   >  1:672006/1‑148 (MQ=255)
                                                                                                                  aCAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAAt                                <  4:319065/149‑1 (MQ=255)
                                                                                                                   cAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCg                                                                                   >  2:402327/1‑97 (MQ=255)
                                                                                                                   cAGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCg                                                                                   <  1:402327/97‑1 (MQ=255)
                                                                                                                    aGCCGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAg                              >  4:198349/1‑149 (MQ=255)
                                                                                                                     gCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGt                             <  4:715360/149‑1 (MQ=255)
                                                                                                                      ccGGTATCGGTGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGt                             >  2:131497/1‑148 (MQ=255)
                                                                                                                      ccGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGtt                            <  4:174922/149‑1 (MQ=255)
                                                                                                                      ccGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGt                             >  2:634560/1‑148 (MQ=255)
                                                                                                                              ggTGAATTTTATGGCGTTACTCAGCAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAACGCCCCTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTg                     <  2:873113/148‑1 (MQ=255)
                                                                                                                               gTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGa                    <  2:711725/148‑1 (MQ=255)
                                                                                                                                tGAATTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAt                   >  1:585838/1‑148 (MQ=255)
                                                                                                                                 gAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAtgt                 <  4:318323/149‑1 (MQ=255)
                                                                                                                                   aTTTTATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGAtgtg                <  1:502894/148‑1 (MQ=255)
                                                                                                                                    ttttATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGt              >  3:739346/1‑149 (MQ=255)
                                                                                                                                       tATGGCGTTACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtca           <  3:198349/149‑1 (MQ=255)
                                                                                                                                           gCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtcatca        >  1:594111/1‑148 (MQ=255)
                                                                                                                                           gCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtcat          <  2:206673/146‑1 (MQ=255)
                                                                                                                                           gCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTtcat          >  1:206673/1‑146 (MQ=255)
                                                                                                                                              ttACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAg                              >  1:394930/1‑123 (MQ=255)
                                                                                                                                              ttACTGAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAg                              <  2:394930/123‑1 (MQ=255)
                                                                                                                                                aCTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCAc  >  4:395948/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GACCTGGAAGAAGGGATCAAACAAGCGCACCACTTCTTTCGCCGGTATCCCCACGCCGGTATCGCGAACACGGATAGAGAGATAATCGCCATCCGCGCGAACATGCAAAACTATACAGCCGGTATCGGTGAATTTTATGGCGTTACTCAACAGGTTGGAGATGACCTGCTGTAAACGCATCGGGTCGCCATTTAAGGCCACTGGCACATCCGGTTCAATAAAGCAGTACAAGCCTAACTGCTTGCGTACCACCAGCGGTAAATAGTTGGCGGTGATGTGGTTCATCACTTCAC  >  NC_000913/2317948‑2318240

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