Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,565,106 T→G 6.2% Q630P (CAG→CCG) ‡ pctB ← methyl‑accepting chemotaxis protein PctB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,565,1060TG 6.2% 284.2 / 9.5 98Q630P (CAG→CCG) ‡pctBmethyl‑accepting chemotaxis protein PctB
Reads supporting (aligned to +/- strand):  ref base T (48/43);  new base G (2/4);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

TTGCTGGCCTGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGTTCTTCGGTCGCTGTCGCCACCGACTGGTTCAT  >  NC_002947/2564959‑2565254
                                                                                                                                                    |                                                                                                                                                    
ttGCTGGCCTGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCctgc                                                                                                                                                    <  2:236189/150‑1 (MQ=255)
   cTGGCCTGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTgg                                                                                                                                                 >  1:1691559/1‑150 (MQ=255)
   cTGGCCTGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTgg                                                                                                                                                 <  1:940872/150‑1 (MQ=255)
    tGGCCTGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTggg                                                                                                                                                >  1:883159/1‑150 (MQ=255)
     ggCCTGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGc                                                                                                                                               <  1:1830006/150‑1 (MQ=255)
        cTGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCAGGCTGGGcgtt                                                                                                                                            <  2:605551/150‑4 (MQ=255)
        cTGGGTAAGGCCGGGCAGGGCGGCTAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAg                                                                                                                                            >  2:594628/1‑150 (MQ=255)
         tGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCAGGCTGGGcgttt                                                                                                                                           <  1:1147793/150‑5 (MQ=255)
            gTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTgg                                                                                                                                        >  2:333418/1‑150 (MQ=255)
             tACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGt                                                                                                                                       <  1:1369739/150‑1 (MQ=255)
                ggcggGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCAGGCTGGGcgtttga                                                                                                                                         <  1:1894918/145‑7 (MQ=255)
                ggcggGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCAGGCTGGGcgtttga                                                                                                                                         >  2:1894918/1‑139 (MQ=255)
                 gcggGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCa                                                                                                                                   <  1:1167984/150‑1 (MQ=255)
                      gCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGGACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGgc                                                                                                                              >  1:1948959/1‑150 (MQ=255)
                         gggCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTCCAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCa                                                                                                                           <  2:307425/150‑1 (MQ=255)
                                      ggTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCGATTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAg                                                                                                              <  1:2246568/150‑1 (MQ=255)
                                       gTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAgg                                                                                                             >  1:848751/1‑150 (MQ=255)
                                          gcggcATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAg                                                                                                              >  1:847888/1‑146 (MQ=255)
                                          gcggcATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAg                                                                                                              <  2:847888/146‑1 (MQ=255)
                                            ggcATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGTTGGGCTAGCTGTTCTCAATGTCGCAGGGTCTCGTTGAGGTtctc                                                                                                        <  2:2082231/150‑1 (MQ=255)
                                            ggcATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTtctc                                                                                                        >  2:1786635/1‑150 (MQ=255)
                                               aTAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCAGGCTGGGcgtt                                                                                                                                            >  2:832482/1‑108 (MQ=255)
                                               aTAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCAGGCTGGGcgtt                                                                                                                                            <  1:832482/111‑4 (MQ=255)
                                                tAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCAcc                                                                                                    >  2:2006052/1‑150 (MQ=255)
                                                    gggCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGAAGCACTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCct                                                                                                >  1:263500/1‑150 (MQ=255)
                                                      gCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCctgc                                                                                              <  1:1786635/150‑1 (MQ=255)
                                                        tGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCctgctg                                                                                            <  1:2006052/150‑1 (MQ=255)
                                                              aGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAg                                                                                      <  1:425360/150‑1 (MQ=255)
                                                                gCGACAACAGGCATAGGCCGTTGCGGCCTGAAACCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAA‑GCCCAGCCTGCTGGGCCTGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAggg                                                                                    >  1:1046067/1‑150 (MQ=255)
                                                                gCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGc                                                                                                                                               >  2:1127693/1‑91 (MQ=255)
                                                                gCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGc                                                                                                                                               <  1:1127693/91‑1 (MQ=255)
                                                                  gACAACAGGCATAGGGCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAggg                                                                                    <  2:460210/148‑1 (MQ=255)
                                                                  gACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATNCTGANGCTGCCAACCACCTCATTCCCTCTTAGCAGGCCGC‑CGAACAGGCGTTGACAATGGCGCAGGGTCTCGTTGAGGTTCTACCCCCACCGCGGGGTCAGGGCg                                                                                  >  1:172756/1‑150 (MQ=255)
                                                                  gACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAggg                                                                                    >  1:460210/1‑148 (MQ=255)
                                                                   acaacaGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGt                                                                                 <  1:1422982/150‑1 (MQ=255)
                                                                     aacaGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCATGGTCTCGTTGa                                                                                                               <  2:1908821/118‑1 (MQ=255)
                                                                     aacaGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGa                                                                                                               >  1:1908821/1‑118 (MQ=255)
                                                                     aacaGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTg                                                                               >  1:1927487/1‑150 (MQ=255)
                                                                       caGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTTAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGAt                                                                             <  2:883159/150‑1 (MQ=255)
                                                                       caGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGc                                                                                   >  2:86689/1‑144 (MQ=255)
                                                                       caGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGc                                                                                   <  1:86689/144‑1 (MQ=255)
                                                                        aGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATc                                                                            >  2:1487408/1‑150 (MQ=255)
                                                                            aTAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTggg                                                                        >  2:2043155/1‑150 (MQ=255)
                                                                            aTAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTggg                                                                        <  1:594628/150‑1 (MQ=255)
                                                                             tAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGt                                                                       <  2:1948959/150‑1 (MQ=255)
                                                                             tAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGt                                                                       >  1:392521/1‑150 (MQ=255)
                                                                               ggCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGa                                                                     >  2:829062/1‑150 (MQ=255)
                                                                                gCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTc                                                                                                                                      <  2:1336126/84‑1 (MQ=255)
                                                                                gCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTc                                                                                                                                      >  1:1336126/1‑84 (MQ=255)
                                                                                gCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGAt                                                                    >  1:1366697/1‑150 (MQ=255)
                                                                                gCCGTTGCCGCCGGAACCCGGCGTTTAAATGCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCACCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCGGCTGGTTCAGGGCGTTGATCTGGGTGAt                                                                    <  2:1014256/150‑1 (MQ=255)
                                                                                accGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGAt                                                                    >  2:1902315/2‑150 (MQ=255)
                                                                                    ttGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTcc                                                                <  1:1487408/150‑1 (MQ=255)
                                                                                      gccgccTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAg                                                              <  1:518069/150‑1 (MQ=255)
                                                                                           cTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGa                                                         <  1:2043155/150‑1 (MQ=255)
                                                                                            tGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTgg                                                                         >  1:1784103/1‑133 (MQ=255)
                                                                                            tGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTgg                                                                         <  2:1784103/133‑1 (MQ=255)
                                                                                             gAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAgg                                                       <  2:848751/150‑1 (MQ=255)
                                                                                               aCCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGct                                                     >  2:1941085/1‑150 (MQ=255)
                                                                                                 ccGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGtc                                                                                                                      <  1:1698285/83‑1 (MQ=255)
                                                                                                 ccGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGtc                                                                                                                      >  2:1698285/1‑83 (MQ=255)
                                                                                                  cGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCtctc                                                  >  2:1880861/1‑150 (MQ=255)
                                                                                                   ggCGTTTAAATCCTGAAGCTGCCCACCAACTTCTTCAAAC‑GCCCAGCCTACTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCg                                                 <  2:1046067/150‑1 (MQ=255)
                                                                                                   ggCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATATTGAGGTTGAGGCTCTCg                                                 >  1:2044737/1‑150 (MQ=255)
                                                                                                    gCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGa                                                >  2:2281011/1‑150 (MQ=255)
                                                                                                        ttAAATCCTGACGCTGCGCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCAc                                            <  1:1062559/150‑1 (MQ=255)
                                                                                                             tCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTtctc                                                                                                        <  2:82057/85‑1 (MQ=255)
                                                                                                             tCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTtctc                                                                                                        >  1:82057/1‑85 (MQ=255)
                                                                                                               cTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCCGCATTCGGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGATCGCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTc                                     >  1:1398266/1‑150 (MQ=255)
                                                                                                               cTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCa                                                                                                                           <  2:1544296/64‑1 (MQ=255)
                                                                                                               cTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCa                                                                                                                           >  1:1544296/1‑64 (MQ=255)
                                                                                                               cTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTc                                     >  2:1056845/1‑150 (MQ=255)
                                                                                                               cTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTc                                     <  1:1902315/150‑1 (MQ=255)
                                                                                                                   aGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGtt                                 <  1:333418/150‑1 (MQ=255)
                                                                                                                       gCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGGGCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGttcttc                             <  2:586126/150‑1 (MQ=255)
                                                                                                                        cccACCAACTGCTTCAAAC‑GCCCAGCCTGCTTGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTggg                                                                        <  1:898579/106‑1 (MQ=255)
                                                                                                                        cccACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTggg                                                                        >  2:898579/1‑106 (MQ=255)
                                                                                                                        cccACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAgg                                                             <  1:1013280/117‑1 (MQ=255)
                                                                                                                        cccACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAgg                                                             >  2:1013280/1‑117 (MQ=255)
                                                                                                                         ccaccaACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCAcggcgg                                       >  2:2168680/1‑138 (MQ=255)
                                                                                                                         ccaccaACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCAcggcgg                                       <  1:2168680/138‑1 (MQ=255)
                                                                                                                                  gCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTtctc                                                                                                        <  1:1504264/64‑1 (MQ=255)
                                                                                                                                  gCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTtctc                                                                                                        >  2:1504264/1‑64 (MQ=255)
                                                                                                                                   cTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCg                                                 >  2:177210/1‑118 (MQ=255)
                                                                                                                                   cTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCg                                                 <  1:177210/118‑1 (MQ=255)
                                                                                                                                       aaaC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGtt                                                           <  2:503561/104‑1 (MQ=255)
                                                                                                                                       aaaC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGtt                                                           >  1:503561/1‑104 (MQ=255)
                                                                                                                                       aaaC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGa                                                         <  2:2343799/106‑1 (MQ=255)
                                                                                                                                       aaaC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGa                                                         >  1:2343799/1‑106 (MQ=255)
                                                                                                                                       aaaC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGttctt                              >  1:2242674/1‑133 (MQ=255)
                                                                                                                                       aaaC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGttctt                              <  2:2242674/133‑1 (MQ=255)
                                                                                                                                       aaaC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGGTCTCCACCCCCTGCTGGGTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGTTCTTCGGTCGCTGTCGCCAcc             >  1:2348931/1‑150 (MQ=255)
                                                                                                                                       aaaC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTATCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGCTCTGGGTGATGTCAAGGTTGAGGCTCTCGACCACGGCGGTAGGATCGTCGGTATCTGTCGCGAcc             >  2:2316720/1‑150 (MQ=255)
                                                                                                                                        ttgcgaccAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGTTCTTCGGTCGCTGTCGCCAcc             >  1:2371497/7‑150 (MQ=255)
                                                                                                                                             cccAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGTTCTTCGGTCGCTGTCGCCACCGACTg        >  2:1240461/1‑150 (MQ=255)
                                                                                                                                              ccAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGTTCTTCGGTCGCTGTCGCCACCGACTgg       >  2:1522881/1‑150 (MQ=255)
                                                                                                                                               cAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGTTCTTCGGTCGCTGTCGCCACCGACTGGt      >  2:1711457/1‑150 (MQ=255)
                                                                                                                                                aGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGt                                                                                          >  2:2070436/1‑65 (MQ=255)
                                                                                                                                                aGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGt                                                                                          <  1:2070436/65‑1 (MQ=255)
                                                                                                                                                   ctgctgGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGTTCTTCGGTCGCTGTCGCCACCGACTGGTTCAt  <  1:599871/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGCTGGCCTGGGTACGGCGGGGCAGGGCTGCGAAGGCGGTGGCGGCATAAGGGGCTGGTACAGGCGACAACAGGCATAGGCCGTTGCCGCCTGAACCCGGCGTTTAAATCCTGAAGCTGCCCACCAACTGCTTCAAAC‑GCCCAGCCTGCTGGGCCAGCTGGTCGCAATGGCGCAGGGTCTCGTTGAGGTTCTCCACCCCCTGCTGGTTCAGGGCGTTGATCTGGGTGATGTCCAGGTTGAGGCTCTCGACCACGGCGGTCTGTTCTTCGGTCGCTGTCGCCACCGACTGGTTCAT  >  NC_002947/2564959‑2565254

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