Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,862 T→C 6.5% intergenic (+113/‑68) 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,8620TC 6.5% 215.0 / 4.1 77intergenic (+113/‑68)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (36/36);  new base C (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 = 6.08e-01

TGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGT  >  NC_002947/2456715‑2457007
                                                                                                                                                    |                                                                                                                                                  
tGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTa                                                                                                                                                   >  1:91600/1‑149 (MQ=255)
   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                                                                                                                             <  1:154487/137‑1 (MQ=255)
                                 gATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCgg                                                                                                                             >  2:154487/1‑137 (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                                                                                                                <  5:151730/143‑1 (MQ=255)
                                        ttGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                >  6:151730/1‑143 (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                                                                                           >  6:204892/1‑76 (MQ=255)
                                                                                                                                aaaCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCga                                                                                           <  5:204892/76‑1 (MQ=255)
                                                                                                                                 aaCCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                 >  7:135918/1‑149 (MQ=255)
                                                                                                                                    ccGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAg              >  5:258136/1‑149 (MQ=37)
                                                                                                                                    ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                 <  5:264374/146‑1 (MQ=255)
                                                                                                                                    ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                 >  6:264374/1‑146 (MQ=255)
                                                                                                                                            aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg      >  8:142231/1‑149 (MQ=255)
                                                                                                                                            aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg      >  7:38934/1‑149 (MQ=255)
                                                                                                                                            aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg      >  2:159448/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               <  6:20812/136‑1 (MQ=255)
                                                                                                                                                 gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGa               >  5:20812/1‑136 (MQ=255)
                                                                                                                                                 gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtg    >  6:288463/1‑147 (MQ=255)
                                                                                                                                                 gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtg    <  5:288463/147‑1 (MQ=255)
                                                                                                                                                 gATTATCGCTGCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGggtggt  <  3:294619/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGT  >  NC_002947/2456715‑2457007

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: