Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 681,972 C→A 5.3% S434I (AGC→ATC)  PP_0584 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947681,9720CA 5.3% 276.6 / 3.2 85S434I (AGC→ATC) PP_0584methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (44/36);  new base A (3/2);  total (47/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTGCCGAGGATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTC  >  NC_002947/681824‑682115
                                                                                                                                                    |                                                                                                                                                
ggTGCCGAGGATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGc                                                                                                                                                  >  6:167950/1‑149 (MQ=255)
   gCCGAGGATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTAt                                                                                                                                               >  5:74041/1‑149 (MQ=255)
      gAGGATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCgg                                                                                                                                            >  2:29497/1‑149 (MQ=255)
         gATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCg                                                                                                                                        <  1:286763/149‑1 (MQ=255)
         gATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGATAGCGGcc                                                                                                                                         <  4:345608/148‑1 (MQ=255)
                   aTGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGcgccgc                                                                                                                              >  3:140050/1‑149 (MQ=255)
                            cTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGc                                                                                                                      <  7:23974/148‑1 (MQ=255)
                            cTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCa                                                                                                                     <  3:26819/149‑1 (MQ=255)
                                cACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGAc                                                                                                                 >  7:71773/1‑149 (MQ=255)
                                 aCGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGAcc                                                                                                                >  5:26351/1‑149 (MQ=255)
                                 aCGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGATAGCGGccagca                                                                                                                                     >  8:49661/1‑124 (MQ=255)
                                 aCGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGATAGCGGccagca                                                                                                                                     <  7:49661/128‑5 (MQ=255)
                                    cTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGTCCGACTGCGCCGCCTGGCTGCAGGACCGt                                                                                                              <  1:29497/149‑1 (MQ=255)
                                       cAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGc                                                                                                            <  3:10077/147‑1 (MQ=255)
                                       cAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGc                                                                                                            <  2:353163/147‑1 (MQ=255)
                                       cAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTACTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCa                                                                                                          >  5:165659/1‑149 (MQ=255)
                                        aGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCAc                                                                                                         >  7:257601/1‑149 (MQ=255)
                                          cTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACt                                                                                                        >  3:232532/1‑148 (MQ=255)
                                             ttGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCa                                                                                                     <  3:287606/148‑1 (MQ=255)
                                                gATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTgg                                                                                                 >  1:294745/1‑149 (MQ=255)
                                                 aTCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGc                                                                                                <  5:100791/149‑1 (MQ=255)
                                                  tCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCg                                                                                               >  1:69344/1‑149 (MQ=255)
                                                      cTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTgg                                                                                           >  3:167992/1‑149 (MQ=255)
                                                        gggCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTCGCTGCAGGACCGTGCCACTTCATTGGCGGTGGc                                                                                          >  1:224624/1‑148 (MQ=255)
                                                         ggCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGa                                                                                        >  2:340379/1‑149 (MQ=255)
                                                          gCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGAc                                                                                       <  6:74041/149‑1 (MQ=255)
                                                          gCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGAc                                                                                       <  6:80877/149‑1 (MQ=255)
                                                           ccGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGAc                                                                                       >  6:99505/1‑148 (MQ=255)
                                                             gATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCAt                                                                                    >  3:66791/1‑149 (MQ=255)
                                                                 gCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGt                                                                                             <  3:236369/136‑1 (MQ=255)
                                                                 gCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGt                                                                                             >  4:236369/1‑136 (MQ=255)
                                                                  cTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGt                                                                               >  3:140401/1‑149 (MQ=255)
                                                                  cTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGt                                                                               <  5:167950/149‑1 (MQ=255)
                                                                  cTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGt                                                                               <  2:230819/149‑1 (MQ=255)
                                                                  cTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGt                                                                               <  5:99505/149‑1 (MQ=255)
                                                                   tCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCgg                                                                                              >  5:285821/1‑133 (MQ=255)
                                                                   tCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCgg                                                                                              <  6:285821/133‑1 (MQ=255)
                                                                   tCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTg                                                                              >  5:110447/1‑149 (MQ=255)
                                                                        ttCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATg                                                                         >  7:252220/1‑149 (MQ=255)
                                                                            tGGCTAACGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGt                                                                     <  6:193354/149‑1 (MQ=255)
                                                                                tAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGa                                                                  >  7:138468/1‑148 (MQ=255)
                                                                                tAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGa                                                                  >  2:279569/1‑148 (MQ=255)
                                                                                   ggCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCACCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAAc                                                               >  3:326834/1‑148 (MQ=255)
                                                                                    gCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCa                                                             >  4:15356/1‑149 (MQ=255)
                                                                                      gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGc                                                                                          >  1:57334/1‑118 (MQ=255)
                                                                                      gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGc                                                                                          <  2:57334/118‑1 (MQ=255)
                                                                                      gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATg                                                           >  8:76798/1‑149 (MQ=255)
                                                                                      gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATg                                                           >  6:316544/1‑149 (MQ=255)
                                                                                      gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATg                                                           >  3:266516/1‑149 (MQ=255)
                                                                                           aCGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCAc                                                      <  7:76798/149‑1 (MQ=255)
                                                                                                 tgtACCGCCGCATCGATCTGCTGCTGCCCTTGGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCtt                                                 <  2:224624/148‑1 (MQ=255)
                                                                                                   tACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTcgc                                              <  4:140401/149‑1 (MQ=255)
                                                                                                    aCCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTcgcg                                             >  8:317/1‑149 (MQ=255)
                                                                                                     ccgccgCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGTAACCATGTCCACGGCTTCgcgc                                            >  3:99259/1‑149 (MQ=255)
                                                                                                     ccgccgCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCgcgc                                            <  4:140050/149‑1 (MQ=255)
                                                                                                     ccgccgCATCGATCTGCTGCTGCCCTTCGCGGGCCTGATGCTGGCCGCTATCGGCCGACTCCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAACCCATGTCCACGGCttacccg                                            >  5:145345/1‑144 (MQ=255)
                                                                                                     ccgccgCATCGATCTGCTGCTGCCCTTCGCGGGCCTACTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCgcgc                                            <  6:165659/149‑1 (MQ=255)
                                                                                                       gccgcATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTg                                          <  1:279569/149‑1 (MQ=255)
                                                                                                        ccgcCTCGATGTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTgg                                         <  8:139907/149‑1 (MQ=255)
                                                                                                        ccgcATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGt                                                                                                              <  1:353441/80‑1 (MQ=255)
                                                                                                        ccgcATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGt                                                                                                              >  2:353441/1‑80 (MQ=255)
                                                                                                        ccgcATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTgg                                         <  1:195299/149‑1 (MQ=255)
                                                                                                              cGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGcggccg                                   >  7:246861/1‑149 (MQ=255)
                                                                                                                 tCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCggccggcc                                >  4:41101/1‑149 (MQ=255)
                                                                                                                    gctgctgcCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCCCGGCTTCGCGCTGGCGGCCGgccgcc                             >  2:323459/1‑149 (MQ=255)
                                                                                                                    gctgctgcCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGgccgcc                             >  7:155401/1‑149 (MQ=255)
                                                                                                                    gctgctgcCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGgccgcc                             >  2:289156/1‑149 (MQ=255)
                                                                                                                       gctgcCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGCCCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCg                          >  1:244087/1‑149 (MQ=255)
                                                                                                                       gctgcCCTTCGCGGGCCTGCTGCTGGCCCCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCg                          >  1:213501/1‑149 (MQ=255)
                                                                                                                        ctgcCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCgg                         <  8:133494/149‑1 (MQ=255)
                                                                                                                          gcCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGcc                       <  8:252220/149‑1 (MQ=255)
                                                                                                                            ccTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCAt                     >  1:58050/1‑149 (MQ=255)
                                                                                                                              ttCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGt                   >  8:60233/1‑149 (MQ=255)
                                                                                                                                cgcgGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGt                   <  8:71352/147‑1 (MQ=255)
                                                                                                                                  cgGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACt               <  4:232532/149‑1 (MQ=255)
                                                                                                                                    ggCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGGGCCGACCATTTCGTGGAATGCGGTGCAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTgg             <  6:26351/149‑1 (MQ=255)
                                                                                                                                      ccTGCTGCTTGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTg              >  2:297570/1‑146 (MQ=255)
                                                                                                                                        tgctgctgGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGAc         >  6:246409/1‑149 (MQ=255)
                                                                                                                                          ctgctgGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGAcgc       <  4:167992/149‑1 (MQ=255)
                                                                                                                                          ctgctgGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGAcgc       <  5:316544/149‑1 (MQ=255)
                                                                                                                                          ctgctgGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGAcgc       <  2:315885/149‑1 (MQ=255)
                                                                                                                                          ctgctgGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGAcg        <  1:206114/148‑1 (MQ=255)
                                                                                                                                           tgctgGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGAcgcg      >  6:213490/1‑149 (MQ=255)
                                                                                                                                            gtcggCCGATAGCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCgg     >  4:267834/4‑149 (MQ=255)
                                                                                                                                             tcggCCGATAGCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCggg    >  3:60448/3‑149 (MQ=255)
                                                                                                                                                gCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTc  <  8:246861/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
GGTGCCGAGGATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTC  >  NC_002947/681824‑682115

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: