Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,419,540 C→T 6.0% P231S (CCG→TCG)  ctpH → methyl‑accepting chemotaxis protein CtpH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,419,5400CT 6.0% 205.3 / 2.6 67P231S (CCG→TCG) ctpHmethyl‑accepting chemotaxis protein CtpH
Reads supporting (aligned to +/- strand):  ref base C (26/37);  new base T (2/2);  total (28/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.54e-01

GCAACTGCGCGCCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGTGACTTCCTTGCCGATACCGTCGGCCGG  >  NC_002947/2419392‑2419686
                                                                                                                                                    |                                                                                                                                                  
gCAACTGCGCGCCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGc                                                                                                                                                    >  8:61646/1‑149 (MQ=255)
gCAACTGCGCGCCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGc                                                                                                                                                    <  1:255485/149‑1 (MQ=255)
gCAACTGCGCGCCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGc                                                                                                                                                    >  2:207676/1‑149 (MQ=255)
gCAACTGCGCGCCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCACACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGc                                                                                                                                                    <  2:336496/149‑1 (MQ=255)
       cgcgcCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGTCGATCCg                                                                                                                                             >  7:228661/1‑149 (MQ=255)
       cgcgcCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGTCGATCCg                                                                                                                                             >  7:222727/1‑149 (MQ=255)
          gcCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCa                                                                                                                                          <  6:259301/149‑1 (MQ=255)
            cGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGt                                                                                                                                        >  7:372771/1‑149 (MQ=255)
              aTGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGc                                                                                                                                            <  4:334001/143‑1 (MQ=255)
              aTGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGc                                                                                                                                            >  3:334001/1‑143 (MQ=255)
                gCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGAt                                                                                                                                    >  5:301170/1‑149 (MQ=255)
                 cTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGTGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGTCGATCCGCCAGTTGATc                                                                                                                                   <  8:222727/149‑1 (MQ=255)
                 cTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGTCGATCCGCCAGTTGATc                                                                                                                                   <  8:228661/149‑1 (MQ=255)
                 cTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATc                                                                                                                                   <  7:366810/149‑1 (MQ=255)
                 cTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATc                                                                                                                                   <  6:301170/149‑1 (MQ=255)
                 cTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATc                                                                                                                                   <  3:156372/149‑1 (MQ=255)
                  tCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGa                                                                                                                                     >  5:213704/1‑146 (MQ=255)
                    gCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGa                                                                                                                                 <  8:56005/148‑1 (MQ=255)
                       aGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTAc                                                                                                                             <  1:209657/149‑1 (MQ=255)
                           tGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGTTGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTcg                                                                                                                         <  6:213704/149‑1 (MQ=255)
                           tGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTcg                                                                                                                         <  8:347108/149‑1 (MQ=255)
                               gccagcATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCa                                                                                                                     <  7:302661/149‑1 (MQ=255)
                                    cATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCATTTCATCGAGTACGTCGCGCAGCTCa                                                                                                                <  7:61646/149‑1 (MQ=255)
                                      tCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGTGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGc                                                                                                              >  2:352315/1‑149 (MQ=255)
                                       cAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCACCCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGc                                                                                                              >  4:334662/1‑148 (MQ=255)
                                             cAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCgcgc                                                                                                                      <  2:365410/134‑1 (MQ=255)
                                             cAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCgcgc                                                                                                                      >  1:365410/1‑134 (MQ=255)
                                              aGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAg                                                                                                      <  3:90030/149‑1 (MQ=255)
                                                  ccGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCttt                                                                                                  >  2:29101/1‑149 (MQ=255)
                                                     aGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTtgct                                                                                               >  5:2213/1‑149 (MQ=255)
                                                     aGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTtgcg                                                                                               >  8:159162/1‑148 (MQ=255)
                                                      gcgcACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGTGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTtgctg                                                                                              <  1:352315/149‑1 (MQ=255)
                                                      gcgcACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTtgctg                                                                                              >  1:42419/1‑149 (MQ=255)
                                                         cACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGccc                                                                                           <  4:150368/149‑1 (MQ=255)
                                                            ggTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCggg                                                                                        <  5:406768/149‑1 (MQ=255)
                                                             gTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGt                                                                                       <  6:134148/149‑1 (MQ=255)
                                                              tATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTc                                                                                      <  3:267646/149‑1 (MQ=255)
                                                                  gTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGTGGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACa                                                                                  <  4:287365/149‑1 (MQ=255)
                                                                  gTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGa                                                                                    >  3:260223/1‑147 (MQ=255)
                                                                  gTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACa                                                                                  <  6:2213/149‑1 (MQ=255)
                                                                  gTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACa                                                                                  <  1:207676/149‑1 (MQ=255)
                                                                      ggCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGcc                                                                               <  8:372771/148‑1 (MQ=255)
                                                                        ccTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTc                                                                                                                 <  8:276196/112‑1 (MQ=255)
                                                                        ccTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTc                                                                                                                 >  7:276196/1‑112 (MQ=255)
                                                                            ctgGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCgg                                                                          >  3:346120/1‑147 (MQ=255)
                                                                            ctgGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGgga                                                                        <  2:42419/149‑1 (MQ=255)
                                                                            ctgGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGgga                                                                        <  6:150254/149‑1 (MQ=255)
                                                                               gTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGggagga                                                                     >  8:225970/1‑149 (MQ=255)
                                                                                    gCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAgg                                                                                                         >  8:247479/1‑108 (MQ=255)
                                                                                    gCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAgg                                                                                                         <  7:247479/108‑1 (MQ=255)
                                                                                      tGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTg                                                              >  7:113469/1‑149 (MQ=255)
                                                                                      tGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTg                                                              >  4:147600/1‑149 (MQ=255)
                                                                                        gCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTg                                                              <  4:346120/147‑1 (MQ=255)
                                                                                         cTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTgggg                                                           >  6:126130/1‑149 (MQ=255)
                                                                                              gCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTg                                                              >  8:239098/1‑141 (MQ=255)
                                                                                              gCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTg                                                              <  7:239098/141‑1 (MQ=255)
                                                                                              gCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCt                                                      >  5:246126/1‑149 (MQ=255)
                                                                                               cGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTg                                                     <  1:29101/149‑1 (MQ=255)
                                                                                                       gtggGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCGTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCgcgg                                               <  3:227952/147‑1 (MQ=255)
                                                                                                          gggCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGggcggc                                          >  7:376267/1‑149 (MQ=255)
                                                                                                          gggCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGggcggc                                          <  7:225970/149‑1 (MQ=255)
                                                                                                                       tgtgGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGt                             <  4:260223/149‑1 (MQ=255)
                                                                                                                         tgGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCgcgggc                                             <  5:176680/131‑1 (MQ=255)
                                                                                                                         tgGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCgcgggc                                             >  6:176680/1‑131 (MQ=255)
                                                                                                                         tgGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGt                             <  6:246126/147‑1 (MQ=255)
                                                                                                                          ggCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGTGAc                          >  3:383608/1‑149 (MQ=255)
                                                                                                                           gCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGTGACt                         <  8:113469/149‑1 (MQ=255)
                                                                                                                                tCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGTGACTTCCtt                    <  3:334662/149‑1 (MQ=255)
                                                                                                                                            tGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGTGACTTCCTTGCCGATACCGTc        >  3:341793/1‑149 (MQ=255)
                                                                                                                                              aTCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGTGACTTCCTTGCCGATACCGTcg       <  5:126130/148‑1 (MQ=255)
                                                                                                                                                  aGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGTGACTTCCTTGCCGATACCGTcggccgg  >  2:33142/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCAACTGCGCGCCGATGCTCGCCAGCGTGCAGCCAGCATCAATGCCAGTGCCGAGCGCACGGTATGGTTGGGCCTGCTGGTGATGCTGGCTTCAGCGGTGCTGGTGGGGCTGCTCAGCCTGTGGCTGGTCAACCGTAGCCTGATCGAGCCGATCCGCCAGTTGATCGAGTACGTCGCGCAGCTCAGCCAAGGCAGCTTTGCTGCCCGGGTCGACAGCCGTCGGGAGGATGAGCTGGGGCGCCTGGCGCGGGCGGCCAACACCCTGCGTGACTTCCTTGCCGATACCGTCGGCCGG  >  NC_002947/2419392‑2419686

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