Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,007,677 C→T 5.1% intergenic (+31/‑59) 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,007,6770CT 5.1% 1088.4 / 27.6 335intergenic (+31/‑59)allE/PP_3531S‑ureidoglycine aminohydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (184/134);  new base T (10/7);  total (194/141)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

TGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGC  >  NC_002947/4007528‑4007825
                                                                                                                                                     |                                                                                                                                                    
  tgGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTg                                                                                                                                                    >  2:2902267/1‑150 (MQ=255)
      tGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCgg                                                                                                                                                >  1:1021616/1‑150 (MQ=255)
      tGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCgg                                                                                                                                                >  1:6563925/1‑150 (MQ=255)
       gcgcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCggg                                                                                                                                               <  2:4237580/150‑1 (MQ=255)
        cgcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGc                                                                                                                                              <  1:4508533/150‑1 (MQ=255)
         gcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCt                                                                                                                                             >  2:2510167/1‑150 (MQ=255)
          cgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCtt                                                                                                                                            <  1:378481/150‑1 (MQ=255)
           gcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttg                                                                                                                                           <  1:5190156/150‑2 (MQ=255)
            ccTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTc                                                                                                                                          <  2:149053/150‑1 (MQ=255)
             cTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCt                                                                                                                                         <  2:6869407/150‑1 (MQ=255)
             cTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCt                                                                                                                                         >  2:489263/1‑150 (MQ=255)
              ttCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTg                                                                                                                                        <  1:2902267/150‑1 (MQ=255)
               tCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGt                                                                                                                                       <  1:1370311/150‑1 (MQ=255)
                            cTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCATCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATAc                                                                                                                          <  2:1714420/150‑1 (MQ=255)
                                tACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGaa                                                                                                                      >  1:5673575/1‑150 (MQ=255)
                                tACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGaa                                                                                                                      >  1:3595524/1‑150 (MQ=255)
                                 aCGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGaaa                                                                                                                     <  1:3125822/150‑1 (MQ=255)
                                   gCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATa                                                                                                                   <  1:4669254/150‑1 (MQ=255)
                                   gCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATa                                                                                                                   <  1:1188005/150‑1 (MQ=255)
                                    ccGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATac                                                                                                                  <  2:2253153/150‑1 (MQ=255)
                                      ggcggcCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATaccc                                                                                                                >  2:6452936/1‑150 (MQ=255)
                                             cAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTg                                                                                                         <  1:79304/150‑1 (MQ=255)
                                               ggCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCa                                                                                                       <  1:6369579/150‑1 (MQ=255)
                                               ggCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCa                                                                                                       <  2:265196/150‑1 (MQ=255)
                                                 cAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAgg                                                                                                     <  2:1613247/150‑1 (MQ=255)
                                                 cAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAgg                                                                                                     <  1:5518907/150‑1 (MQ=255)
                                                 cAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAgg                                                                                                     <  1:1543799/150‑1 (MQ=255)
                                                  aaCTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAggg                                                                                                    >  1:4925432/1‑150 (MQ=255)
                                                     tttAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTAc                                                                                                 >  2:1252270/1‑150 (MQ=255)
                                                     tttAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTAc                                                                                                 >  1:4790209/1‑150 (MQ=255)
                                                       tAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCg                                                                                               >  2:1368277/1‑150 (MQ=255)
                                                           tACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCa                                                                                           <  1:1957533/150‑1 (MQ=255)
                                                           tACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCa                                                                                           >  1:2355693/1‑150 (MQ=255)
                                                            aCCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCAt                                                                                          >  1:5890362/1‑150 (MQ=255)
                                                            aCCTGCTGCACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCAt                                                                                          >  2:3458081/1‑150 (MQ=255)
                                                             ccTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATg                                                                                         >  2:1860750/1‑150 (MQ=255)
                                                              ctgctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGa                                                                                        <  1:5542607/150‑1 (MQ=255)
                                                               tgctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGaa                                                                                       >  2:6267190/1‑150 (MQ=255)
                                                                gctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAAc                                                                                      >  1:111307/1‑150 (MQ=255)
                                                                       aaGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcc                                                                               <  1:698661/150‑1 (MQ=255)
                                                                          gACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcctgc                                                                            <  1:2987630/150‑1 (MQ=255)
                                                                          gACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTGCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcctgc                                                                            >  2:5496612/1‑150 (MQ=255)
                                                                                   cGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTg                                                                   >  1:1640696/1‑150 (MQ=255)
                                                                                   cGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTg                                                                   <  1:2024055/150‑1 (MQ=255)
                                                                                     tCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGcc                                                                 >  2:6658303/1‑150 (MQ=255)
                                                                                     tCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGcc                                                                 >  1:285063/1‑150 (MQ=255)
                                                                                       aCGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCtt                                                               >  2:2840989/1‑150 (MQ=255)
                                                                                       aCGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCtt                                                               <  2:3231335/150‑1 (MQ=255)
                                                                                       aCGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCtt                                                               <  2:6106065/150‑1 (MQ=255)
                                                                                               cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                       >  2:263912/1‑150 (MQ=255)
                                                                                                tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTAc                                                                                                 <  1:147398/107‑1 (MQ=255)
                                                                                                tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgtg                                                      >  1:6478522/1‑150 (MQ=255)
                                                                                                    gCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCAcc                                                  >  1:5799513/1‑150 (MQ=255)
                                                                                                    gCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCAcc                                                  >  2:3705266/1‑150 (MQ=255)
                                                                                                    gCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCAcc                                                  >  2:4841975/1‑150 (MQ=255)
                                                                                                      tGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGTCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACctc                                                >  2:505407/1‑150 (MQ=255)
                                                                                                       gAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACctct                                               <  1:3458081/150‑1 (MQ=255)
                                                                                                       gAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACctct                                               <  2:524336/150‑1 (MQ=255)
                                                                                                       gAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACctct                                               >  2:235607/1‑150 (MQ=255)
                                                                                                         aCCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGc                                             <  2:5417045/150‑1 (MQ=255)
                                                                                                           cccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGccc                                           >  2:4119885/1‑150 (MQ=255)
                                                                                                           cccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGccc                                           >  1:540115/1‑150 (MQ=255)
                                                                                                               aGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGgc                                       <  2:5855494/150‑1 (MQ=255)
                                                                                                                  gcTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGgccgc                                    <  2:1021616/150‑1 (MQ=255)
                                                                                                                   cTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGgccgcc                                   >  2:80156/1‑150 (MQ=255)
                                                                                                                    tAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGccgccg                                  >  2:1684244/1‑150 (MQ=255)
                                                                                                                     aaGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGa                                 <  1:6729702/150‑1 (MQ=255)
                                                                                                                      aGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGAt                                <  1:4631805/150‑1 (MQ=255)
                                                                                                                       gCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATg                               >  1:4199246/1‑150 (MQ=255)
                                                                                                                        cccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGa                              >  2:3775391/1‑150 (MQ=255)
                                                                                                                         ccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGAc                             >  1:3292825/1‑150 (MQ=255)
                                                                                                                         ccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGAc                             >  1:3832742/1‑150 (MQ=255)
                                                                                                                          cAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACt                            >  1:2272410/1‑150 (MQ=255)
                                                                                                                            ggCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTgg                          <  1:2554597/150‑1 (MQ=255)
                                                                                                                               tGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGgcgc                       >  2:4654096/1‑150 (MQ=255)
                                                                                                                               tGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGgcgc                       >  1:6150137/1‑150 (MQ=255)
                                                                                                                                  tGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTTTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTg                    <  1:1055590/150‑1 (MQ=255)
                                                                                                                                  tGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGc                                                                  >  1:3155102/1‑104 (MQ=255)
                                                                                                                                  tGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGc                                                                  <  2:3154763/104‑1 (MQ=255)
                                                                                                                                  tGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTg                    <  2:1579422/150‑1 (MQ=255)
                                                                                                                                    cccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGc                  >  1:2146317/1‑150 (MQ=255)
                                                                                                                                     ccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGcc                 <  1:180914/150‑1 (MQ=255)
                                                                                                                                      cGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGccc                                           >  1:4013447/1‑123 (MQ=255)
                                                                                                                                       gCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGccgc               >  2:3906270/1‑148 (MQ=255)
                                                                                                                                       gCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTa               >  2:1642555/1‑150 (MQ=255)
                                                                                                                                       gCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTa               >  1:4152898/1‑150 (MQ=255)
                                                                                                                                       gCAAGCCCGCCAACCTGGCGGGCTTGCTGGTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGGACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCGGTGCACCTCTGCCCTGGCCGCCGATGACAGGCGCCTGGCCtc               >  2:4580208/1‑149 (MQ=255)
                                                                                                                                        aaaGCCCGCCAGGTTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTAc              >  1:641989/2‑150 (MQ=255)
                                                                                                                                        aaaGCCCGCCAGGTAGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGATGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTAc              >  2:1172412/2‑150 (MQ=255)
                                                                                                                                         aaGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACg             >  1:965723/1‑150 (MQ=255)
                                                                                                                                          aGCCCGCCAGGTTGGTGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGa            >  2:4455423/1‑150 (MQ=255)
                                                                                                                                          aGCCCGCCAGGTTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGa            >  2:3399106/1‑150 (MQ=255)
                                                                                                                                          aGCCCGCCAGGTTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGa            >  1:2712331/1‑150 (MQ=255)
                                                                                                                                          aGCCCGCCAGGTTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGa            >  1:2613902/1‑150 (MQ=255)
                                                                                                                                           gCCCGCCAGGTTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGAc           >  1:6033900/1‑150 (MQ=255)
                                                                                                                                           gCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGAc           >  1:5054331/1‑150 (MQ=255)
                                                                                                                                             ccGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGa                                 <  2:2879011/126‑1 (MQ=255)
                                                                                                                                             ccGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACcg         >  1:1338141/1‑150 (MQ=255)
                                                                                                                                              cGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACcgc        >  2:2003167/1‑150 (MQ=255)
                                                                                                                                               gCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACcgcg       >  2:3593652/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGC  >  NC_002947/4007528‑4007825

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