Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,727 C→G 8.5% 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 8.5% 361.5 / 17.5 115intergenic (+30/‑60)allE/PP_3531S‑ureidoglycine aminohydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (57/48);  new base G (5/5);  total (62/53)
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 116 total aligned reads displayed.

ATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGC  >  NC_002947/4002578‑4002870
                                                                                                                                                       |                                                                                                                                               
aTGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAAc                                                                                                                                                 >  1:2162812/1‑150 (MQ=255)
       tGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCgg                                                                                                                                          >  2:1932388/1‑150 (MQ=255)
        gcgcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCggg                                                                                                                                         >  2:2150375/1‑150 (MQ=255)
         cgcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGc                                                                                                                                        <  1:822985/150‑1 (MQ=255)
          gcgcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCt                                                                                                                                       <  2:1338693/150‑1 (MQ=255)
            gcCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCttt                                                                                                                                     <  2:135437/150‑1 (MQ=255)
             ccTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTc                                                                                                                                    <  2:1656536/150‑1 (MQ=255)
             ccTTCTGACCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGCCGGTGCAACCACGACGTGAAGCGCATCCGGATGCTCGCAAGCCCTCCAGGTTGGCGGGCttgc                                                                                                                                    >  2:2330244/1‑148 (MQ=255)
              cTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgcg                                                                                                                                   <  1:2330244/150‑4 (MQ=255)
               ttCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgcgg                                                                                                                                  <  2:1251057/150‑5 (MQ=255)
               ttCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTg                                                                                                                                  >  1:1528814/1‑150 (MQ=255)
                tCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGt                                                                                                                                 <  1:631383/150‑1 (MQ=255)
                 cTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGtt                                                                                                                                <  2:179764/150‑1 (MQ=255)
                 cTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAACGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGtt                                                                                                                                <  1:1345340/150‑1 (MQ=255)
                    cccGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATTCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGc                                                                                                                             >  1:21882/1‑150 (MQ=255)
                       gCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCAAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCcgc                                                                                                                          <  1:1159905/150‑1 (MQ=255)
                          aGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGCCACTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCa                                                                                                                       <  1:1967174/150‑1 (MQ=255)
                            ccTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGtt                                                                                                                                <  1:804211/139‑1 (MQ=255)
                            ccTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGAAAGCCCGCCAACCTGGCGGGCTTTCTGtt                                                                                                                                >  2:804211/1‑139 (MQ=255)
                                  aCGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGaaa                                                                                                               <  2:2162812/150‑1 (MQ=255)
                                    gCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATa                                                                                                             <  2:179800/150‑1 (MQ=255)
                                         cggcCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATacccac                                                                                                        <  2:2226168/150‑1 (MQ=255)
                                          ggcCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATacccacc                                                                                                       <  1:1088208/150‑1 (MQ=255)
                                            cccAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCt                                                                                                     >  2:2185812/1‑150 (MQ=255)
                                                 gCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAg                                                                                                <  1:681242/150‑1 (MQ=255)
                                                   aaCTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAggg                                                                                              >  1:558482/1‑150 (MQ=255)
                                                     cTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTa                                                                                            >  2:1414124/1‑150 (MQ=255)
                                                      tttAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTc                                                                                                                                    >  1:2085296/1‑109 (MQ=255)
                                                      tttAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTc                                                                                                                                    <  2:2085296/109‑1 (MQ=255)
                                                      tttAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTAc                                                                                           >  2:654209/1‑150 (MQ=255)
                                                           cTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGcc                                                                                      >  2:559730/1‑150 (MQ=255)
                                                           cTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGATAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCGTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGcc                                                                                      >  2:175891/1‑150 (MQ=255)
                                                             aCCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgcgg                                                                                                                                  >  2:180484/1‑100 (MQ=38)
                                                             aCCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgcgg                                                                                                                                  <  1:180484/104‑5 (MQ=38)
                                                              ccTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATg                                                                                   >  1:1391644/1‑150 (MQ=255)
                                                              ccTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATg                                                                                   >  2:117795/1‑150 (MQ=255)
                                                                tgctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgc                                                                                                                                    >  2:2128750/1‑97 (MQ=39)
                                                                 gctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGTTCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAAc                                                                                <  1:25099/150‑1 (MQ=255)
                                                                  ctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCt                                                                                                                                       <  1:2071003/94‑1 (MQ=255)
                                                                  ctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCt                                                                                                                                       >  2:2071003/1‑94 (MQ=255)
                                                                     tACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGct                                                                            >  2:1385831/1‑150 (MQ=255)
                                                                     tACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGAAAGCCCGCCAGGTTGGCGGGCttgcgg                                                                                                                                  <  2:1369219/96‑5 (MQ=37)
                                                                        aaGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAg                                                                                                <  1:945819/127‑1 (MQ=255)
                                                                           gACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctgcctgc                                                                      <  2:1391644/150‑1 (MQ=255)
                                                                              gTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGCAGCAGCACATCCCGCCAGCAGATGTTGCTGTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCtgctgc                                                                   >  1:1498988/1‑150 (MQ=255)
                                                                                   ccGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAg                                                                                                >  2:159864/1‑116 (MQ=255)
                                                                                   ccGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAg                                                                                                <  1:159864/116‑1 (MQ=255)
                                                                                    cGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTg                                                             <  1:887617/150‑1 (MQ=255)
                                                                                     gTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGc                                                            >  1:1492346/1‑150 (MQ=255)
                                                                                      tCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCa                                                                                     >  2:1722981/1‑124 (MQ=255)
                                                                                      tCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCa                                                                                     <  1:1722981/124‑1 (MQ=255)
                                                                                          gTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTtgc                                                       >  2:259872/1‑150 (MQ=255)
                                                                                          gTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCACGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCTTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTtgc                                                       >  1:1945152/1‑150 (MQ=255)
                                                                                              acCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGctg                                                                           <  1:683073/126‑1 (MQ=255)
                                                                                              acCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTgctgct                                                   >  1:1427131/1‑150 (MQ=255)
                                                                                                cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                 >  1:906784/1‑150 (MQ=255)
                                                                                                cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                 <  1:117795/150‑1 (MQ=255)
                                                                                                cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAACCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGACCATGAACCGCTGCCTGCTGCTCGCTCCCTTgctgctat                                                 >  2:1761852/1‑148 (MQ=255)
                                                                                                 tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGttt                                                                                                                               >  1:1914219/1‑71 (MQ=255)
                                                                                                 tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGttt                                                                                                                               <  2:1914219/71‑1 (MQ=255)
                                                                                                 tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCAt                                                                                    >  2:744188/1‑114 (MQ=255)
                                                                                                 tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCAt                                                                                    <  1:744188/114‑1 (MQ=255)
                                                                                                 tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGAGTCTCCTTGAATACCCACCCTTGCAGGGTACCGGCCATGAAACGCTGCCTGGTGCTTCCTTCCTGGCTTCtgtg                                                >  1:1133701/1‑150 (MQ=255)
                                                                                                  gACGCTGAACCCCACGACCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGc                                               <  2:558482/150‑1 (MQ=255)
                                                                                                  gACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGc                                               >  1:479582/1‑150 (MQ=255)
                                                                                                   aCGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAgg                                                                                               <  1:1899612/101‑1 (MQ=255)
                                                                                                     gCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCAcc                                            >  1:1898281/1‑150 (MQ=255)
                                                                                                      cTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACct                                           <  1:259872/150‑1 (MQ=255)
                                                                                                       tGAACCCCAAGCTCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCt                                                                                                                                   >  2:1899974/1‑61 (MQ=255)
                                                                                                       tGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCt                                                                                                                                   <  1:1899974/61‑1 (MQ=255)
                                                                                                       tGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCAGTCAACCTAGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACctc                                          >  2:1892412/1‑150 (MQ=255)
                                                                                                             cccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGACAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGccc                                     >  1:1947051/1‑150 (MQ=255)
                                                                                                              ccAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCt                                    >  2:1420595/1‑150 (MQ=255)
                                                                                                                aaGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTgg                                  <  1:1385831/150‑1 (MQ=255)
                                                                                                                    gcTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGgccgc                              <  2:1427131/150‑1 (MQ=255)
                                                                                                                    gcTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGgccgc                              >  2:1249563/1‑150 (MQ=255)
                                                                                                                     cTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGcc                                      >  1:698626/1‑141 (MQ=255)
                                                                                                                     cTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGcc                                      <  2:698626/141‑1 (MQ=255)
                                                                                                                      tAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGATGGTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGccgccg                            >  1:22668/1‑150 (MQ=255)
                                                                                                                        aGCCCAGGCTGATGCCCGCAAGCCCGCCACCCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGAt                          <  2:1498988/150‑1 (MQ=255)
                                                                                                                          cccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGa                        >  1:2139672/1‑150 (MQ=255)
                                                                                                                           ccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGAc                       <  2:1898281/150‑1 (MQ=255)
                                                                                                                           ccAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGAc                       >  2:1869923/1‑150 (MQ=255)
                                                                                                                            cAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCa                                                                                                 <  2:299537/75‑1 (MQ=255)
                                                                                                                            cAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCa                                                                                                 >  1:299537/1‑75 (MQ=255)
                                                                                                                             aGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTTTTTGCCGCGCATACt                                                                                                                   <  2:908514/56‑1 (MQ=255)
                                                                                                                             aGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTTTTTGCCGCGCATACt                                                                                                                   >  1:908514/1‑56 (MQ=255)
                                                                                                                              ggCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTgg                    >  1:2198123/1‑150 (MQ=255)
                                                                                                                                 tGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAg                                                                                                <  2:2203273/71‑1 (MQ=255)
                                                                                                                                 tGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAg                                                                                                >  1:2203273/1‑71 (MQ=255)
                                                                                                                                   aTGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCt               <  2:2139672/150‑1 (MQ=255)
                                                                                                                                    tGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCCGACTCCCCAGCCTACTGAAAGGCCCACCCTTGCAGGGAAACCGCTATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTg              <  1:541219/150‑1 (MQ=255)
                                                                                                                                       ccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGcc           >  2:1115001/1‑150 (MQ=255)
                                                                                                                                       ccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGcc           >  1:688472/1‑150 (MQ=255)
                                                                                                                                       ccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGcc           >  2:684297/1‑150 (MQ=255)
                                                                                                                                        agaaaGCCCGCCAGGTTGGCGGGCTTGCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCt          >  1:1321271/4‑150 (MQ=255)
                                                                                                                                              gCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGAc     >  2:2104130/1‑150 (MQ=255)
                                                                                                                                                ccGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAACTACCCACCCTTGCAGGGTACCGGACATGAACCGCTTCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACcg   >  2:847730/1‑150 (MQ=255)
                                                                                                                                                ccGCCAACCTGGCGGGCTCTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCCACGACcg   >  1:104371/1‑150 (MQ=255)
                                                                                                                                                 cgcgAACCTTGCGGGCTCTCGGTTTGCCCCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTCCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCCCCGGTGACTGCCGCCTGGCCTACGACagc  >  2:1757495/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                               
ATGTGGCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGC  >  NC_002947/4002578‑4002870

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