Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,772,156 G→C 6.0% intergenic (+75/‑232) PP_3331 → / → PP_3332 membrane protein/cytochrome c‑type protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,772,1560GC 6.0% 242.0 / 4.0 78intergenic (+75/‑232)PP_3331/PP_3332membrane protein/cytochrome c‑type protein
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base C (3/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

TGGTGGTTTCACTCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTAGATGAAGCCCT  >  NC_002947/3772008‑3772297
                                                                                                                                                        |                                                                                                                                             
tggtggTTTCACTCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATg                                                                                                                                               >  7:7406/1‑149 (MQ=255)
  gtggtTTCACTCAGCGGCGTCGTACTAGGCTGGAGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGATTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATgag                                                                                                                                             >  2:263177/1‑147 (MQ=255)
      tttCACTCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGcccc                                                                                                                                         <  6:43598/149‑1 (MQ=255)
       ttCACTCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGTGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCg                                                                                                                                        >  1:72946/1‑149 (MQ=255)
        tCACTCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCg                                                                                                                                        >  1:314546/1‑149 (MQ=255)
         cACTCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCgg                                                                                                                                       <  4:8059/148‑1 (MQ=255)
          aCTCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCt                                                                                                                                     <  2:206653/149‑1 (MQ=255)
            tCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGAAGGTCTCAAAAAACAAAAG‑‑CCGGGACGTTAATCCGCCCCGGAtt                                                                                                                                    >  1:132205/1‑149 (MQ=255)
            tCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCttt                                                                                                                                   >  6:37950/1‑149 (MQ=255)
              aGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCtttt                                                                                                                                  <  8:144741/148‑1 (MQ=255)
              aGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCtttt                                                                                                                                  <  4:277159/148‑1 (MQ=255)
              aGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCtttt                                                                                                                                  <  4:9873/148‑1 (MQ=255)
              aGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCtt                                                                                                                                    >  5:227518/1‑146 (MQ=255)
              aGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCtt                                                                                                                                    <  6:227518/146‑1 (MQ=255)
              aGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTtttg                                                                                                                                 >  4:153664/1‑149 (MQ=255)
               gcggcgTCGTACTAGTCTGGAGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCATCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTtttgt                                                                                                                                <  3:64069/149‑1 (MQ=255)
                ctgcgTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTtttgtt                                                                                                                               <  8:255455/147‑1 (MQ=255)
                cggcgTCGTACTTGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCTCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCCGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTtttgtt                                                                                                                               <  8:255453/149‑1 (MQ=255)
                cggcgTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCGCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGTGACGGTAATGCGCCCCGGCTtttgtt                                                                                                                               <  5:30643/149‑1 (MQ=255)
                    gtcgtACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAg                                                                                                                           <  2:369680/149‑1 (MQ=255)
                       gtACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGGCGCATTACCGTCCCGGCTtttgttttttca                                                                                                                         <  1:230568/148‑6 (MQ=37)
                        tACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGTGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGa                                                                                                                       <  2:72946/149‑1 (MQ=255)
                        tACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CAGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGa                                                                                                                       >  1:160071/1‑149 (MQ=255)
                         aCTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTCTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGAt                                                                                                                      >  7:196487/1‑149 (MQ=255)
                          cTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGAt                                                                                                                      <  5:37950/148‑1 (MQ=255)
                              gCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGaaaa                                                                                                                 >  2:111704/1‑149 (MQ=255)
                                  gCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTtttgtt                                                                                                                               <  4:93802/131‑1 (MQ=255)
                                  gCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTtttgtt                                                                                                                               >  3:93802/1‑131 (MQ=255)
                                  gCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAAccc                                                                                                              <  4:20831/148‑1 (MQ=255)
                                   cGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTg                                                                                                            >  6:326970/1‑149 (MQ=255)
                                   cGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATACGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTg                                                                                                            <  2:376772/149‑1 (MQ=255)
                                    gAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACTAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAg                                                                                                                           <  1:179749/133‑1 (MQ=255)
                                    gAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACTAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAg                                                                                                                           >  2:179749/1‑133 (MQ=255)
                                     aaGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGc                                                                                                           >  5:211375/1‑148 (MQ=255)
                                          tACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAagag                                                                                                     >  6:260396/1‑149 (MQ=255)
                                            cGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGt                                                                                                    >  3:41962/1‑148 (MQ=255)
                                             gCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTaa                                                                                                  <  5:260396/149‑1 (MQ=255)
                                             gCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTaa                                                                                                  <  8:7406/149‑1 (MQ=255)
                                                 aaTGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAaca                                                                                              >  8:351511/1‑149 (MQ=255)
                                                  aTGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAGTTCCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTaaaa                                                                                                <  1:111704/148‑1 (MQ=255)
                                                      ccGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGa                                                                                         <  3:153664/149‑1 (MQ=255)
                                                       cGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGAt                                                                                        >  3:179644/1‑149 (MQ=255)
                                                               ttCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTc                                                                                 >  4:85366/1‑148 (MQ=255)
                                                                  gcgGTTGAGTCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGGCCCGGCTTTTGTTTGAGCGGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAAc                                                                              <  4:27614/148‑1 (MQ=255)
                                                                          aGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCAc                                                                       <  4:41962/148‑1 (MQ=255)
                                                                          aGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACg                                                                      >  5:167040/1‑149 (MQ=255)
                                                                           gCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACgg                                                                     >  1:230118/1‑149 (MQ=255)
                                                                             ctcctTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAGAACACCGGATCAGCTTCAACGTAGCACGGCa                                                                   >  6:182768/1‑149 (MQ=255)
                                                                              tcctTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCAc                                                                  >  1:383716/1‑149 (MQ=255)
                                                                                ctTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCAc                                                                  <  5:177296/147‑1 (MQ=255)
                                                                                ctTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCt                                                                <  5:326970/149‑1 (MQ=255)
                                                                                     gtgtTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAg                                                           >  5:105854/1‑149 (MQ=255)
                                                                                          tCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAgccgc                                                       <  4:179644/148‑1 (MQ=255)
                                                                                           cGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTg                                                               <  1:210945/139‑1 (MQ=255)
                                                                                           cGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTg                                                               >  2:210945/1‑139 (MQ=255)
                                                                                             cAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGTTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCTGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCAcc                                                   <  2:140694/149‑1 (MQ=255)
                                                                                                       aGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGt                                         <  6:211375/149‑1 (MQ=255)
                                                                                                         cGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGGTTTTGTTTGAGCCGATGAAAACCGTGCCAAGAGTAAAACACCGGATCCGCTTCAACGTAGCACGGCACATGCAAGCCGCCACCCTGCGACGGTTc                                       >  1:83776/1‑149 (MQ=255)
                                                                                                         cGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTc                                       <  3:85366/149‑1 (MQ=255)
                                                                                                         cGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTc                                       <  2:230118/149‑1 (MQ=255)
                                                                                                           aGTGCCGCTGATCACACAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGATATTGTTTCAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGAAACTGCAAGCAGCAACACTGCGACGGTTcgc                                     >  8:267421/1‑149 (MQ=255)
                                                                                                             tGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGATGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTcgccg                                   >  8:333099/1‑149 (MQ=255)
                                                                                                             tGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTcgccg                                   <  2:160071/149‑1 (MQ=255)
                                                                                                                         aCTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTgccag                        <  8:116121/149‑1 (MQ=255)
                                                                                                                          cTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTgccagc                       >  3:186847/1‑149 (MQ=255)
                                                                                                                            gAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTgccagc                       >  5:365204/1‑147 (MQ=255)
                                                                                                                             ctcaaaCAAAAG‑‑CCGGGGCGCATTACCGTCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCg                                              <  6:138576/120‑1 (MQ=25)
                                                                                                                             ctcaaaCAAAAG‑‑CCGGGGCGCATTACCGTCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCg                                              >  5:138576/4‑123 (MQ=25)
                                                                                                                               aaaaCAAACG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGc                                                                                    <  4:75537/83‑1 (MQ=255)
                                                                                                                               aaaaCAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGc                                                                                    >  3:75537/1‑83 (MQ=255)
                                                                                                                                aaaCAAAAG‑‑CCGGGGCGCATTACCGTCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGAcc                 >  3:252491/1‑149 (MQ=38)
                                                                                                                                    aaaaG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTa             <  8:196487/149‑1 (MQ=255)
                                                                                                                                       aG‑‑CCGGTACGGTAATGCGCCCCGGCTTTTGTTTTAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGTATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTAGAt          <  1:263177/149‑1 (MQ=255)
                                                                                                                                            cGGGATGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTAGATGaa       <  7:333099/149‑1 (MQ=255)
                                                                                                                                            cGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCCCCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGAGCGGTAGATg         >  4:248099/1‑147 (MQ=255)
                                                                                                                                              ggACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTAGa           <  3:203074/143‑1 (MQ=255)
                                                                                                                                              ggACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTAGa           >  4:203074/1‑143 (MQ=255)
                                                                                                                                                 cGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGGTCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTAGATGAAGCCCt  >  3:356045/1‑149 (MQ=255)
                                                                                                                                                 cGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCTCACCCTGCCAGCCAGACCGGTAGATGAAGCCCt  <  2:314546/149‑1 (MQ=255)
                                                                                                                                                 cGGTAATGCGCCCCGGATTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTAGATGAAGCCCt  <  2:383716/149‑1 (MQ=255)
                                                                                                                                                        |                                                                                                                                             
TGGTGGTTTCACTCAGCGGCGTCGTACTAGGCTGGCGAAGGCTACGCAGAATGACCGCTCGAGTTCGCGGTTGAGCCTCCTTTCGTGTTTCGCAGCGGCCTGAGCGAGTGCCGCTGATCACTGAAAAAACAAAAG‑‑CCGGGACGGTAATGCGCCCCGGCTTTTGTTTGAGCCGATGAAAACCCTGCCAAGAGTAAAACACCGGATCAGCTTCAACGTAGCACGGCACCTGCAAGCCGCCACCCTGCGACGGTTCGCCGCACCCTGCCAGCCAGACCGGTAGATGAAGCCCT  >  NC_002947/3772008‑3772297

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