Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 681,975 T→G 5.9% D433A (GAT→GCT)  PP_0584 ← methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947681,9750TG 5.9% 265.5 / 6.5 85D433A (GAT→GCT) PP_0584methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (44/36);  new base G (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 = 8.63e-01

GCCGAGGATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTCGAGCTGCT  >  NC_002947/681827‑682123
                                                                                                                                                    |                                                                                                                                                     
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                                                                                       <  5:99505/149‑1 (MQ=255)
                                                               cTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGt                                                                                       <  5:167950/149‑1 (MQ=255)
                                                               cTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGt                                                                                       >  3:140401/1‑149 (MQ=255)
                                                               cTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGt                                                                                       <  2:230819/149‑1 (MQ=255)
                                                                tCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCgg                                                                                                      <  6:285821/133‑1 (MQ=255)
                                                                tCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCgg                                                                                                      >  5:285821/1‑133 (MQ=255)
                                                                tCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTg                                                                                      >  5:110447/1‑149 (MQ=255)
                                                                     ttCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATg                                                                                 >  7:252220/1‑149 (MQ=255)
                                                                         tGGCTAACGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGt                                                                             <  6:193354/149‑1 (MQ=255)
                                                                             tAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGa                                                                          >  2:279569/1‑148 (MQ=255)
                                                                             tAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGa                                                                          >  7:138468/1‑148 (MQ=255)
                                                                                ggCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCACCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAAc                                                                       >  3:326834/1‑148 (MQ=255)
                                                                                 gCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCa                                                                     >  4:15356/1‑149 (MQ=255)
                                                                                   gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGc                                                                                                  <  2:57334/118‑1 (MQ=255)
                                                                                   gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGc                                                                                                  >  1:57334/1‑118 (MQ=255)
                                                                                   gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATg                                                                   >  6:316544/1‑149 (MQ=255)
                                                                                   gATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATg                                                                   >  8:76798/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                                                                                                                      >  2:353441/1‑80 (MQ=255)
                                                                                                     ccgcATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGt                                                                                                                      <  1:353441/80‑1 (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                                     >  2:289156/1‑149 (MQ=255)
                                                                                                                 gctgctgcCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGgccgcc                                     >  7:155401/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               <  2:315885/149‑1 (MQ=255)
                                                                                                                                       ctgctgGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGAcgc               <  5:316544/149‑1 (MQ=255)
                                                                                                                                       ctgctgGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGAcgc               <  4:167992/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)
                                                                                                                                                  tATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTCGAgctg    >  4:371071/1‑149 (MQ=255)
                                                                                                                                                    tCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTCGAgctgct  <  3:41101/149‑1 (MQ=255)
                                                                                                                                                    tCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTCGAgctgc   <  7:60233/148‑1 (MQ=255)
                                                                                                                                                    tCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTCGAgctg    <  2:213501/147‑1 (MQ=255)
                                                                                                                                                    tCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATATCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTCGAgctgct  >  6:384400/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCCGAGGATCGACTGAATGGCGTTGCTGTCACGCTCCAGCTGTTGGATCGACTGGGCCGATTGCTCGATTTCGTGGCTAAGGCGATCAACGCTGTGTACCGCCGCATCGATCTGCTGCTGCCCTTCGCGGGCCTGCTGCTGGCCGCTATCGGCCGACTGCGCCGCCTGGCTGCAGGACCGTGCCACTTCATTGGCGGTGGCGACCATTTCGTGGAATGCGGTGGAAACCATGTCCACGGCTTCGCGCTGGCGGCCGGCCGCCTCGGCCATGTCACTGGAGACGCGGGTCGAGCTGCT  >  NC_002947/681827‑682123

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: