Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,860 A→G 6.4% intergenic (+111/‑70) 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,8600AG 6.4% 206.0 / 4.6 78intergenic (+111/‑70)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (37/36);  new base G (3/2);  total (40/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 = 1.48e-01

ACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGT  >  NC_002947/2456712‑2457007
                                                                                                                                                     |                                                                                                                                                    
aCATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGa                                                                                                                                                      >  3:42629/1‑149 (MQ=255)
  aTGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGAt                                                                                                                                                     >  7:306672/1‑148 (MQ=255)
   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                                                                                                                             <  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                                                                                                                             >  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                                                                                                                <  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                 >  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                             >  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)
                                                                                                                                                     |                                                                                                                                                    
ACATGGCTGGTGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGT  >  NC_002947/2456712‑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: