Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,456,870 C→A 7.4% intergenic (+121/‑60) 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,8700CA 7.4% 220.5 / 6.5 77intergenic (+121/‑60)gapB/PP_2150glyceraldehyde‑3‑phosphate dehydrogenase/thiamine biosynthesis lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (36/34);  new base A (4/2);  total (41/36)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.74e-01

TGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGG  >  NC_002947/2456722‑2457018
                                                                                                                                                      |                                                                                                                                                    
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                                                                                                                           <  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                                                                                                                                  >  6:213896/1‑77 (MQ=255)
                                                                                            gtgtTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTAc                                                                                                                                  <  5:213896/77‑1 (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                                        <  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                          <  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)
                                                                                                                                              aTCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGa         >  4:38520/1‑149 (MQ=255)
                                                                                                                                               tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAc        <  8:135918/149‑1 (MQ=255)
                                                                                                                                               tCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGAc        >  8:246213/1‑149 (MQ=255)
                                                                                                                                                   cTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTa                                                                                             <  7:17110/61‑1 (MQ=255)
                                                                                                                                                   cTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTa                                                                                             >  8:17110/1‑61 (MQ=255)
                                                                                                                                                      ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg  <  1:159448/149‑1 (MQ=255)
                                                                                                                                                      ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg  <  3:38520/149‑1 (MQ=255)
                                                                                                                                                      ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATggg  >  6:87516/1‑149 (MQ=255)
                                                                                                                                                      ccccTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTACAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATg    <  5:140707/147‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
TGTGAACCCGCCAGCGTTTCCACGCTGATTCGCTTGTAAGGCAATGAAAAAACGGGAACCTTCGGGTTCCCGTTTTTTTTTATTCTGAATGTGTGTTGTCAGTGCCGGCCCTTTCGCGGGTAAACCCGCTCCTACAGGATTATCGCTGCCCCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCAGCACCGGCAACCGAGAAACCCTACTGTGCGCACAGCCCTACTGTTCATTCTGCTGCTACTCACCGCCTGCAACCAAGGCCCCACCCTGGAGCGCCTGGGTGGTCCGACCATGGG  >  NC_002947/2456722‑2457018

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: