Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,497,921 C→G 7.2% L332V (CTG→GTG)  cobL → precorrin‑6Y C(5,15)‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,497,9210CG 7.2% 205.8 / 4.9 67L332V (CTG→GTG) cobLprecorrin‑6Y C(5,15)‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base G (2/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GCCGCCAGGGTTTGATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGG  >  NC_002947/5497775‑5498069
                                                                                                                                                    |                                                                                                                                                       
gccgccAGGGTTTGATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTg                                                                                                                                                       <  4:125186/149‑1 (MQ=255)
     cAGGGTTTGATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCt                                                                                                                                                  <  6:246431/149‑1 (MQ=255)
        ggTTTGATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCtgtg                                                                                                                                               <  6:58994/149‑1 (MQ=255)
         gTTTGATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTCCAGTACACCTGTGc                                                                                                                                              <  5:134525/149‑1 (MQ=255)
          tttGATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCt                                                                                                                                             >  5:125452/1‑149 (MQ=255)
            tGATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTCCAGTACACCTGTGcgg                                                                                                                                            <  7:42339/148‑3 (MQ=255)
             gATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑gg                                                                                                                                          <  1:6560/148‑1 (MQ=255)
               tCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGc                                                                                                                                        >  4:348904/1‑148 (MQ=255)
                 gAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAg                                                                                                                                     >  4:39835/1‑149 (MQ=255)
                  aGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCt                                                                                                                                                  >  8:250283/1‑136 (MQ=255)
                  aGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCt                                                                                                                                                  <  7:250283/136‑1 (MQ=255)
                  aGCACAATCGCGATGCACTGGGTGTACCCGCCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGc                                                                                                                                    >  1:32468/1‑149 (MQ=255)
                       aaTCGCGCTGCACTGGGTGTAACCGGCCTGCAACTGGTTCGCGGCAATGCGCCAGCGGCGGTGCAATGCC‑TGGAACAACCGGCTGCGATCTTTAATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTgc                                                                                                                                <  1:104480/149‑1 (MQ=255)
                         tCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCCCGCCTGCGcc                                                                                                                             >  5:207182/1‑149 (MQ=255)
                         tCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGcc                                                                                                                             >  8:313311/1‑149 (MQ=255)
                                 cACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCggcgg‑c                                                                                                                    <  7:90229/149‑1 (MQ=255)
                                  aCTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGg‑cc                                                                                                                   >  5:269783/1‑149 (MQ=255)
                                  aCTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTCCAGTACACCTGTGcg                                                                                                                                             <  3:212867/125‑2 (MQ=255)
                                  aCTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTCCAGTACACCTGTGcg                                                                                                                                             >  4:212867/1‑124 (MQ=255)
                                      ggTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGg‑ccg‑cc                                                                                                               <  7:313311/148‑1 (MQ=255)
                                      ggTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCt                                                                                                              >  8:258535/1‑149 (MQ=255)
                                      ctggtACCCGGACTGCAACTTGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCt                                                                                                              <  4:190163/146‑1 (MQ=255)
                                        tgtACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCtgg                                                                                                            >  1:199030/1‑149 (MQ=255)
                                            cccGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGc                                                                                                        >  8:131057/1‑149 (MQ=255)
                                             ccGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCCGGGGACCACCCGGATGGGAGC‑TTCATCGGGGGCGCCGTTCCCCGCACAGGTGTAGGGCCCCTGTGGT‑GGGCAAGCCTCCGCCCCGGCGG‑CCG‑CCtggtgcc                                                                                                        >  1:278074/1‑147 (MQ=255)
                                                     cAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCA‑CCAGGCGGccgccggg                                                                                                 <  6:125452/147‑7 (MQ=255)
                                                           gTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCt                                                                                         <  2:52921/149‑1 (MQ=255)
                                                           gTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCt                                                                                         >  8:12651/1‑149 (MQ=255)
                                                                cgGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCaaa                                                                                    <  3:39835/149‑1 (MQ=255)
                                                                cgGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCaaa                                                                                    <  4:204671/149‑1 (MQ=255)
                                                                cgGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCa                                                                                      >  3:204671/1‑147 (MQ=255)
                                                                 ggCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGcc                                                                                                 <  5:273425/135‑1 (MQ=255)
                                                                 ggCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGcc                                                                                                 >  6:273425/1‑135 (MQ=255)
                                                                  gCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATTCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGTTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGc                                                                                  <  4:80518/149‑1 (MQ=255)
                                                                    aaGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCg                                                                                 <  4:112485/148‑1 (MQ=255)
                                                                    aaGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCg                                                                                 >  3:112485/1‑148 (MQ=255)
                                                                    aaGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGa                                                                                <  6:269783/149‑1 (MQ=255)
                                                                       gcgcCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACt                                                                             >  2:98281/1‑149 (MQ=255)
                                                                          ccAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGc                                                                          >  1:205847/1‑149 (MQ=255)
                                                                           cAGCGGCGCTGCAAGGCC‑TCGAACAACCGGATGCGATC‑TTCATCGGAGCCGGCGTTACCCGCACAGGTGTACTGGACCTCTGCT‑GGGCAAGCCAGCGCCCCCGCGGACAGTGCAGGTCGCCAACGCCGTCACCCAGCAAAGCCAACTgg                                                                           >  1:316275/1‑149 (MQ=255)
                                                                            aGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGccc                                                                        <  5:138548/149‑1 (MQ=255)
                                                                             gcggcgCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCt                                                                       >  8:157414/1‑149 (MQ=255)
                                                                                  gcTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGcac                                                                  <  8:356657/149‑1 (MQ=255)
                                                                                   cTGCAAGGCC‑TGTAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACTGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGcaca                                                                 <  2:48304/149‑1 (MQ=255)
                                                                                      cAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGc                                                                                  <  5:366024/129‑1 (MQ=255)
                                                                                      cAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGc                                                                                  >  6:366024/1‑129 (MQ=255)
                                                                                            c‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTcgcg                                                         >  1:352363/1‑148 (MQ=255)
                                                                                            c‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTcgcg                                                         >  7:138898/1‑148 (MQ=255)
                                                                                            c‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGa                                                        <  7:211822/149‑1 (MQ=255)
                                                                                                gAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGaa                                                       <  7:131057/147‑1 (MQ=255)
                                                                                                gAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACa                                                     <  2:32468/149‑1 (MQ=255)
                                                                                                 aacaacCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCcggcg                                                                                                                       <  8:299777/84‑1 (MQ=255)
                                                                                                 aacaacCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCcggcg                                                                                                                       >  7:299777/1‑84 (MQ=255)
                                                                                                      ccGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACgg                                               >  2:296865/1‑149 (MQ=255)
                                                                                                      ccGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGCCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACgg                                               >  5:327701/1‑149 (MQ=255)
                                                                                                      ccGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGCACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCT‑CCTGGAGGCCAACGCCGTAACCCTGCAAACCGAACTGCCCCTGGCACATTTTCGCGAACAGCACgg                                               >  7:59166/1‑149 (MQ=255)
                                                                                                            gCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACggtggt                                           >  7:317338/1‑147 (MQ=255)
                                                                                                             ctaTC‑TTCATCGGATGCTGCGTTACCCGCACAGTTGTACTTGACCTGTGCT‑GGGCAAGCATGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAg                                        <  8:166378/147‑1 (MQ=255)
                                                                                                             cGATC‑TTCCTCGGAGGCGGCGTTACCCGCACACGTGTACTTGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAg                                        <  6:207182/149‑1 (MQ=255)
                                                                                                                     cATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGGGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGa                                    <  6:21964/146‑1 (MQ=255)
                                                                                                                       tCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCTGGGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCg                                >  1:263534/1‑149 (MQ=255)
                                                                                                                          gAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATc                            >  4:283741/1‑149 (MQ=255)
                                                                                                                             gcggcgTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCAt                         <  2:205847/149‑1 (MQ=255)
                                                                                                                                     aCCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCa                 >  4:156631/1‑149 (MQ=255)
                                                                                                                                       ccGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCT‑GCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCAt                             >  1:212351/1‑134 (MQ=255)
                                                                                                                                       ccGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCT‑GCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCAt                             <  2:212351/134‑1 (MQ=255)
                                                                                                                                               gtgtACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGccc       <  2:223462/149‑1 (MQ=255)
                                                                                                                                                    cTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGGACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTgg   <  2:263534/148‑1 (MQ=255)
                                                                                                                                                    cTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTgg   <  7:157414/148‑1 (MQ=255)
                                                                                                                                                    cTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTgg   <  2:199030/148‑1 (MQ=255)
                                                                                                                                                    cTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTgg   <  3:348904/148‑1 (MQ=255)
                                                                                                                                                    cTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTgg   <  7:12651/148‑1 (MQ=255)
                                                                                                                                                    cTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑ACTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTgg   <  8:59166/148‑1 (MQ=255)
                                                                                                                                                    cTGGACCTGTGCT‑GGGCAAGCCTGCGCACCGGCGG‑CCG‑CCTGGTGGCCAACGCCGACACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTggg  >  6:323928/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                       
GCCGCCAGGGTTTGATCGAGCACAATCGCGATGCACTGGGTGTACCCGGCCTGCAACTGGTTCGCGGCAAGGCGCCAGCGGCGCTGCAAGGCC‑TGGAACAACCGGATGCGATC‑TTCATCGGAGGCGGCGTTACCCGCACAGGTGTACTGGACCTGTGCT‑GGGCAAGCCTGCGCCCCGGCGG‑CCG‑CCTGGTGGCCAACGCCGTCACCCTGCAAAGCGAACTGGCCCTGGCACATTTTCGCGAACAGCACGGTGGTGAGCTGACCCGCATCCATGTCGCGCATGCCCAGCCCTTGGG  >  NC_002947/5497775‑5498069

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: ‑