Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,682 A→G 8.5% intergenic (+29/‑61) 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,6820AG 8.5% 195.0 / 7.4 72intergenic (+29/‑61)allE/PP_3531S‑ureidoglycine aminohydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (39/25);  new base G (3/3);  total (43/29)
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.65e-01

TGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCG  >  NC_002947/4002541‑4002827
                                                                                                                                             |                                                                                                                                                 
tGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCg                                                                                                                                            >  6:32357/1‑149 (MQ=255)
  cgcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCg                                                                                                                                            <  4:247914/147‑1 (MQ=255)
    cgcCTTCTGCCCGCAAGCCGGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCt                                                                                                                                        <  8:66209/149‑1 (MQ=255)
       cTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTc                                                                                                                                     >  6:311473/1‑149 (MQ=255)
        ttCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCttt                                                                                                                                      >  5:330771/1‑147 (MQ=255)
        ttCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCttt                                                                                                                                      <  6:330771/147‑1 (MQ=255)
          cTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCt                                                                                                                                    <  2:137395/147‑1 (MQ=255)
           tGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGtt                                                                                                                                 >  6:355904/1‑149 (MQ=255)
             cccGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGAAAGCCCGCCAGGTTGGCGGGCttgcggg                                                                                                                                  <  7:210289/146‑6 (MQ=255)
              ccGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACCCCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCGTTCTGTTTGc                                                                                                                              >  7:311421/1‑149 (MQ=255)
               cGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCtt                                                                                                                                       >  5:243044/1‑139 (MQ=255)
               cGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCtt                                                                                                                                       <  6:243044/139‑1 (MQ=255)
                  aaGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCAGCCACCCTGGCGGGCTTTCTGTTTGCcgcg                                                                                                                          >  8:257781/1‑149 (MQ=255)
                      cTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATa                                                                                                                      >  6:28441/1‑149 (MQ=255)
                       tGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGttt                                                                                                                                >  8:91163/1‑138 (MQ=255)
                       tGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGttt                                                                                                                                <  7:91163/138‑1 (MQ=255)
                         cTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACt                                                                                                                    <  5:32357/148‑1 (MQ=255)
                         cTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACt                                                                                                                    >  3:96931/1‑148 (MQ=255)
                         cTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGAAAGCCCGCCAGGTTGGCGGGc                                                                                                                                         >  2:67990/1‑127 (MQ=255)
                         cTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGAAAGCCCGCCAGGTTGGCGGGc                                                                                                                                         <  1:67990/127‑1 (MQ=255)
                                        aGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTg                                                                                                    <  4:356801/149‑1 (MQ=255)
                                            aaCTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCACGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAgg                                                                                                <  4:196175/149‑1 (MQ=255)
                                             aCTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAggg                                                                                               >  5:197422/1‑149 (MQ=255)
                                                   gCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCgg                                                                                         <  4:212197/149‑1 (MQ=255)
                                                    cTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGc                                                                                        <  5:311473/149‑1 (MQ=255)
                                                    cTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGACCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGc                                                                                        >  5:716/1‑149 (MQ=255)
                                                    cTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATCCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCAAGCAGGGTACCgga                                                                                        >  5:373975/1‑148 (MQ=255)
                                                      aCCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCa                                                                                      >  2:30607/1‑149 (MQ=255)
                                                       ccTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCAt                                                                                     >  6:307006/1‑149 (MQ=255)
                                                         tgctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATacccac                                                                                                         >  1:44960/1‑127 (MQ=255)
                                                         tgctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATacccac                                                                                                         <  2:44960/127‑1 (MQ=255)
                                                              tACAAGGACGTGAACCGTCACGTACACCTGACGCTCACCCCCAAGCGCTTAGCCCAGGCTGATGCCCGCAAGCCCGCCACCCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAAcc                                                                                <  8:311421/147‑1 (MQ=255)
                                                                 aaGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCttt                                                                                                                                      >  4:352620/1‑90 (MQ=255)
                                                                 aaGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCttt                                                                                                                                      <  3:352620/90‑1 (MQ=255)
                                                                 aaGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgc                                                                           >  4:323928/1‑149 (MQ=255)
                                                                  aGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcc                                                                          >  3:198419/1‑149 (MQ=255)
                                                                   ggACTTGAACCGTCACGTAAAACTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcct                                                                         <  5:378535/149‑1 (MQ=255)
                                                                        ttcagcGTCAGGTGTACGTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCtgctgc                                                                    >  3:171546/6‑149 (MQ=255)
                                                                             cGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCt                                                               <  5:28441/149‑1 (MQ=255)
                                                                                   gTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTtg                                                         >  7:206297/1‑149 (MQ=255)
                                                                                        ccTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTgctgct                                                    >  4:320626/1‑149 (MQ=255)
                                                                                          tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                  >  3:85963/1‑149 (MQ=255)
                                                                                          tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                  <  4:171546/149‑1 (MQ=255)
                                                                                            aCGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                  <  5:355904/147‑1 (MQ=255)
                                                                                            aCGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                  <  4:96931/147‑1 (MQ=255)
                                                                                                    ccccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGc                                        >  7:26328/1‑149 (MQ=255)
                                                                                                        aaGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTg                                    >  1:158834/1‑149 (MQ=255)
                                                                                                        aaGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTg                                    >  5:236569/1‑149 (MQ=255)
                                                                                                          gcgcTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCttt                                                                                                                                      <  1:316280/49‑1 (MQ=255)
                                                                                                          gcgcTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCttt                                                                                                                                      >  2:316280/1‑49 (MQ=255)
                                                                                                           cgcTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGgcc                                 >  6:220263/1‑149 (MQ=255)
                                                                                                             cTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTGTGTTTGCCGCGCATACTGAAATACCCCCCCTTGCGGGGTACCGGCCATGAACCGCTGCATGCTGCTCGCTGGCTTGCTGCTGTGCACCTCTGGCCTGgccgc                               >  8:336922/1‑149 (MQ=255)
                                                                                                              tAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTgctgct                                                                   >  5:82306/1‑112 (MQ=255)
                                                                                                              tAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTgctgct                                                                   <  6:82306/112‑1 (MQ=255)
                                                                                                              tAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTCCTCGCTGCCTTCCTGCTGTGCACCTCTGCCCTGgccgcc                              >  5:191703/1‑149 (MQ=255)
                                                                                                               aaGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGgccgcc                              >  6:265769/1‑148 (MQ=255)
                                                                                                                aGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGa                            <  6:197422/149‑1 (MQ=255)
                                                                                                                    cAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGAc                        >  5:297759/1‑149 (MQ=255)
                                                                                                                     aGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACt                       >  3:3954/1‑149 (MQ=255)
                                                                                                                       gCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTgg                     >  3:339287/1‑149 (MQ=255)
                                                                                                                           aTGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGcc                 >  2:118863/1‑149 (MQ=255)
                                                                                                                            tGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCt                <  8:253996/149‑1 (MQ=255)
                                                                                                                            aaacagaaaGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCt                >  5:120218/8‑149 (MQ=255)
                                                                                                                              cccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTgg              <  4:3954/149‑1 (MQ=255)
                                                                                                                               ccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGc             >  3:334769/1‑149 (MQ=255)
                                                                                                                               cagaaaGCCCGCCAGGTTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGc             >  8:132785/5‑149 (MQ=255)
                                                                                                                                 gCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCt           <  8:319867/149‑1 (MQ=255)
                                                                                                                                  cAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTa          <  3:320626/149‑1 (MQ=255)
                                                                                                                                     gCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCACCGATGACTGGCGCCTGGCCTACGa       >  2:9785/1‑149 (MQ=255)
                                                                                                                                     ccccGACAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGCTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTCCTGTGGACCTCTGCCCTGGCCGCCGATGACTGGCGCCTCGCCTACGa       >  7:112375/2‑149 (MQ=255)
                                                                                                                                        cGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGGCGCCGATGACTGGCGCCTGGCCTACGACcg    <  1:30607/149‑1 (MQ=255)
                                                                                                                                           cAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACcgcg  >  7:199161/1‑148 (MQ=255)
                                                                                                                                             |                                                                                                                                                 
TGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCG  >  NC_002947/4002541‑4002827

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