Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,864 T→G 6.6% intergenic (+115/‑66) gapB → / → PP_2150 glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,456,8640TG 6.6% 205.8 / 7.4 77intergenic (+115/‑66)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base G (3/2);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.44e-01

CTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAC  >  NC_002947/2456718‑2457012
                                                                                                                                                   |                                                                                                                                                     
cTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCg                                                                                                                                                    <  2:91600/149‑1 (MQ=255)
   gtgtTAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGc                                                                                                                                                   <  1:215758/147‑1 (MQ=255)
    tgtgAACCCGCCAGCGTTTCGACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGc                                                                                                                                                <  2:247975/149‑1 (MQ=255)
      tgAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGccc                                                                                                                                              <  8:306672/149‑1 (MQ=255)
           ccGCCAGCGTTTCCACGCTGATTCGCTTGTAATGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTa                                                                                                                                         <  7:34797/149‑1 (MQ=255)
           ccGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTg                                                                                                                                           >  2:213838/1‑147 (MQ=255)
            cGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAg                                                                                                                                        <  4:294121/149‑1 (MQ=255)
                   gTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCgg                                                                                                                                  >  3:188059/1‑148 (MQ=37)
                   gTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCgg                                                                                                                                  <  4:188059/148‑1 (MQ=37)
                         aCGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAcc                                                                                                                           <  4:42629/149‑1 (MQ=255)
                         aCGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTTcc                                                                                                                           >  7:258941/1‑149 (MQ=255)
                              gATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                  >  2:154487/1‑137 (MQ=255)
                              gATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                                  <  1:154487/137‑1 (MQ=255)
                                ttCGCTTGTAAGGCAATGAAAAACCGGCAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                     <  7:319158/148‑1 (MQ=255)
                                ttCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCCCTACTACAGGATTATAGCTGCCCCTGTAGGAGCGGGTTTACCcgcg                                                                                                                       >  8:319158/1‑146 (MQ=255)
                                 tCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTGTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTTTGGCTGCCCCGGTAGGAGCGGGTTTTCCCGCGaa                                                                                                                     >  1:276806/1‑147 (MQ=38)
                                     ttGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGGATTATCGCTGACCCTGTAGGAGCGGGTTTTCCCGCGAAGGGGc                                                                                                                >  8:325942/1‑148 (MQ=255)
                                     ttGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                     >  6:151730/1‑143 (MQ=255)
                                     ttGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                     <  5:151730/143‑1 (MQ=255)
                                       gTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCa                                                                                                              <  2:261981/148‑1 (MQ=255)
                                          aGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCCGCTGCCCCTGTAGGAGCGGGTTTACCcgcg                                                                                                                       <  1:212791/137‑1 (MQ=255)
                                          aGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTAAAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCcgcg                                                                                                                       >  2:212791/1‑136 (MQ=255)
                                           ggCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTTTCGGTGCCCCGGTTGGGGCGGGTTTTCCCGCGAAGAGGCCAGGAcc                                                                                                         >  3:93264/1‑149 (MQ=38)
                                              aaTGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCgg                                                                                                       >  1:15818/1‑148 (MQ=255)
                                               aTGAAAACCCGGTACCCTTGGGGTTCCCGTTTTTTTTTTTTCTTTATGTGTGTTGGCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCGCATTCTCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                <  6:142816/138‑1 (MQ=25)
                                               aTGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTACTACAGGAATATAGCTGCCACTGTAGGGGCGGGTTTAACCGGGAAGAGGCCatatccg                                                                                                        >  5:142816/1‑140 (MQ=38)
                                               aTGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCGCCTACAGGATTATCGCTGGCCCGGTAGGAGCGGGTTTTGCCGCGAAGAGGCCAGCCCCGGc                                                                                                      >  3:104637/1‑148 (MQ=38)
                                                     aaaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTTAAAACGCTCCTAGAGGATTATAGGTGCCCCTGTAGGAGCCGGTTTACCCCCCCAGAGGCCAGCACCGGCAAcc                                                                                                  >  5:99770/1‑146 (MQ=38)
                                                      aaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACCAGGGGCAGCGATAATCCTGTAGGAGCggg                                                                                                                                 <  1:26119/115‑1 (MQ=21)
                                                      aaCGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCggg                                                                                                                                 >  2:26119/1‑114 (MQ=21)
                                                          ggAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAAcc                                                                                          >  2:91108/1‑149 (MQ=255)
                                                             aCCTTGGGGTTCCGGTTTTTTTTTATTCTGATTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTa                                                                                       <  8:258941/149‑1 (MQ=255)
                                                             aCCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTa                                                                                       <  7:325942/149‑1 (MQ=255)
                                                                   gggTTCCGTTTTTTTTTTATTCTTAATGTGGGTTGCCAGTGCCGGCCCTTTCGGGGGTAACCCCGCTCCTACAGCATTTTCGATCCCCCTGTTGTTGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTgc                                                                                 <  6:99770/149‑1 (MQ=25)
                                                                    ggTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTCCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTgcg                                                                                <  2:15818/149‑1 (MQ=255)
                                                                         ccGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCGCACCTAAAGGATTATCGCTGCCCCTGTTGGAGCGGGTTTTCCCGCGAAGAGGCCCGCACCGGCAACCGCGAAACCAAAAAGTGCGCCCAg                                                                           >  6:140707/1‑149 (MQ=37)
                                                                           gTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAAACCCCTCCTAAAGGATTATAGCTGCCCCTGTAGGAGCGGGTTTAACCGCGAAGAGGCCAGCACCGGCCACCGAGAAACCCTACTGTGCGCACAGcc                                                                         >  3:318659/1‑149 (MQ=39)
                                                                           gTTTTTTTTTATTCTCAATGGGTGTGGTCCGTGCCGCCCCTTTCGCGGGAAACCCCGCCCCTACCGGATTATCGCTCCCCCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGcc                                                                         <  4:93264/149‑1 (MQ=37)
                                                                               ttttttATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCTACCGAGAAACCCTACTGTGCGCACAGCCCTa                                                                      <  1:213838/148‑1 (MQ=255)
                                                                                  tatATTCTTAATGTGTGTTGTCTGGGCGGGCCCTTTCGCGGTTAAACCCGCTCATACAGGATTATCGCTGCCCCTTTAGGAGGGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGt                                                                  <  1:309074/147‑1 (MQ=255)
                                                                                          gAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCGCCGGTAGGAGCGGGGTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTctg                                                          >  7:224047/1‑149 (MQ=255)
                                                                                           aaTGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTTTTGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACCCCCCAACTGTTCATTctgc                                                         >  2:80400/1‑149 (MQ=38)
                                                                                           aaTGTGTGTTGTCAGTGCCGGCCCTTTCGCAGGTAAACCCCCACCTACAGGATTATAGCTGCCCCAGTAGGAGCCGGTTTAACCGCGAAGAGGCCCGCACCGGCAACCCACAAACCCAACAGTGCGCACCGCCCGACAGTTCATTctgc                                                         >  4:168010/1‑149 (MQ=255)
                                                                                            ttttgtgtTGTATGTGCCGGCCCTTTCGCGTGTAACCCCTCTCCTACATTATTATCGCTGCCCCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTctgc                                                         <  5:213099/145‑1 (MQ=38)
                                                                                                gtgtTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAc                                                                                                                            <  5:213896/77‑1 (MQ=255)
                                                                                                gtgtTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAc                                                                                                                            >  6:213896/1‑77 (MQ=255)
                                                                                                   ttgtCAGTGCCGGCCCTTTCGCGGGTAAACCCCCTCCTACAGGATTATAGCTGCCCCTGTAGGAGCCGGTTTACCCGCGAAGAGGCCAGCACCCGCAACCGAGAAACCCTACAGTGCGCACAGCCCTACAGTTCATTATGCTGCTACt                                                  >  2:73039/1‑148 (MQ=38)
                                                                                                       ccgTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCg                                             <  4:318659/147‑1 (MQ=255)
                                                                                                         gTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCAACTCTGCGCACAGCCCAACTGTTCATTCTGCTGCTACTAACCGc                                            >  8:271611/1‑148 (MQ=255)
                                                                                                          tGCCGGCCCTTTCGCGGGTAAACCGGCTCCTACTGGATTATCGCTGCCCCTGTAGGAGGGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCt                                          <  7:271611/149‑1 (MQ=255)
                                                                                                              ggCCCTTTCGCGGGTAAAACCCCTACAACAGGATTATACCATCCCCAGTAGGAACCGGTTTAACCCCCAACAGGCCACCACCGGCAACCCACAAACCCAACAGTGCGCACACCCCAAAAGTTCATTATGATGGTACTCACCGCCAGCaa                                      >  6:189237/1‑149 (MQ=255)
                                                                                                               gCCATTTCGCGGTTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAAc                                     <  7:27727/149‑1 (MQ=255)
                                                                                                                       gcgGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                    >  3:127269/1‑110 (MQ=255)
                                                                                                                       gcgGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                    <  4:127269/110‑1 (MQ=255)
                                                                                                                       gcgGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGccc                             >  4:294619/1‑149 (MQ=255)
                                                                                                                             aaaCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCga                                                                                                <  5:204892/76‑1 (MQ=255)
                                                                                                                             aaaCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCga                                                                                                >  6:204892/1‑76 (MQ=255)
                                                                                                                              aaCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                      >  7:135918/1‑149 (MQ=255)
                                                                                                                                 ccGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAg                   >  5:258136/1‑149 (MQ=37)
                                                                                                                                 ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                      >  6:264374/1‑146 (MQ=255)
                                                                                                                                 ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                      <  5:264374/146‑1 (MQ=255)
                                                                                                                                         aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg           >  2:159448/1‑149 (MQ=255)
                                                                                                                                         aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg           >  8:142231/1‑149 (MQ=255)
                                                                                                                                         aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg           >  7:38934/1‑149 (MQ=255)
                                                                                                                                           cAGGATTTACGCTGCCCCTGTAGGAGCTGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg            <  5:189237/147‑1 (MQ=255)
                                                                                                                                           cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                  >  7:228471/1‑125 (MQ=255)
                                                                                                                                           cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                  <  8:228471/125‑1 (MQ=255)
                                                                                                                                           cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg           <  8:38934/148‑1 (MQ=255)
                                                                                                                                            aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGcccc                            <  6:302229/130‑1 (MQ=255)
                                                                                                                                            aGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGcccc                            >  5:302229/1‑130 (MQ=255)
                                                                                                                                             ggATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtgg        <  1:91108/149‑1 (MQ=255)
                                                                                                                                              gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                    >  5:20812/1‑136 (MQ=255)
                                                                                                                                              gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa                    <  6:20812/136‑1 (MQ=255)
                                                                                                                                              gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtg         <  5:288463/147‑1 (MQ=255)
                                                                                                                                              gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtg         >  6:288463/1‑147 (MQ=255)
                                                                                                                                              gATTATCGCTGCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtggt       <  3:294619/149‑1 (MQ=255)
                                                                                                                                                  aTCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGa   >  4:38520/1‑149 (MQ=255)
                                                                                                                                                   tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAc  <  8:135918/149‑1 (MQ=255)
                                                                                                                                                   tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAc  >  8:246213/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAC  >  NC_002947/2456718‑2457012

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: