Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,560,597 A→C 6.6% N498K (AAT→AAG) ‡ pctA ← methyl‑accepting chemotaxis protein PctA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,560,5970AC 6.6% 189.9 / 4.5 61N498K (AAT→AAG) ‡pctAmethyl‑accepting chemotaxis protein PctA
Reads supporting (aligned to +/- strand):  ref base A (29/28);  new base C (2/2);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

GCTCACCCGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCTTGC  >  NC_002947/1560449‑1560745
                                                                                                                                                    |                                                                                                                                                    
gCTCACCGGCGTGGTGGGCAAGCTGCATGCGGTCCTGGCTGTGGCGCTGGCTCTGCGCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGGCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGa                                                                                                                                                      <  5:290614/149‑1 (MQ=255)
 cTCACCCGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGAt                                                                                                                                                     >  1:130712/1‑149 (MQ=255)
    aCCCGCCTGGGTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTgc                                                                                                                                                  >  6:120173/1‑149 (MQ=255)
      ccGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTgcgc                                                                                                                                                <  2:99789/149‑1 (MQ=255)
       cGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCa                                                                                                                                               <  6:63093/149‑1 (MQ=255)
        gCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCAc                                                                                                                                              >  4:180275/1‑149 (MQ=255)
         ccTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCAcc                                                                                                                                             >  3:329398/1‑149 (MQ=255)
         ccTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCAc                                                                                                                                              >  5:156897/1‑148 (MQ=255)
          cTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCAATCTGGCGCAccg                                                                                                                                            <  2:314329/149‑2 (MQ=255)
              ttGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATc                                                                                                                                        <  4:184531/149‑1 (MQ=255)
                 gCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGc                                                                                                                                     <  5:270571/149‑1 (MQ=255)
                   aaTCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGa                                                                                                                                   <  2:130712/149‑1 (MQ=255)
                       tGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCAc                                                                                                                               >  2:15502/1‑149 (MQ=255)
                         cATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACgg                                                                                                                             <  5:341833/149‑1 (MQ=255)
                            gCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCa                                                                                                                           >  1:309588/1‑148 (MQ=255)
                             cTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCaaa                                                                                                                         <  1:13528/149‑1 (MQ=255)
                                  cTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACgg                                                                                                                             >  1:138358/1‑140 (MQ=255)
                                  cTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACgg                                                                                                                             <  2:138358/140‑1 (MQ=255)
                                  cTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCAc                                                                                                                    <  4:329398/149‑1 (MQ=255)
                                     gCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACg                                                                                                                   <  7:186992/147‑1 (MQ=255)
                                     gCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGcc                                                                                                                 >  7:17122/1‑149 (MQ=255)
                                                 gCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCgg                                                                                                     >  4:136843/1‑149 (MQ=255)
                                                 gCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCgg                                                                                                     >  2:259045/1‑149 (MQ=255)
                                                   tctGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGgcg                                                                                                   >  3:261220/1‑149 (MQ=255)
                                                   tctGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGgcg                                                                                                   >  7:11716/1‑149 (MQ=255)
                                                   tctGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGgc                                                                                                    >  6:290372/1‑148 (MQ=255)
                                                     tGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGgc                                                                                                    <  5:290372/146‑1 (MQ=255)
                                                     tGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGcgcg                                                                                                 >  6:88037/1‑149 (MQ=255)
                                                       cTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCggg                                                                                               <  6:156897/149‑1 (MQ=255)
                                                          cATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGc                                                                                              <  2:364882/147‑1 (MQ=255)
                                                                 tCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGAt                                                                                     <  1:259045/149‑1 (MQ=255)
                                                                         cTCGCGGGCACCGCCTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCGTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCCAAGCCACGCCCGGCCTCCCCGGCGCGGGCCGCGTCGATGgcggcgg                                                                             >  8:244889/1‑148 (MQ=25)
                                                                            gcgGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGtt                                                                            >  5:358707/1‑147 (MQ=255)
                                                                                gCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAgcgc                                                                      <  5:88037/149‑1 (MQ=255)
                                                                                  aCCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCca                                                                    <  1:15502/149‑1 (MQ=255)
                                                                                  aCCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACAACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAgcgc                                                                      >  6:292528/1‑147 (MQ=255)
                                                                                    cGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCcagc                                                                  <  7:239211/149‑1 (MQ=255)
                                                                                    cGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGAGGTGCGCAATCTGGCGCACCGCa                                                                                                                                          <  3:329958/77‑1 (MQ=25)
                                                                                    cGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGAGGTGCGCAATCTGGCGCACCGCa                                                                                                                                          >  4:329958/1‑77 (MQ=25)
                                                                                        ttGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAggt                                                              >  5:349607/1‑149 (MQ=255)
                                                                                           cAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAggttg                                                            <  2:309588/148‑1 (MQ=255)
                                                                                              ctctTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGttt                                                        >  8:208988/1‑149 (MQ=255)
                                                                                               tctTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTg                                                       >  5:83322/1‑149 (MQ=255)
                                                                                                ctTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGc                                                      >  7:31812/1‑149 (MQ=255)
                                                                                                   cGATCATGGTCTTCACTTGCTGCGCCGACTCCTGGGTGCGGTTCGCCAGATTGCGCACCGCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCt                                                     <  8:175025/147‑1 (MQ=255)
                                                                                                   cGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCt                                                     <  8:31812/147‑1 (MQ=255)
                                                                                                    gATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTgg                                                   <  1:201378/148‑1 (MQ=255)
                                                                                                            tCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTg                                          <  3:180275/149‑1 (MQ=255)
                                                                                                             cTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTgg                                         >  3:138128/1‑149 (MQ=255)
                                                                                                                 aCTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGAt                                     <  8:11716/149‑1 (MQ=255)
                                                                                                                   ttGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATcc                                   >  3:280821/1‑148 (MQ=255)
                                                                                                                     gctgcGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCAc                                  <  4:261220/148‑1 (MQ=255)
                                                                                                                     gctgcGCCGACTCCTGGGTGCGGTGCGCAAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACt                                 >  3:99710/1‑149 (MQ=255)
                                                                                                                         cgcCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCa                             <  8:17122/149‑1 (MQ=255)
                                                                                                                           ccGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAg                            >  3:349635/1‑148 (MQ=255)
                                                                                                                           ccGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAg                            >  7:348833/1‑148 (MQ=255)
                                                                                                                              aCTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATc                        >  6:349376/1‑149 (MQ=255)
                                                                                                                                  cTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTgg                     <  3:136843/148‑1 (MQ=255)
                                                                                                                                     ggtgcggtgcGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGcc                   <  1:372082/147‑1 (MQ=255)
                                                                                                                                          ggtgcgCAATCTGGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATAtt             >  2:315305/1‑148 (MQ=255)
                                                                                                                                                 cAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCt     <  7:208988/149‑1 (MQ=255)
                                                                                                                                                 cAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCt     >  2:368458/1‑149 (MQ=255)
                                                                                                                                                    aTTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCTTGc  >  4:139506/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCTCACCCGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCTTGC  >  NC_002947/1560449‑1560745

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: