Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,727 C→G 6.4% intergenic (+30/‑60) allE → / → PP_3531 S‑ureidoglycine aminohydrolase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,002,7270CG 6.4% 294.0 / 8.0 90intergenic (+30/‑60)allE/PP_3531S‑ureidoglycine aminohydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (47/36);  new base G (2/4);  total (50/40)
Fisher's exact test for biased strand distribution p-value = 4.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAG  >  NC_002947/4002583‑4002874
                                                                                                                                                |                                                                                                                                                    
gCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGc                                                                                                                                                 >  1:363511/1‑150 (MQ=255)
   gcgcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCggg                                                                                                                                              <  2:372896/150‑1 (MQ=255)
     gcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACACGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCt                                                                                                                                            >  2:144057/1‑150 (MQ=255)
         cTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAATGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCt                                                                                                                                        >  2:608704/1‑150 (MQ=255)
          ttCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTg                                                                                                                                       >  2:1127367/1‑150 (MQ=255)
          ttCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTg                                                                                                                                       >  2:848812/1‑150 (MQ=255)
           tCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGCCAACTTTAGCTACCGGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGACCCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGt                                                                                                                                      <  2:643284/150‑1 (MQ=255)
            cTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGtt                                                                                                                                     <  2:584568/150‑1 (MQ=255)
               cccGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCGTTCTGTTTGc                                                                                                                                  >  2:586395/1‑150 (MQ=255)
                  gCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGATGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCcgc                                                                                                                               >  1:351194/1‑150 (MQ=255)
                       ccTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATa                                                                                                                          <  2:329231/150‑1 (MQ=255)
                       ccTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATa                                                                                                                          >  1:1608015/1‑150 (MQ=255)
                         tGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGTGGGCTTTCTGTTTGCCGCGCa                                                                                                                            <  2:888642/146‑1 (MQ=255)
                         tGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGTGGGCTTTCTGTTTGCAGCGCa                                                                                                                            >  1:888642/1‑146 (MQ=255)
                            tACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGc                                                                                                                                             <  1:857045/126‑1 (MQ=255)
                             aCGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGaaa                                                                                                                    <  2:1354674/150‑1 (MQ=255)
                                   gcggcCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATaccca                                                                                                              >  1:138882/1‑150 (MQ=255)
                                       cccAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCt                                                                                                          >  1:1116598/1‑150 (MQ=255)
                                         cAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGATAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTg                                                                                                        >  2:875414/1‑150 (MQ=255)
                                           ggCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCa                                                                                                      <  2:363511/150‑1 (MQ=255)
                                           ggCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCa                                                                                                      <  2:138882/150‑1 (MQ=255)
                                             cAACTTTAGCTACCTGCTGTACAAGGAGCTTCGCCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAgg                                                                                                    >  1:41914/1‑150 (MQ=255)
                                                 tttAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTAc                                                                                                <  2:1627347/150‑1 (MQ=255)
                                                         ccTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTgg                                                                                                                                                  <  2:573004/92‑1 (MQ=255)
                                                         ccTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTgg                                                                                                                                                  >  1:573004/1‑92 (MQ=255)
                                                           tgctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGaa                                                                                      >  1:1495508/1‑150 (MQ=255)
                                                             ctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGAACCggttt                                                                                          >  1:522552/1‑141 (MQ=255)
                                                             ctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGAACCggttt                                                                                          <  2:522552/144‑4 (MQ=255)
                                                               gTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGc                                                                                 >  1:295863/1‑150 (MQ=255)
                                                                  cAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgc                                                                              <  1:881533/150‑1 (MQ=255)
                                                                  cAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAAGACCCACCCTTGCAGGGTACCGGCCATGAACCGctgc                                                                              <  2:1446018/150‑1 (MQ=255)
                                                                   aaGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcc                                                                             >  1:839487/1‑150 (MQ=255)
                                                                   aaGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcc                                                                             >  1:144966/1‑150 (MQ=255)
                                                                        cGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGTGCTAAGCCCNGGCTGATGCCCGCAAGCCCGCCAAAATGGCTGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGCCCCTCTGCctgctg                                                                        >  1:108859/1‑150 (MQ=255)
                                                                          tGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTgctgct                                                                      >  2:964080/1‑150 (MQ=255)
                                                                           gAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTc                                                                     <  1:144057/150‑1 (MQ=255)
                                                                              ccGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgcgg                                                                                                                                       <  2:104544/82‑5 (MQ=37)
                                                                              ccGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgcgg                                                                                                                                       >  1:104544/1‑78 (MQ=37)
                                                                               cGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATac                                                                                                                 <  1:429349/103‑1 (MQ=255)
                                                                               cGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATac                                                                                                                 >  2:429349/1‑103 (MQ=255)
                                                                               cGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATac                                                                                                                 >  1:227759/1‑103 (MQ=255)
                                                                               cGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATac                                                                                                                 <  2:227759/103‑1 (MQ=255)
                                                                                     gTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATTAGCCGCTGCCTGCTGCTCGCTGCGTtgg                                                           >  1:343039/1‑149 (MQ=255)
                                                                                        cacCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcctgc                                                                          >  2:1531869/1‑132 (MQ=255)
                                                                                        cacCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcctgc                                                                          <  1:1531869/132‑1 (MQ=255)
                                                                                         acCTGACGCTGAACCCCAAACGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGCGTACCGACCATGAACCGCTGCCTGCTGCTCGCTGCCTTcctgct                                                       >  2:212237/1‑150 (MQ=255)
                                                                                           cTGACGTTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                     >  2:1228803/1‑150 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATa                                                                                                                          >  1:886323/1‑82 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATa                                                                                                                          <  2:886323/82‑1 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATaccca                                                                                                              >  1:361998/1‑94 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATaccca                                                                                                              <  2:361998/94‑1 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCg                                                                                              <  1:1549804/110‑1 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGCAc                                                                                    <  2:1258555/119‑1 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAAc                                                                                    >  1:1258555/1‑119 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                     <  1:1127367/150‑1 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                     <  1:964080/150‑1 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTGGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCg                                                                                              >  2:1549804/1‑110 (MQ=255)
                                                                                            tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgtg                                                    >  1:1344051/1‑150 (MQ=255)
                                                                                            tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgtg                                                    <  2:1608015/150‑1 (MQ=255)
                                                                                                gCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCGCTGCTCGCTGCCTTGCTGCTGTGCACct                                               >  2:1366783/1‑150 (MQ=255)
                                                                                                   gAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACctct                                             <  1:1228803/150‑1 (MQ=255)
                                                                                                     aCCCCAAGGGCTGAGCCCAGGCTGATGGCCGCAAGCCCGCCAGGCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGc                                           <  1:586395/150‑1 (MQ=255)
                                                                                                      ccccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGcc                                          >  2:1020488/1‑150 (MQ=255)
                                                                                                       cccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttg                                                                                                                                          <  1:1384531/54‑2 (MQ=21)
                                                                                                       cccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttg                                                                                                                                          >  2:1384531/1‑53 (MQ=21)
                                                                                                       cccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGccc                                         >  2:1642165/1‑150 (MQ=255)
                                                                                                       cccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTACCGGCGACTCGCTGCCTTGCTTCTGTGCACCTCTGccc                                         >  2:1319293/1‑150 (MQ=255)
                                                                                                         cAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGctgctg                                                      >  2:452961/1‑135 (MQ=255)
                                                                                                         cAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGctgctg                                                      <  1:452961/135‑1 (MQ=255)
                                                                                                         cAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTg                                       <  1:875414/150‑1 (MQ=255)
                                                                                                          aaGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGt                                                                                                                                      >  1:687262/1‑55 (MQ=255)
                                                                                                          aaGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGt                                                                                                                                      <  2:687262/55‑1 (MQ=255)
                                                                                                              gcTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGgccgc                                  <  2:839487/150‑1 (MQ=255)
                                                                                                                 aaGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGa                               <  1:212237/150‑1 (MQ=255)
                                                                                                                  aGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGAt                              <  2:1495508/150‑1 (MQ=255)
                                                                                                                   acccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATg                             >  2:33414/2‑150 (MQ=255)
                                                                                                                    cccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGa                            >  1:1530314/1‑150 (MQ=255)
                                                                                                                     ccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCt                                                                  <  1:572738/111‑1 (MQ=255)
                                                                                                                     ccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCt                                                                  >  2:572738/1‑111 (MQ=255)
                                                                                                                     ccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGAc                           >  1:477775/1‑150 (MQ=255)
                                                                                                                      cAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACt                          <  2:1530314/150‑1 (MQ=255)
                                                                                                                              tGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTg                  >  1:362847/1‑150 (MQ=255)
                                                                                                                                 ccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGcc               >  2:456053/1‑150 (MQ=255)
                                                                                                                                    cAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGATCCCTGCCTTGCGGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTAc            <  1:1020488/150‑1 (MQ=255)
                                                                                                                                      aGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGCACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGa          <  2:327045/150‑1 (MQ=255)
                                                                                                                                       gCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGAc         >  1:664117/1‑150 (MQ=255)
                                                                                                                                         ccGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACcg       >  1:783724/1‑150 (MQ=255)
                                                                                                                                           gCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACcgcg     >  1:157309/1‑150 (MQ=255)
                                                                                                                                           gCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCAACTATGCCATGGCCGCCGATGACTGGCGCCTGGCCTACGACcgcg     >  1:506331/1‑150 (MQ=255)
                                                                                                                                              aaCCTGGCGGGCTTTCTGTTTGCCTCGCATACGTAAATACCCACCCTTGCAGGGTACCGGCCATGAAGGCGCGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAg  <  1:1179428/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAG  >  NC_002947/4002583‑4002874

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