Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,782,835 G→C 6.5% S437T (AGC→ACC)  kefA → intermediate mechanosensitive channel protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,782,8350GC 6.5% 201.0 / 2.9 61S437T (AGC→ACC) kefAintermediate mechanosensitive channel protein
Reads supporting (aligned to +/- strand):  ref base G (32/25);  new base C (2/2);  total (34/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

TGACACCGCAGCTGCGCAAGGCCCTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGCCGATGGC  >  NC_002947/5782688‑5782983
                                                                                                                                                    |                                                                                                                                                    
tgACACCGCAGCTGCGCAAGGCCCTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGc                                                                                                                                                     <  4:272408/149‑1 (MQ=255)
 gACACCGCAGCTGCGCAAGGCCCTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGcc                                                                                                                                                    >  2:310740/1‑149 (MQ=255)
      cGCAGCTGCGCAAGGCCCTGCTGGAAGTGGCCATCACGCGCAGCGACCTGGTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCTATCACCTTGCAGCTCAACCAGAAGCATCTGCTGGGCACCGCGCAGAGCCTgcgc                                                                                                                                               <  5:26462/149‑1 (MQ=255)
         aGCTGCGCAAGGCCCTGCTGGAAGTGGTCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGGCTCTgcgcggt                                                                                                                                            <  6:1850/149‑4 (MQ=255)
               gcAAGGCCCTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTc                                                                                                                                      >  2:15036/1‑149 (MQ=255)
                   tctcgtgcttGTAGTGGCCATCAACCTCAGCGACCTGTTGTAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTcgac                                                                                                                                   <  6:28271/143‑1 (MQ=255)
                    gCCCTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTcgacga                                                                                                                                 <  1:120622/149‑1 (MQ=255)
                    gCCCTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTcgacga                                                                                                                                 <  3:260354/149‑1 (MQ=255)
                     cccTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACAAg                                                                                                                                <  1:310740/149‑1 (MQ=255)
                      ccTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGGCTCTgcgcggt                                                                                                                                            <  1:318117/136‑4 (MQ=255)
                      ccTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGGCTCTgcgcggt                                                                                                                                            >  2:318117/1‑133 (MQ=255)
                      ccTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAg                                                                                                                                >  4:333259/1‑148 (MQ=255)
                      ccTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAg                                                                                                                                >  2:371313/1‑148 (MQ=255)
                      ccTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGc                                                                                                                               <  3:148424/149‑1 (MQ=255)
                       ctgctgGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTcgac                                                                                                                                   >  5:28271/1‑144 (MQ=255)
                              aaGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTc                                                                                                                       <  6:125895/149‑1 (MQ=255)
                               aGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCt                                                                                                                      <  1:371313/149‑1 (MQ=255)
                                     cATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCAGCTCAACGAATCGAGAAACCTTGCAGCTCAACCACAAGCAGCTGCTGGGCACCGCGCAGACCCTGCGCACCACCGTCCACGAGCAAATGTTCTGGATc                                                                                                                 >  1:110602/1‑149 (MQ=255)
                                         aCCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCa                                                                                                             <  6:294503/148‑1 (MQ=255)
                                           ccGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGc                                                                                                           >  4:330397/1‑148 (MQ=255)
                                           ccGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGc                                                                                                           >  5:92204/1‑148 (MQ=255)
                                               aGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGcaac                                                                                                        >  5:224908/1‑147 (MQ=255)
                                                   aCCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCg                                                                                                  >  5:99530/1‑149 (MQ=255)
                                                     cTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCt                                                                                                >  2:263183/1‑149 (MQ=255)
                                                     cTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCt                                                                                                >  4:16304/1‑149 (MQ=255)
                                                           gAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactg                                                                                          >  3:271379/1‑149 (MQ=255)
                                                           gAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactg                                                                                          >  6:216904/1‑149 (MQ=255)
                                                            aGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactgg                                                                                         <  3:16304/149‑1 (MQ=255)
                                                                gCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactgggact                                                                                     >  5:182008/1‑149 (MQ=255)
                                                                  tGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGCGCCTGCGCACCACCCTCTACGAGCAAATGTTCTGGATCCCCAGCAACCCGCCGCGGgactgggactgg                                                                                   >  8:337404/1‑149 (MQ=255)
                                                                  tGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGCTCCCCAGCAACAAGCCGCTGgactgggactgg                                                                                   >  3:74508/1‑149 (MQ=255)
                                                                   gAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactgg                                                                                         >  3:366257/1‑142 (MQ=255)
                                                                   gAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactgg                                                                                         <  4:366257/142‑1 (MQ=255)
                                                                        gcgAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTAGACGAGCAAATGTTCTGTATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCaa                                                                             <  2:110602/149‑1 (MQ=255)
                                                                          gAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATa                                                                           <  1:15036/149‑1 (MQ=255)
                                                                             cTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGt                                                                        <  8:191888/149‑1 (MQ=255)
                                                                                tCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactgg                                                                                         <  7:175024/129‑1 (MQ=255)
                                                                                tCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactgg                                                                                         >  8:175024/1‑129 (MQ=255)
                                                                                tCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGcc                                                                     >  3:332498/1‑149 (MQ=255)
                                                                                   gcgcTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCg                                                                    >  1:255702/1‑147 (MQ=255)
                                                                                    cgcTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAg                                                                 <  6:182008/149‑1 (MQ=255)
                                                                                     gctgctCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAgc                                                                >  4:297508/1‑149 (MQ=255)
                                                                                         ctCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCt                                                            >  8:277716/1‑149 (MQ=255)
                                                                                         ctCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCt                                                            >  6:276153/1‑149 (MQ=255)
                                                                                            aaCGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTg                                                           <  6:224908/147‑1 (MQ=255)
                                                                                            aaCGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGc                                                         >  4:67715/1‑149 (MQ=255)
                                                                                             aCGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTgtactgg                                                                                         <  7:335351/116‑1 (MQ=255)
                                                                                             aCGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTggactgg                                                                                         >  8:335351/1‑116 (MQ=255)
                                                                                             aCGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGcc                                                        <  6:99530/149‑1 (MQ=255)
                                                                                                  tCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGAc                                                     >  5:290266/1‑147 (MQ=255)
                                                                                                  tCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCa                                                   <  3:333259/149‑1 (MQ=255)
                                                                                                                     aaCCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATg                                 >  2:112492/1‑149 (MQ=255)
                                                                                                                         aGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGACGCCAACTTGCCATGGGGc                             >  5:127948/1‑149 (MQ=255)
                                                                                                                          gAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCt                            >  1:299084/1‑149 (MQ=255)
                                                                                                                             gcagcTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTc                           >  3:277618/1‑147 (MQ=255)
                                                                                                                                      gggCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGctcttt                  >  1:247452/1‑142 (MQ=255)
                                                                                                                                       ggCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGaaa                 <  2:299084/147‑1 (MQ=255)
                                                                                                                                       ggCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGaaa                 <  7:277716/147‑1 (MQ=255)
                                                                                                                                       ggCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGaaa                 <  5:276153/147‑1 (MQ=255)
                                                                                                                                          aCCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCt            <  4:74508/149‑1 (MQ=255)
                                                                                                                                             gcgcAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGc         >  1:20260/1‑149 (MQ=255)
                                                                                                                                                 agagCCTGCGAACCACCCTCGACGAGCAAATGTTCTGGAGCCCCCGCAACTAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGCCGAt     <  4:271379/149‑1 (MQ=255)
                                                                                                                                                  gagCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGCCGATg    >  1:206085/1‑149 (MQ=255)
                                                                                                                                                   agCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGCCGATgg   >  4:315130/1‑149 (MQ=255)
                                                                                                                                                   agCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCCGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGCCg       <  1:161560/145‑1 (MQ=255)
                                                                                                                                                   agCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCAGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGCCg       >  2:161560/1‑145 (MQ=255)
                                                                                                                                                    gCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGCCGATGGc  <  2:255702/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TGACACCGCAGCTGCGCAAGGCCCTGCTGGAAGTGGCCATCACCCGCAGCGACCTGTTGGAGCGGCTGAACCGCGAGCTGTCGGCGCTGCTCAACGAATCGATCACCTTGCAGCTCAACCAGAAGCAGCTGCTGGGCACCGCGCAGAGCCTGCGCACCACCCTCGACGAGCAAATGTTCTGGATCCCCAGCAACAAGCCGCTGGACTGGGACTGGCTGCAATACGTGCCGAAGCGCCTGGCCGACCAGGTCGCCAACTTGCCATGGGGCTCCAGCCTGAAAGAGCTGGCCGATGGC  >  NC_002947/5782688‑5782983

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: