Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,063,466 G→A 7.1% G119S (GGC→AGC)  ubiC → chorismate pyruvate‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,063,4660GA 7.1% 210.2 / 4.6 70G119S (GGC→AGC) ubiCchorismate pyruvate‑lyase
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base A (2/3);  total (37/33)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

ATCGCTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACCTGCTGGT  >  NC_002947/6063322‑6063614
                                                                                                                                                |                                                                                                                                                    
aTCGCTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGa                                                                                                                                                  <  5:51012/149‑1 (MQ=255)
 tCGCTGCCGGCTCCGAAGGCTGGGTGCTCGAAGTGTATCTGCGCGGCCATGGCAAACCTTGTGTATTCTCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTTGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAg                                                                                                                                                 <  3:52939/149‑1 (MQ=255)
    cTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGGCAACCTTGTGTATTCGCACGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAgctg                                                                                                                                              <  8:258954/149‑1 (MQ=255)
      gCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCCCAGCGAgcggct                                                                                                                                            <  5:43090/149‑1 (MQ=255)
           cgcCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCt                                                                                                                                       <  5:2664/149‑1 (MQ=255)
           cgcCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCt                                                                                                                                       <  5:272622/149‑1 (MQ=255)
             ccGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCAGCTCGCCCAGCGAGcggctgcccaag                                                                                                                                     <  3:238678/149‑7 (MQ=255)
               gAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGa                                                                                                                                   >  8:296578/1‑149 (MQ=255)
                  ggCTTGGTGCGGGAAGTTTATCTGCGCGGCCATGGCCAACCTTGTGTATTCGCCCTCATCGTTGCCTGCCGCAGCGCCCTTGAACGTGGTGGCCTTGACCTGGAAACCTTTGGCAGCCGGTCGCTGGGCGTGCTGCTGTTCTGCGTcc                                                                                                                                 <  4:15625/148‑1 (MQ=255)
                       ggTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGAc                                                                                                                                  >  8:245578/1‑142 (MQ=255)
                       ggTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGAc                                                                                                                                  <  7:245578/142‑1 (MQ=255)
                         tGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGAGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGtt                                                                                                                          >  6:85319/1‑148 (MQ=255)
                              gAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCg                                                                                                                    >  7:269532/1‑149 (MQ=255)
                              gAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCg                                                                                                                    <  7:73597/149‑1 (MQ=255)
                                aGTGTATCTGCGCGGCCATGGACAACCTTGGGTATTCGCCCGGAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTc                                                                                                                  <  8:247672/149‑1 (MQ=255)
                                     aTCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCa                                                                                                             >  5:301317/1‑149 (MQ=255)
                                      tCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGt                                                                                                                           >  1:329495/1‑134 (MQ=255)
                                      tCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGt                                                                                                                           <  2:329495/134‑1 (MQ=255)
                                      tCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCAGCTCGCCCAGCGAgcggct                                                                                                                                            <  8:132171/117‑1 (MQ=255)
                                      tCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCAGCTCGCCCAGCGAgcggct                                                                                                                                            >  7:132171/1‑117 (MQ=255)
                                        tGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAACCCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGa                                                                                                                                   >  7:119576/1‑124 (MQ=255)
                                         gcgcGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTc                                                                                                          >  5:183170/1‑148 (MQ=255)
                                         gcgcGGCCATGGCCAACCTTGGGTATTCGCCCGCAACGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTc                                                                                                          >  4:276106/1‑148 (MQ=255)
                                          cgcgGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCAc                                                                                                            <  7:296578/145‑1 (MQ=255)
                                          cgcgGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTc                                                                                                          <  2:229532/147‑1 (MQ=255)
                                             ggCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAgt                                                                                                     >  3:71633/1‑149 (MQ=255)
                                             ggCCATGGCCAACCTTGGGTATTCGCCCGCAACGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAgt                                                                                                     <  3:276106/149‑1 (MQ=255)
                                                cATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGtgtg                                                                                                  >  1:42224/1‑149 (MQ=255)
                                                 aTGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGGCCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGtgtg                                                                                                  >  7:66671/1‑148 (MQ=255)
                                                   ggCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAgt                                                                                                     >  8:294496/1‑143 (MQ=255)
                                                   ggCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAgt                                                                                                     <  7:294496/143‑1 (MQ=255)
                                                   ggCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCAc                                                                                               >  3:33843/1‑149 (MQ=255)
                                                         ccTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATc                                                                                          >  5:131414/1‑148 (MQ=255)
                                                          cTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCca                                                                                        >  7:89974/1‑149 (MQ=255)
                                                           ttGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATc                                                                                                               <  6:57483/125‑1 (MQ=255)
                                                           ttGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATc                                                                                                               >  5:57483/1‑125 (MQ=255)
                                                              ggTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAgg                                                                                    >  6:155378/1‑149 (MQ=255)
                                                                tATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGTGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGcc                                                                                  >  8:285817/1‑149 (MQ=255)
                                                                    cGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTgg                                                                               <  4:33843/148‑1 (MQ=255)
                                                                      cccGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTg                                                                            <  1:59784/149‑1 (MQ=255)
                                                                      cccGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTg                                                                            <  1:54864/149‑1 (MQ=255)
                                                                        cGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGcc                                                                          >  6:184152/1‑149 (MQ=255)
                                                                        cGCAGAGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGcc                                                                          <  5:85319/149‑1 (MQ=255)
                                                                         gCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTAt                                                                                           >  8:196920/1‑131 (MQ=255)
                                                                         gCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTAt                                                                                           <  7:196920/131‑1 (MQ=255)
                                                                         gCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCg                                                                         <  5:90321/149‑1 (MQ=255)
                                                                           aGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTc                                                                       <  6:301317/149‑1 (MQ=255)
                                                                             cGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGGCCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCg                                                                     <  8:66671/149‑1 (MQ=255)
                                                                               tGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGaa                                                                   <  5:155378/149‑1 (MQ=255)
                                                                                  ccAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCt                                                                <  6:131414/149‑1 (MQ=255)
                                                                                      ccGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGcac                                                            <  4:325441/149‑1 (MQ=255)
                                                                                        gCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACAt                                                          >  3:137241/1‑149 (MQ=255)
                                                                                          aGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATgc                                                        <  6:183170/149‑1 (MQ=255)
                                                                                           gcgcCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGCCCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCAAATgcg                                                       >  2:294851/1‑149 (MQ=255)
                                                                                               ccTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGgcgc                                                   >  7:220127/1‑149 (MQ=255)
                                                                                               ccTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGgcgc                                                   >  8:62439/1‑149 (MQ=255)
                                                                                               ccTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGgcgc                                                   >  7:312666/1‑149 (MQ=255)
                                                                                                 tGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCtt                                                 >  8:298256/1‑149 (MQ=255)
                                                                                                    aaCGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTgg                                              <  5:184152/149‑1 (MQ=255)
                                                                                                    aaCGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTgg                                              <  4:71633/149‑1 (MQ=255)
                                                                                                       gtggtgGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGgc                                           >  6:99397/1‑149 (MQ=255)
                                                                                                                           aCCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAgcg                       >  8:102200/1‑149 (MQ=255)
                                                                                                                              ttGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCaa                    >  2:82029/1‑149 (MQ=255)
                                                                                                                              ttGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCaa                    >  7:68162/1‑149 (MQ=255)
                                                                                                                              cagaacagcagCTCGCCCAGCGAGCGGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCaa                    >  5:134217/6‑149 (MQ=255)
                                                                                                                                      ccGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTg             >  5:57172/1‑148 (MQ=255)
                                                                                                                                        gctcgctGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGAc          >  2:79166/1‑149 (MQ=255)
                                                                                                                                            gctGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACctgc      <  8:220127/149‑1 (MQ=255)
                                                                                                                                              tGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACctgctg    >  6:294007/1‑149 (MQ=255)
                                                                                                                                              tGGGCGAGCTGCTGTTCTGCGACAAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTCGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACctgctg    >  2:153279/1‑149 (MQ=255)
                                                                                                                                                ggCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACCTGCtggt  <  1:82029/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
ATCGCTGCCGGCGCCGAAGGCTGGGTGCGCGAAGTGTATCTGCGCGGCCATGGCCAACCTTGGGTATTCGCCCGCAGCGTGGCCAGCCGCAGCGCCCTGGAACGTGGTGGCCTGGACCTGGAAACCTTGGGCAGCCGCTCGCTGGGCGAGCTGCTGTTCTGCGACCAGGCGTTCATCCGTCATCCACTCGAAGTGTGCACTTATCCACAGGCCTGGCTGCCGTCCGAAGCTGCACATGCGGCGCTTTGGGGCCGCCGCTCGCGCTTCGAGCGCAACGGCCTGGACCTGCTGGT  >  NC_002947/6063322‑6063614

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: