Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,871 C→T 6.6% intergenic (+122/‑59) 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,8710CT 6.6% 223.0 / 3.9 77intergenic (+122/‑59)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base T (3/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

TGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGC  >  NC_002947/2456724‑2457019
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                                         <  4:188059/148‑1 (MQ=37)
             gTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGGGCAGCGATAATCCTGTAGGAGCgg                                                                                                                                         >  3:188059/1‑148 (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                                                                           <  4:127269/110‑1 (MQ=255)
                                                                                                                 gcgGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACt                                                                           >  3:127269/1‑110 (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                             >  6:264374/1‑146 (MQ=255)
                                                                                                                           ccGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTg                             <  5:264374/146‑1 (MQ=255)
                                                                                                                                   aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                  >  7:38934/1‑149 (MQ=255)
                                                                                                                                   aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                  >  8:142231/1‑149 (MQ=255)
                                                                                                                                   aCAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTggg                  >  2:159448/1‑149 (MQ=255)
                                                                                                                                     cAGGATTTACGCTGCCCCTGTAGGAGCTGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTgg                   <  5:189237/147‑1 (MQ=255)
                                                                                                                                     cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                         <  8:228471/125‑1 (MQ=255)
                                                                                                                                     cAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCaa                                         >  7:228471/1‑125 (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                >  6:288463/1‑147 (MQ=255)
                                                                                                                                        gATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCAGGCGCTCCAGggtg                <  5:288463/147‑1 (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)
                                                                                                                                                 cTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTa                                                                                              >  8:17110/1‑61 (MQ=255)
                                                                                                                                                 cTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTa                                                                                              <  7:17110/61‑1 (MQ=255)
                                                                                                                                                    ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg   <  3:38520/149‑1 (MQ=255)
                                                                                                                                                    ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg   <  1:159448/149‑1 (MQ=255)
                                                                                                                                                    ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg   >  6:87516/1‑149 (MQ=255)
                                                                                                                                                    ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTACAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATg     <  5:140707/147‑1 (MQ=255)
                                                                                                                                                     cccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGgc  <  5:87516/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGGC  >  NC_002947/2456724‑2457019

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: