Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,560,596 G→T 6.4% L499M (CTG→ATG) ‡ 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,5960GT 6.4% 201.5 / 3.9 62L499M (CTG→ATG) ‡pctAmethyl‑accepting chemotaxis protein PctA
Reads supporting (aligned to +/- strand):  ref base G (29/29);  new base T (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.86e-01

GCTCACCCGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCT  >  NC_002947/1560449‑1560742
                                                                                                                                                   |                                                                                                                                                  
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                                                                                                                          <  2:138358/140‑1 (MQ=255)
                                  cTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACgg                                                                                                                          >  1:138358/1‑140 (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  >  2:368458/1‑149 (MQ=255)
                                                                                                                                                 cAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCt  <  7:208988/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTCACCCGCCTGGTTGGCAATCTGCATGCTGTCCTGGCTGTGGCGCTGGCTCTGCTCCATGGTGTCGACCGACTCGCGGGCACCGACTTGCAGCTCTTCGATCATGGTCTGCACTTGCTGCGCCGACTCCTGGGTGCGGTGCGCCAGATTGCGCACCTCATCGGCGACCACGGCAAAGCCACGCCCGGCCTCACCGGCGCGGGCCGCTTCGATGGCGGCGTTCAGCGCCAGCAGGTTGGTTTGCTGGGAAATGCTGGTGATCACTTCCAGGATCTGGCCGATATTCACGGTCT  >  NC_002947/1560449‑1560742

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: