Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,560,599 T→G 6.6% N498H (AAT→CAT) ‡ 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,5990TG 6.6% 195.7 / 4.6 61N498H (AAT→CAT) ‡pctAmethyl‑accepting chemotaxis protein PctA
Reads supporting (aligned to +/- strand):  ref base T (29/28);  new base G (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.94e-01

ACCCGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCTTGC  >  NC_002947/1560453‑1560745
                                                                                                                                                  |                                                                                                                                                  
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                                                                                                     >  2:259045/1‑149 (MQ=255)
                                             gCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCgg                                                                                                     >  4:136843/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                                                                                                                                          >  4:329958/1‑77 (MQ=25)
                                                                                cGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGAGGTGCGCAATCTGGCGCACCGCa                                                                                                                                          <  3:329958/77‑1 (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)
                                                                                                                                                  |                                                                                                                                                  
ACCCGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCTTGC  >  NC_002947/1560453‑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: