Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,728 C→T 8.7% 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,002,7280CT 8.7% 324.4 / 15.7 115intergenic (+31/‑59)allE/PP_3531S‑ureidoglycine aminohydrolase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (56/48);  new base T (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 = 9.55e-01

Only 100 of 116 total aligned reads displayed.

GCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCA  >  NC_002947/4002583‑4002877
                                                                                                                                                   |                                                                                                                                                     
gCTGCGCGCCTTCTGCCCGCAACCCCGCTACGCCGGCGGCCCAGGGAACTTTAGCTACCTGCTGTACAAGTACGTGATCCGTCACGTACTCCTGACGATGATCCCCAAGCGCTAAGCCCAGGCCGATGCCCGCAAGCCCGCCAACCTGGc                                                                                                                                                   <  1:747344/150‑1 (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)
            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)
                       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)
                                                       tACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCa                                                                                            <  1:2368457/150‑1 (MQ=255)
                                                        aCCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgcgg                                                                                                                                         >  2:180484/1‑100 (MQ=38)
                                                        aCCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgcgg                                                                                                                                         <  1:180484/104‑5 (MQ=38)
                                                         ccTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATg                                                                                          >  2:117795/1‑150 (MQ=255)
                                                         ccTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATg                                                                                          >  1:1391644/1‑150 (MQ=255)
                                                           tgctgTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAGGTTGGCGGGCttgc                                                                                                                                           >  2:2128750/1‑97 (MQ=39)
                                                             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                                                                                                       <  1:159864/116‑1 (MQ=255)
                                                                              ccGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAg                                                                                                       >  2:159864/1‑116 (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:117795/150‑1 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                        >  2:1008797/1‑150 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCtgt                                                        >  1:906784/1‑150 (MQ=255)
                                                                                           cTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAACCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGACCATGAACCGCTGCCTGCTGCTCGCTCCCTTgctgctat                                                        >  2:1761852/1‑148 (MQ=255)
                                                                                            tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGttt                                                                                                                                      <  2:1914219/71‑1 (MQ=255)
                                                                                            tGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGttt                                                                                                                                      >  1:1914219/1‑71 (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)
                                                                                                      aCCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGc                                              <  1:627585/150‑1 (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                                             <  2:698626/141‑1 (MQ=255)
                                                                                                                cTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGcc                                             >  1:698626/1‑141 (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:684297/1‑150 (MQ=255)
                                                                                                                                  ccGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGcc                  >  1:688472/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)
                                                                                                                                                   cTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTAAGACCGCGAAGGCa  >  2:1339625/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GCTGCGCGCCTTCTGCCCGCAAGCCTGCTACGCCGGCGGCCCAGGCAACTTTAGCTACCTGCTGTACAAGGACGTGAACCGTCACGTACACCTGACGCTGAACCCCAAGCGCTAAGCCCAGGCTGATGCCCGCAAGCCCGCCAACCTGGCGGGCTTTCTGTTTGCCGCGCATACTGAAATACCCACCCTTGCAGGGTACCGGCCATGAACCGCTGCCTGCTGCTCGCTGCCTTGCTGCTGTGCACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCA  >  NC_002947/4002583‑4002877

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: